TRACKING ERROR BASED SLIDING MODE 
CONTROLLER DESIGN WITH APPLICATION 

TO FLIGHT CONTROL 


by 

GIRESH KUMAR SINGH 


ee/n«/p 

SiS^t 



DEPARTMENT OF ELECTRICAL ENGINEERING 

Indian Institute of Technology Kanpur 


FEBRUARY, 1998 





TRACKING ERfedR BASED SLIDING MODE 
CONTROLLER DESIGN WITH APPLICATION 

TO FLIGHT CONTROL 


A Thesis Submitted 

in Partial Fulfillment of the Requirements 
for the Degree of 
Doctor of Philosophy 


by 

GIRESH KUMAR SINGH 



DEPARTMENT OF ELECTRICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

February, 1998 



1 3 JUN ZOOO/Et- 

JENTRAL UURAB> 

i.i.7>, mmm 



To 

Mini, Baby, Papa, 

& 

Sweet Memories of Mummy 



CERTIFICATE 


It is certified that the work contained in this thesis entitled Tracking 
Error Based Sliding Mode Controller Design with Application to Flight 
Control by Giresh Kumar Singh has been carried out under my su- 
pervision and that this work has not been submitted elsewhere for a 
degree. 


Professor K. E. Hole 
Department of Electrical Engineering 
Indian Institute of Technology 
Kanpur 


February, 1998 


Acknowledgments 


I would like to take this opportunity to firstly thank all those who have taught me 
and from whom I have learnt. It is mainly because of the efforts of so many gurus 
that I have reached this stage. It is a great pleasure for me to express my deep sense 
of gratitude to Prof. K.E. Hole for his sincere and constant guidance throughout the 
present work. I have always found him very approachable, communicative, patient, 
and always willing to help in all aspects of life. His gentle proddings have been greatly 
responsible for bringing the present work to a logical conclusion. 

I am thankful to Profs. S.C.Raisinghani and C.V.R. Murthy for having introduced 
the fascinating area of aircraft dynamics and control and Dr. R.K. Bera for inspiring 
me to work on it. I am indebted to Profs. A.K. Raina, V.R. Sule, P.K.Kalra, and 
my guide for guiding me through the maze of control theory and for the valuable 
suggestions and advice they gave me from time to time. 

The time spent in the Control System Laboratory has been a very memorable one. 
It is only during my involvement in the setting up of laboratory experiments along with 
Dr. A.K. Raina, Shri Gurucharan Singh and Shri Uday Majumdar that I realized the 
difficulties in instrumentation and implementation of control systems. The company 
of Uday, Qadeer, Sashmita and Ramnath in the laboratory was always refreshing and 
encouraging. 

I would like to thank IIT Kanpur for the various facilities including the scholarship, 
the computer centre and the central library. 

\ 

I would like to thank my seniors Drs. IN. Kar, Deepak Murthy, and Shiv Narayan, 
and friends Ryali Venkat, Vydyanathan, Rajiv Kumar, Ashish Balaya and many others, 
who were always ready for a discussion or an argument on various technical aspects of 


IV 



V 


the work. These long talks have helped me in developing a better appreciation of the 
area. They have always been a source of inspiration. 

All work and no play! It is difficult to think of life without friends. During my stay 
in IIT Kanpur, I was fortunate enough to have many close friends, whose company will 
be remembered forever. The pleasure of being in the company of Swapan, Alok, Sujoy 
Fatta, Dipak, Raghu, Rahul, Abir, Him, Shivgovind, Khanna, Ajju, Bpandey, Kak, 
Kaki, Jauhari, Anjani, Shankha, Pankha, Medh, Sharma, Suman, Kishore, Mama, PK 
Dada, Bogada, Jawed, Immi, Kashi, Ksri, Palani, Govind, Mogambo, Tapo, Shanti, 
Mangsuli, Joydeep, Rogers, and others, can hardly be described in words - it is best 
experienced. I have always found them by my side when need be. The canteens of the 
various halls have been quite instrumental in nurturing our friendships, and their role 
is gratefully acknowledged. 

I can hardly ever forget the staff, hall4 and the faculty cricket teams and the matches 
between the three, and the discussions following them. These have gone a long way in 
making my stay here pleasant and memorable. I would like to thank the members of 
the family of Shri Uday Majumdar for extending warm hospitality during my stay in 
the campus. 

Closure: Last, but not the least, it is impossible for me to express in words my deep 
feelings towards my family. Looking back, I find myself lacking in the roles of a son 
as well as a brother. It is my parents’ and sisters’ whole-hearted love, support, and 
encouragement which has kept me going. Also the affection and support extended by 
friends and relatives back home during some of the turbulent periods of life has been 
very helpful. It saddens me to think of how my mother would have reacted on the 
completion of this thesis. To her, papa and my sisters, I dedicate this work. 

G.K. Singh 



Synopsis 


The Flight Control System (FCS) design for modern fighter aircrafts is one of the major 
problems encountered in its development phase. This design problem has assumed 
more importance with unstable airframes, ever expanding flight envelopes and involved 
technological developments. These developments are directed to enhance the combat 
agility and survivability. In fact, the development of new tactical maneuvers and 
strategies now hinges on the capability of the FCS. Conventionally, gain-scheduling has 
been the most widely used method to design the FCS. However, the method is rather 
ad hoc in nature and requires extensive simulation before it can be applied in flight. 
Recent developments in nonlinear control theory, namely the feedback linearization or 
the dynamic inversion, have also been applied to FCS design. However, guaranteeing 
stability and performance robustness in the presence of uncertainty and disturbances 
is not very straightforward and has been tackled using the approaches developed in the 
linear control theory. 

The present work is a study of the feasibility and applicability of the Slid ing Mode 
Control (SMC) method for the FCS design. The SMC method uses a nonlinear discon- 
tinuous control law and offers attractive robustness properties. The method originated 
in the Soviet literature| The method involves a) the selection of a sliding surface (SS), 
and b) the choice of a control law to make the SS attractive. The evolution of trajecto- 
ries may be considered in two distinct phases, a reaching phase, in which a trajectory 
reaches the SS, and a sliding phase where it will evolve according to the dynamics 
specified by the selected SS. Invariance to unknown but bounded matched uncertain- 
ties and disturbances in the sliding phase has been established. This is made possible 
by using a high gain or switching type of control law. However, practical applications 
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of the method are limited by the facts that 

1. It is directly applicable to input afSne systems only. 

2. Robustness in the presence of unmatched uncertainties is not guaranteed. 

3 The inputs to any practical system are the outputs of a dynamic actuation system 
with a finite bandwidth, and thus instantaneous changes in the input cannot be 
achieved ruling out the use of a discontinuous, term as in the conventional SMC 
method. 

These issues are addressed in this work for the aircraft longitudinal FCS design problem. 

The Chapter 2 considers the design of a constant plus proportional rate reaching law 
SMC for second order nonlinear systems. Conventionally, the controller parameters are 
selected so as to minimize the reaching time, i.e the time taken to arrive at the SS, and 
use the guaranteed invariance property offered by the SMC. However, a quick tracking 
performance in terms of the time taken to reach a small output tracking error (OTE) 
is of primary interest and should govern the parameter selection. For this, firstly it is 
established that this composite rate reaching law has a robust performance even before 
the onset of sliding. A new method of selection of the controller parameters based 
on the criterion of minimizing the time to reach a pre-specified OTE is proposed. It 
is established that controllers designed using the criterion of the time to OTE, while 
using similar peak levels of control effort, are expected to take substantially shorter 
times to the OTE as compared to designs based on the reaching time criterion. The 
method is applied to two second order mechanical systems, and the simulation results 
show that control laws designed based on the new criterion result in appreciably faster 
output tracking, verifying the parameter selection approach. 

Although the SMC method has excellent robustness properties, its applications 
to the FCS design problem have been few. The forces and moments generated are 
known to be complicated nonlinear functions of the aerodynamics and control surface 
deflections. Most of the published works use linearized approximations of the nonlinear 
aerodynamics to develop input-affine models for the aircraft dynamics, as the SMC is 
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directly applicable to this class of systems. In this work, the forces and moments are 
directly treated as inputs to the aircraft instead of the conventional control surface 
deflections, with the assumption that appropriate methods of generation of the desired 
forces and moments exist. A bounded multiplicative uncertainty is assumed between 
the command force vector and the actually realized force vector to account for any errors 
that might be involved in generating the desired force-moment vector. An explicit 
model following (MF) structure is used to design a SMC (MF-SMC) based full-envelope 
control law, using full state feedback, for the nonlinear longitudinal dynamics of an 
aircraft for the case of matched uncertainties and disturbances in Chapter 3. Analytical 
control laws for two operational modes of the aircraft, namely i) the velocity tracking 
mode, and ii) the trimming mode are derived. Results of the numerical simulation 
using parameters of the F-16 aircraft conflrm the robustness of the control law and 
establish the applicability of the design method. 

The SMC method is next extended to include the presence of unmatched uncer- 
tainties and disturbances as it would be present in most practical systems. The SMC 
design for this case is developed based on a multi-loop (ML) strategy. With the veloc- 
ities as primary outputs of the aircraft system, the uncertainties/disturbances present 
in the dynamics of the position state variables would be unmatched. It is assumed that 
model responses for the velocity and position state variables are available. The outer 
loop control law is designed so as to track the position state variables byl treating the 
command signals for the velocity variables as the input. The ML strate^ passes on 
the eSect of the unmatched disturbances in the outer loop to the inner loop. A SMC 
based inner loop control law is then designed to track the command velocity variables. 
The difference between the model and the command velocity variables has been de- 
rived and is shown to depend on the magnitude of the unmatched uncertainties. The 
ultimate boundedness of the position and velocity errors is established and their time 
domain bounds are obtained as functions of the controller parameters. These bounds 
can thus serve as a guideline for their selection. The method is applied to the design of 
the velocity tracking mode controller for the aircraft longitudinal dynamics. Numeri- 
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cal results for various choices of the controller parameters are presented in Chapter 4. 
The ML approach is then applied to single input (SI) higher order nonlinear systems 
with unmatched uncertainties. Analysis for a third order system shows that, with the 
proposed method of selection of the switching term gain, the conventional SMC design 
can be applied to tackle unmatched uncertainty. The errors in the various loops are 
shown to be ultimately bounded. Numerical simulation on a nonlinear system confirms 
the developed bounds and verify the applicability of the ML strategy. 

The chattering phenomenon, due to the use of a discontinuous control law, is one 
of the most serious limitations of the SMC approach which restricts its application to 
practical problems. Moreover, the control inputs to any system are generally applied 
through an actuation system with a finite bandwidth, ruling out the possibility of 
instantaneous switching capability. In this work, the input is considered to be the 
output of a specified first order system with disturbances and uncertainties. Available 
studies use higher order sliding modes and a differential algebraic approach to address 
this problem. However, both of these methods involve higher order derivatives of the 
states and thus the propagation of disturbances and uncertainties becomes complicated 
making robustness analysis very difficult. The ML strategy, developed in Chapter 4, 
is extended to the SMC design for such systems in Chapter 5. The propagation of 
disturbances in the ML strategy is better handled and obviates the need of derivatives 
of the uncertainties/disturbances leading to a simpler robustness analysis. The state 
and input errors are shown to be ultimately bounded. The method is first developed for 
a SI nonlinear system and applied to the nonlinear spring mass problem. Simulation 
results verify the approach. The method is then applied to the MI aircraft velocity 
tracking problem, considered earlier, with specified first order input dynamics, and 
the ultimate boundedness of the the velocity, position, and input errors is established. 
A parametric study of the choice of controller parameters is done. Results of the 
numerical simulation in the velocity tracking mode establish the applicability of the 
ML approach to this problem. 
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Chapter 1 
Introduction 


1.1 The Aircraft Flight Control Problem 

The aircraft flight control system (FCS) design is amongst the most critical problems 
faced in an aircraft development project. This design problem haa become more critical 
with unstable nominal aircraft designs, ever expanding flight envelopes, and ongoing 
technological developments. McRuer and Graham [1] chart out the historical develop- 
ments in the controller design since the flrst flight in 1901 by the Wright Brothers. The 
onset of high speed digital computers revolutionized the conventional control schemes 
resulting in the fly-by- wire and the active control technology [2, 3]. Novel concepts 
such as the Control Conflgured Vehicles (CCVs) [2, 4, 5] have been tried out on many 
operational aircrafts with associated beneflts. Another such emerging technology area 
is the thrust vectoring (TV) which will be equipping almost all the future generation 
fighter aircrafts [6-9]. This makes use of the fact that the engine thrust remains al- 
most constant as compared to the flow dependent conventional flight control (CFG) 
surfaces, namely the elevator, aileron and the rudder, and hence can be used even after 
the CFG surfaces are ineffective. To increase combat effectiveness and survivability, 
fighters are being pushed not only to the limits of the conventional flight envelopes 
but also beyond to acquire a tactical edge over the contemporary fighters. So much 
so that the development of new combat tactics and strategies is now governed by the 
capability of the FGS. The Tougachev’s Gobra maneuver’ routinely displayed by the 
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Sukhoi fighters [9] is one of the so called super-maneuvers. Other examples include 
the ’Herbst Maneuver’ in which the aircraft rotates its velocity vector through 180° 
in a minimum possible time. The governing theme guiding these developments is that 
maneuverability of manned flight should be bounded by the limitations of the human 
pilot and not the airframe. 

The aircraft is conventionally modelled as a 6 degree of freedom (dof) rigid body [10] 
( 3 for the translational and 3 for rotational). This results in a highly coupled and 
nonlinear set of dynamic equations. The dynamic modes of the aircraft may be broadly 
classified into two categories, i) Rotational modes, and ii) Flight path modes. The 
responsiveness of an aircraft to maneuvering commands, i.e its agility, is determined 
by the speed of the rotational modes. The natural frequencies of these modes are 
sufficiently high so as to make it impossible for the pilot to control the aircraft if these 
are lightly damped or unstable. Control systems are needed to augment the natural 
modes to have a suitable damping and frequency characteristics. These are popularly 
called the Stability Augmentation Systems. In addition to these, there are the Control 
Augmentation Systems, which help in achieving a desired response to the the pilot 
control inputs as in tracking etc. Autopilots are also used to perform pilot relief 
functions such as attitude hold or for navigational purposes. These control systems are 
intended to perform the functions designated to them over the entire flight envelope. 
The uncertainties and disturbances are an unavoidable feature of any practical dynamic 
system. These arise either because of the inexact knowledge of the parameters in 
the system (structured uncertainty) or because of unmodeled dynamics (unstructured 
uncertainty). A practical control system design besides performing stabilization should 
also provide 

• robustness to uncertainty and disturbances, and 

• some desired response characteristics- popularly known as Flying /Handling qual- 
ities (HQ). 

This necessitates the design and use of a flight control system (FCS) to meet stability 
and performance specifications over an entire flight envelope. 
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Linearized aircraft dynamics has been conventionally used to design linear control 
laws [10]. One approach was to design a single linear controller for the entire flight 
envelope. Hence the variation in aircraft dynamics at different flight conditions was 
modeled as a nominal plant with an additional uncertainty term. Yedavalli [11] con- 
siders robust linear controller design under structured uncertainties. Houpis et.al. [12] 
using quantitative feedback techniques design a robust control law for the F-16 aircraft 
accounting for four different flight conditions and six failure modes, i.e. a total of 24 
possible flight configurations. However, any scheme using a single linear control law for 
the entire flight envelope would achieve stability robustness only at the cost of sacrifice 
in performance. Moreover, the nature of changes involved in the aircraft dynamics 
with varying flight conditions is not completely unknown and may be incorporated in 
the design procedure. 

Thus, a better control strategy was the gain-scheduled (GS) approach. In this 
method, linearized models of the aircraft, which would approximate the actual dynam- 
ics better in some local neighbourhood, at the specified operating points are obtained. 
Separate linear controllers, with desired stability and performance robustness in a local 
sense, are then designed for each linear model obtained earlier. These point designs 
are then ’stitched together’ using a gain schedule to obtain the overall control law. 
Scheduling is heuristically done on parameters which i) vary slowly so that the actual 
plant may be considered quasi-steady, and ii) which capture the plant nonlinearities. 
Gangsaas et.al. [13] report practical applications of the LQG based gain-scheduled 
control law design. Adams et.al. [14] developed a a full envelope nonlinear control law 
for the longitudinal and lateral modes. Blight et.al. [15] design GS control law using 
LQG and structured singular value techniques to design for robust point designs. Shue 
et.al. [16] use a mixed H 2 /H 00 method for the point controller designs and schedule 
the controllers over the mach numbers. The main advantage of GS is that the wealth 
of linear control methods, performance measures, design intuition and computational 
tools may be used on the multi- variable nonlinear systems. However, the choice of the 
number of operating points and the point itself is not very well laid out, and therefore 
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the GS law has to be verified by extensive nonlinear simulation before it can be ac- 
tually applied to an aircraft. Traditionally, this approach has been the most popular 
and widely used control law with both analog and digital versions of the method being 
used in various fighter aircrafts such as the Mirage 2000, and the F-16 among others. 
However, there have been notable failures also, for e.g. the JAS-39 and the YF-22 
crashes [17, 18]. Shamma and Athans [19] discuss the application of this approach to 
non-linear plants. Shamma and Athans [20] point out some pitfalls of the method. 

With the recent developments using differential geometric concepts in nonlinear 
control theory [21, 22], nonlinear inverse dynamics or feedback linearization based 
flight control design has been investigated by many researchers. The main advantage 
of the method was that it directly results in a global/ full-envelope control law. It also 
offers generality of application to different aircraft types. Smith and Meyer [23] used 
model inversion techniques to the aircraft control design problem treating the acceler- 
ations as pseudo-inputs. Lane and Stengel [24] used the method to decouple specific 
state variables that are of some particular interest to the pilots and used it to enhance 
stall prevention capabilities. They observed that ’Flight control systems based upon 
nonlinear inverse dynamics offer the potential of providing improved levels of safety 
and performance over conventional designs developed using linearizing assumptions’. 
Vukobratovic and Stojic [25] discuss the centralized and decentralized control schemes 
using the inverse aircraft model. Singh [26] and Romano and Singh [27] design con- 
trollers for large nonlinear maneuvers of the aircraft. Durham et.al. [28] use a nonlinear 
model following approach to airplane control. Sun and Clarke [29] study auxiliary con- 
trol devices, namely the flaps and spoilers for direct force control, for active control of 
the aircraft. Enns et.al. [17] apply the method to the F-18 HARV executing a Herbst 
maneuver. 

However, this approach does not lead to an easy robustness analysis. This issue 
was addressed by using a multi-loop strategy with the inner loop performing the lin- 
earization of the unperturbed plant, and an outer loop was designed to take care of 
any uncertainties in the actual system. Adams and Banda [30] following the method 
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developed in [31] applied the structured singular value analysis to incorporate robust- 
ness in this approach. In another related work Sparks et.al. [32] use a gain-scheduled 
inner loop to obtain dynamic inversion. Enns et.al. [17] use loop shaping techniques to 
obtain robustness properties. An alternative method which uses a nonlinear discontin- 
uous control law and offers attractive robustness properties, namely the sliding model 
control method, is next considered for this control design problem. 

1.2 The Sliding Mode Control Method 

The Sliding Mode Control (SMC) originated in the Soviet literature [33]. Some of 
the recent survey /review papers on this topic are Utkin [34], DeCarlo et.al. [35], 
and Hung et.al. [36]. They reflect the increasing interest of the control community in 
the method and report the developments and applications of the SMC to a myriad 
of control problems. The SMC design method involves two steps, i) the selection of 
a stable hyperplane in the state/error space on which motion should be restricted, 
called the sliding surface (SS), and ii) the synthesis of a control law which makes the 
selected surface attractive. In this method a trajectory, starting from a non-zero initial 
condition, evolves in two phases : 

• A Reaching phase, in which it reaches the SS, and 

• A Sliding phase, in which the trajectory on reaching the SS, remains there for all 
further times and thus evolves according to the dynamics specified by the SS. 

The sliding surface specifies the dynamics in the second phase. Trajectories from any 
starting point can be shown to arrive at the SS in a finite time [37]. It has been shown 
that the trajectories on the SS exhibit invariance to unknown but bounded matched 
uncertainties and disturbances [21, 38]. This is achieved by using a switching/relay 
type of control law. However, the sliding motion commences only after the trajectory 
reaches the SS. Thus, to reap the benefits of the invariant response on the SS, it has 
been considered imperative to reduce the reaching time, i.e. the time spent in the 
first phase, and this has been traditionally used to select the controller parameters. 
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However, for a tracking task the quickness of response in terms of the time required 
to reach a small output tracking error (OTE) is of paramount interest and should 
be considered for the parameter selection. For the class of second order nonlinear 
systems, it is shown that a constant plus proportional rate reaching law as proposed 
in [37] offers a robust performance in the reaching phase also. Using this result, a 
new method of selection of controller parameters based on the time to OTE as the 
criterion is proposed. Simulations on two second-order mechanical systems show that 
the proposed parameter selection method takes substantially lower times to reach the 
OTE, while using similar levels of control effort, as compared with the selection method 
based on the traditional criterion of the reaching time. 

The SMC method has been used to design controllers for various practical appli- 
cations [36] including the aircraft FCS design. However, the method primarily applies 
to input affine systems, i.e. systems linear in the input, and has been applied to such 
systems. The SMC needs to be extended to tackle the presence of unmatched un- 
certainties and disturbances. Moreover, inputs to any practical system are generally 
generated by a dynamic actuation system and hence a discontinuous control law is not 
readily implementable. These issues need to be properly addressed for the method to 
gain more acceptance among the practising engineers and find more practical applica- 
tions. The present work attempts to study these factors as related to an aircraft FCS 
design problem. 

1.3 Application of SMC to Flight Control Design 

The SMC, with its robustness properties, is an attractive alternate method for the 
aircraft flight control design problem. However, its application to the FCS design has 
been rather slow to pick up as is indicated by the method finding no application in the 
special issue of the International Journal of Control on aircraft flight control design [17, 
39]. Calise and Kramer [40] were among the pioneers to adopt the SMC method for 
designing the FCS. They considered the linearized dynamics of the AV-8A Harrier 
in a hover mode under matched uncertainties. Singh [41] designed control laws for 
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large simultaneous longitudinal and lateral maneuvers of aircraft with linearized aero- 
dynamics and matched uncertainty. Hedrick and Gopalswamy [42] reported the pitch 
axis control system design with nonlinear dynamics and linearized aerodynamics under 
matched uncertainty. They also considered various response criteria as developed by 
HQ engineers. Thukral and Innocenti [43] consider a linear model for the longitudinal 
problem. They consider the stabilization and tracking problems. Fossard [44] coupled 
the SMC with model following for the helicopter problem. 

The above cited works model the nonlinear aircraft at best as an input affine system 
as the SMC technique is easily applied to such systems. However, the aircraft system 
of equations is known to be nonlinear in the states as well as the inputs [10]. Also, the 
number of independent inputs limit the number of outputs that may be exactly tracked. 
Thus the use of CFG surfaces, which is basically intended to control the three moments, 
imposes a fundamental limit on the number of tracked outputs. If, however, indepen- 
dent control of all the accelerations or equivalently the complete force-moment vector 
is possible, then these may be treated as pseudo-inputs, and the assumption of input 
affineness remains valid. It also increases the number of trackable outputs. Smith and 
Meyer [23] considered the accelerations as pseudo-inputs to the aircraft system. How- 
ever, their design did not include robustness properties. Franklin [39] uses the forces 
as input for the low speed longitudinal landing control problem. Since, the generated 
forces depend on the bounded (with estimate of the bounds known) but uncertain time- 
varying inertia properties, the forces and moments have been treated as inputs to the 
aircraft system in this work, and is thus called the force control approach (FCA). To 
realize the forces as desired by the controller, appropriate control effectors are assumed 
to be present. With TV-FCS gaining popularity [6, 9, 45] independent control of the 
forces and moments affecting the motion is a feasible assumption. In fact flight tests on 
unmanned scaled fighter models have shown the practicality of a pure TV-FCS [46]. It 
is however observed that a proper mix of the CFC surfaces and TV leads to maximized 
safety agility and controllability levels. Presence of more alternatives, via the CFC 
surfaces, may be used to advantage by choosing control surface deflections that also 



perform the CCV functions such as maneuver load control [2, 4]. This inverse problem 
of control allocation is being actively pursued in various development projects [7]. The 
constrained control allocation problem has been considered by Durham [47, 48]. The 
control surface deflections thus selected, on application to the airframe, are expected 
to generate the desired forces and moments. However, in practice the actually gener- 
ated force vector would be different from its desired values, because of various reasons 
including the unsteady uncertain aerodynamics and the generally neglected aeroelastic 
effects. This discrepancy between the desired and the actually generated force vector 
has been modeled as a bounded multiplicative uncertainty in this study. 

The explicit model following (MF) approach has been used to design the FCS. It 
offers many advantages. Firstly it facilitates the incorporation of various flight modes 
for different operational requirements. The HQs of an aircraft may be easily optimized 
by varying the models in-flight. One of its biggest advantages is that simulators can be 
widely used to train pilots with obvious advantages. In this work, a MF-SMC based 
full-envelope controller design for the nonlinear longitudinal dynamics of an aircraft 
has been done under matched uncertainties as is considered in almost all of the above 
cited works. A practical system however would not readily satisfy this assumption, 
and this is considered next. 

1.3.1 The Unmatched Uncertainty Case 

As seen earlier, there is a limit on the number of outputs that can be exactly tracked. 
The control of the remaining states is however still desirable. This was attempted by 
making use of the fact that there is a significant difference in time scales between the 
fast and slow states of an aircraft [49, 50]. In this multi-loop (ML) approach, a dynamic 
inversion based control law was synthesized for the slow states treating the fast states 
as the the input for the outer loop. Then an inner loop control law is designed so that 
the fast states actually track the command signals generated in the outer loop. Snell 
et.al. [49] report that better responses are obtained with this ML strategy as compared 
to the conventional gain-scheduled control law. The robustness of the control law was 
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however not dealt with. An integral error feedback in the outer loop was introduced 
by Azam and Singh [50] to obtain robustness while using a similar control structure. 

As seen earlier, the SMC law guarantees invariance against the matched uncer- 
tainties and disturbances only, however, the presence of unmatched disturbances and 
uncertainties is not uncommon in practical systems. This would adversely affect the 
intended system performance. Design of control laws for the unmatched uncertainty 
case has been considered by many researchers [51-53] under various assumptions on the 
types of mismatch. In this work, the ML strategy is extended to include the effects of 
unknown but bounded unmatched uncertainties. The propagation of disturbances and 
uncertainties, from the outer to the inner loops, in this approach is worked out. Using 
the ML approach, a SMC law is developed for a general multi-input (MI) mechanical 
system under unmatched uncertainty and applied to the nonlinear longitudinal aircraft 
dynamics. The ultimate boundedness of the position and velocity tracking errors is 
established and the associated bounds are calculated as functions of the selected con- 
troller parameters. This can thus serve as a guideline for their selection. Application 
of the ML analysis, to a third order nonlinear system with unmatched uncertainties 
and disturbances, shows that conventional SMC design procedure, with a new method 
of selection of the gain of the switching term, may still be applied to tackle unmatched 
uncertainty. The ultimate boundedness of the errors in the various loops is established. 

1.3.2 SMC with Input Dynamics 

The switching control law as used by SMC is not practically feasible as it would re- 
quire an infinite switching frequency to achieve the ’ideal’ sliding. However, in practice 
any actuation system will have a finite bandwidth and delays resulting in ’real’ slid- 
ing [54, 55]. Thus the invariance property is in effect lost and motion within some 
boundary of the SS can be guaranteed. This leads to the high frequency chattering 
effects which may excite the high frequency unmodeled dynamics in the system. To 
circumvent this difficulty, researchers have used two approaches i) use of continuous 
approximations of the discontinuous control law [21, 56], and ii) application of higher 
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order sliding methods [54, 55]. The control law then results in a continuous control 
law with discontinuities in the input rates. Ramirez and coworkers [57-59] using a 
differential algebraic approach construct input dependent sliding surfaces, and obtain 
similar results. The higher order sliding methods require higher order derivatives of the 
states and thus of the disturbances and uncertainties. This complicates the analysis 
of robustness properties, as the propagation of disturbances and uncertainties through 
the time derivatives is not very straight forward [60]. 

In this work, the input is considered to be the output of a first order actuator 
system. The ML strategy as discussed earlier is extended to address this problem. 
The propagation of disturbances and thus the robustness issues are more simply dealt 
with in the developed procedure. The method is first developed for a single input 
(SI) nonlinear system and ultimate boundedness of the state and the input errors 
are established. It is then applied to the MI aircraft velocity tracking problem with 
unmatched uncertainties. The position, velocity and input errors in the various loops 
are shown to be ultimately bounded. A parametric study of the involved feedback 
gains is required for their proper selection. 

1.4 Organization of the Thesis 

In this thesis, the aircraft control design problem using the SMC approach has been 
considered. Chapter 2 considers the design of a constant plus proportional rate reach- 
ing law based SMC design for second order mechanical systems. A method of selec- 
tion of parameters for the SMC which results in lower times to an acceptable level of 
OTE is proposed. Simulation results obtained for a single input (SI) nonlinear spring- 
mass-damper and a two input two-link robot problems verify the method. Chapter 3 
introduces the aircraft longitudinal dynamics and the force control approach. A MF- 
SMC based full-envelope controller design for the nonlinear longitudinal dynamics of 
a fighter aircraft has been done under the assumption of matched uncertainty. The 
desired model responses are obtained from a linear aircraft model given in Appendix 
D. The aircraft parameters, the assumed disturbances and uncertainties used in the 
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simulation are also given there. Analytical control laws for two operational modes, 
namely i) the Velocity Tracking Mode, and ii) the Trimming Mode, are derived for 
the case of matched uncertainty. A ML strategy is used to develop a SMC design 
procedure in the case of unmatched uncertainty in Chapter 4. Results of numerical 
simulation for the aircraft in the velocity tracking mode ascertain the robustness prop- 
erties. This approach is then applied to SI higher order systems. Synthesis of a SMC 
law with first order input dynamics is discussed in Chapter 5. The method is first 
developed for the SI nonlinear spring-mass problem. Simulation results are presented 
to verify the feasibility of the approach. The design procedure is then applied to the 
MI aircraft problem. Simulation results for the velocity tracking mode i;^ith unmatched 
uncertainties are presented. Conclusions and scope of further research form Chapter 6. 



Chapter 2 


Design of Proportional plus 
Constant Rate Reaching Law SMC 
for Second-order Systems 

2.1 Introduction 

Many mechanical systems are modeled as second order systems e.g. as in the trajec- 
tory following robots, or the angular position control in aircrafts as nose pointing or 
pitch hold autopilots. Other examples include the classical inverted pendulum, or the 
antenna pointing of satellites among others. The SMC method has been used for many 
such applications as is evident from [36, 37]. Thus it is of particular interest to study 
the effect of the SMC parameters for this class of systems. As discussed earlier the 
sliding motion, and thus the invariant performance, occurs only after the trajectory 
reaches the SS, while the tracking behaviour in the reaching phase depends on the 
uncertainties. Thus, a smaller reaching time would imply that sliding begins earlier 
and this has been conventionally used to select the SMC parameters. 

Gao and Hung [37] proposed a reaching law method which specifies the error dynam- 
ics in the reaching phase. They suggested three structures- the constant rate reaching, 
constant plus proportional rate, and the power rate reaching. The time taken to reach 
the SS is also evaluated. By appropriately selecting parameters the dynamic quality 
of the SMC system can be controlled. Choi et al. [61] suggest the use of a stepwise 
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time- varying SS, in which the SS is shifted and/or rotated as a function of time. They 
applied this method to a non-linear spring mass damper system and reported a hand- 
some decrease in time to achieve zero tracking error. However, the resulting controller 
is cumbersome to implement because the parameters of the SS need to be calculated 
at different time instants. Bartoszewicz [62] proposed a continuously time-varying SS 
that gives still faster tracking. The choice of the time varying SS is extended in [63] so 
as to minimize the integral of the absolute error. However, the parameter selection is 
dependent on the initial error conditions. 

In this chapter a study on the design of constant plus proportional rate reaching 
law for second order systems is done. Using Lyapunov theory, it is shown that the 
constant plus proportional rate reaching law imposes a second-order dynamics on the 
error variables, and thus guarantees a robust performance in the reaching phase also. 
The selection of controller parameters based on minimizing the reaching time is shown 
to be inadequate with regard to the time to reach a small output tracking error (OTE). 
A new method of the selection of the control law parameters based on the criterion 
of time to a specified OTE is proposed. To illustrate the method, it is applied to a 
nonlinear spring mass problem taken from Choi et al. [61] and a two-input two link 
robotic arm trajectory following problem [37]. Simulation results show that, for similar 
peak levels of control effort, considerably faster tracking performances are obtained for 
control laws designed using the proposed method. 

2.2 The Reaching Law based SMC Design Method 

Consider a single input, second order input-affine non-linear system of the form 

X = f(x,x) + g{x,x)u + df (2.1) 

where u is the control input, and df represents the bounded uncertainties and external 
disturbances, with \df\ < D. Since the uncertainties and/or disturbances which are 
in the range space of the input are called matched, it is obvious that any uncertainty 
entering in the scalar system as given above is matched. The arguments of / and g are 
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dropped from now on for brevity. The uncertainty in the input distribution function g 
is expressed as ^ = ^(1 + g), where g represents its nominal value, and g is bounded, 
i.e., 1^1 < C? < 1.0. 

Defining the state variables as Xi = x and X 2 = x, (2.1) can be written in the state 
space form as 

\^2 J y f + gu + df J 

Let the trajectory to be followed (desired trajectory) be given as xid and X 2 d = iid- 
Then the error between the actual and the desired trajectories can be expressed as 
e = Xi — Xid and e = X 2 — X 2 d- In vector form, 



^ Xi~ Xid 

^ X2— X2d 


(2.3) 


Let y — (e, e)'^, where the superscript T denotes the transpose. Taking the time 
derivative of the above equation and using (2.2) results in the dynamics 


y = 


f + gu + df- Xu 


(2.4) 


The switching function s for second order systems is conventionally defined as a 
combination of the errors in (2.3) [21, 37] as 


s = e + ae, 


(2.5) 


where, a > 0 sets the dynamics in the sliding mode. The SS selected above should be 
made attractive so that trajectories starting from a non-zero value of s approach the 
SS, and then stay on s = 0 for all future time instants. For this, using e from (2.4), 
the time derivative s is obtained from (2.5) as 


s — e + ae, or 

= f ■\-gu + df-Xid-^Oie, 

= f g{l-¥g)u-\- df - xu + oce. 


(2.6) 
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The input u remains to be chosen so as to make s go to zero for sliding. The normal 
procedure for its selection [21] is to express the input u as a sum of two terms. The 
first term is chosen, using the nominal plant parameters, so as to make s = 0, when 
s = 0, and is called the equivalent control. The second term is then chosen to tackle 
the uncertainties in the system. Note that in (2.6) g is the nominal g. Thus, the first 
term of u is chosen to cancel out the known terms in the right hand side of (2.6) and 
is seen to be g~^{xu — f — ae). The second term is selected as g~^{—es — iirsgn[s]) to 
introduce stabilizing terms in the s dynamics (2.6). The complete u thus becomes 

u = g~'^{xici- f - ae- es - Ksgn.[s]), (2.7) 

where, e,K>0 are positive real numbers to be selected, and sgn[-] is the signum 
function. It is established later on in this section that the above control law actually 
imposes a constant plus proportional rate reaching law [37]. Note that the u in the 
above equation may be simplified by using the definition of s in (2.5) and defining a 
and b as 

a = ae, and 

h = or + e, (2-8) 

both of which are obviously positive real numbers. Thus u becomes ' 

u = — ae - 5e - iiT sgn[s]). (2.9) 

Using this u in (2.6) and substituting g = ^(1 + g) gives 

s = / + gg~'^ {xid — f — ae — be — es — K sgn[s]) + dj — xu-¥ ae 
= 9{S:id — f — ae — be) — (ae + be) — (1 + p)jK'sgn[s] + df + ae 
= -ae- (b- a)e+g{xid — f — ae — be) — (l + g)Ksga.[s] + df. (2.10) 

Let 

R = xu- f-ae- be. (2.11) 
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Using R as above, a and b as defined in (2.8), and s from (2.5), s in (2.6) may be 
simplified to get 


s = — es + — (1 + 5)i{' sgn[s] + d/, or (2-12) 

ss = —es^ + s{gR + df) — {l-\-g)K\s\. (2.13) 

The scalar K should be chosen so that s above satisfies a constant plus proportional 

rate reaching law [37] as 

s = — es - r; sgn[s], or (2.14) 

ss = — es^-r7|s|, (2-15) 


where, 77 is a positive real number to be selected later. For ss in (2.13) to satisfy the 
above equation, K should be chosen large enough so that — (l+^)ii[r|s| in (2.13) always 
dominates the s{gR + dj) and in the worst case, —(1 + g)K\s\ 4- s{gR + df) is at least 
— 77|s|. Thus K may be selected [21] as 

G\R\ + D + 'n 




l-G 


(2,16) 


For this choice of K, it may be easily seen that the —(1 + g)K\s\ + s{gR + df) term in 
(2.13) is always < — Tyjsj. The term may thus be equivalently written as — r^'ls] where 
r]' > 77. The ss in (2.13) thus becomes 


ss = —es^ — 77']s|, or (2-17) 

s = —es — 77'sgn[s]. (2-18) 


It may be easily established for the s dynamics above that s = 0 is the only stable 
equilibrium point. The control law (2.9) is thus seen to satisfy the above reaching 
condition, i.e. it makes the s = 0 plane an attractive surface. It can be easily shown 
that the s dynamics as above will have better decay characteristics (See Appendix 
A) than the reaching law specified in (2.14). Furthermore, estimates of the reaching 
time i.e. the time to reach the SS (s = 0) may be evaluated using (A.7). Note that 
77' as defined above can be shown to be bounded with a maximum possible value of 
t]m = r] + 2K. The effect of the controller parameters, namely the positive scalars a, e 
and 77 are next studied. 



17 


2.3 Selection of SMC Parameters 


In this section a method of selection of the parameters in u (2.9) is proposed. Using s 
from (2.5), and s = e + ae from (2.6), (2.18) may be rewritten as 

e = —aee — {a + e)e — risgD.[s], or (2.19) 

e = —ae — be — Ty'sgnfe + ae], (2.20) 

with a and b as defined in (2.8). 

With this, the y dynamics in (2.4) may be rewritten as 


y = Ay + q 


( 2 . 21 ) 


where. 


and, 



( 2 . 22 ) 


(2.23) 


With a and b defined as in (2.8), the A matrix above is seen to have its poles at — a 
and — e. These poles govern the. response of the linear homogeneous part (i.e. the 
A-dynamics). Conventionally, these are selected to achieve a desired error response 
and lower reaching times. Also, the choice of a in (2.5) gives the rate of error decay 
in the SS. Note that K specifies the magnitude of the discontinuity in the control law 
in (2.9). The value of K in (2.16) is seen to increase with an increase in 77. A higher 
value of 77 leads to smaller reaching times as discussed in Appendix A and has been 
shown to give rise to more chattering [37]. 

Next the effect of the a and e parameters on the tracking performance is studied. 
This is done using the Lyapunov approach. A candidate Lyapunov function V for the 
y dynamics as specified in (2.21) is selected as 


V{y)=y'^Py 


(2.24) 
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where, P is a positive definite symmetric matrix and is taken as 


P = 


+ P 
a 


a 

1 


(2.25) 


and, /? > 0 guarantees the positive definiteness of the P matrix. The time derivative 
of V {y) may be obtained as 

V = if^Py + y^Py. (2.26) 

Using (2.21) in the above equation gives on simplification 


V = y^[A^P + PA]y + 2/Pg. 


(2.27) 


Since A in (2.22) is a Hurwitz matrix, it satisfies the Lyapunov equation 


A^P + PA = -<5 


(2.28) 


where, the Q can be derived from the above equation using A and P as given in (2.22) 
and (2.25) as 


Q = 


2a^e 2Q'e — /? 
20:6 — P 2e 
Note that positive definiteness of Q is guaranteed if 


(2.29) 


p < 406. (2.30) 

From q and P as defined in (2.23) and (2.25) respectively, the y'^Pq is seen to be 
equal to — s? 7 'sgn[s]. Using this and Q as obtained above in (2.29), V in (2.27) may be 
simplified to 

V = -y'^Qy - 2s77'sgn[s], or 

= -y'^Qy - 2?7'|s|, (2-31) 

which is seen to be always negative definite. Thus, the stability of e dynamics in (2.21) 
or equivalently in (2.20) is guaranteed. To study the time decay performance of V in 
the above equation, the following result is required. 
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Claim 1 For the second order error dynamics as given in (2.19), a Lyapunov function 
may be selected as in (2.24). For the Q matrix in (2.29), as obtained from the Lyapunov 
equation (2.28), there exists a positive real number 5 > 0 such that 

Q = 5P + H, (2.32) 

where H is at least a positive semi-definite matrix. 


Proof: It can be easily shown by a constructive proof that for the system under 
consideration, there is always a range of values that <5 can take. The detailed proof is 
presented in the Appendix B. In particular, it is established that permissible range of 
values of 8 lies in an open set and is given by 5 G (0, min[2Q;, 2e]), where min[-, •] gives 
the minimum of its arguments. ■ 

Using the decomposition of <5 as in (2.32), the V in (2.31) may be rewritten as 

V = + - 2s?7'sgn[s], or 

= -5V - y'^Hy - 2t]'\s\. (2.33) 

From the developments in Appendix A, it follows that the V decay response of the 
above system will be better than the response of the system described by 

V = -6V. (2.34) 


It is interesting to find out the Af, matrix with the dynamics specified as 


y = Aby, 


(2.35) 


that satisfies (2.3-11 The Aj should thus satisfy the Lyapunov equation (2.28) with the 
P as selected in (2.25), i.e. AjP + PAb = —SP, for some feasible choice of 5 > 0 as in 
(2.34). For this, the general form of the Ab matrix is taken as 


Ab = 


0 (f 


(2.36) 
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where, a', 6', c' are to be determined. The Lyapunov equation for Ab matrix as given 
above is solved for these unknowns for a general value of 5. These may be easily 
obtained by the solution of a set of simultaneous linear equations in a', b', d as 

, _ <5(q:^ + /3) 

2a 

h' = 6 


It is known that (see Appendix B) 5 can take values as large as 2e or 2q;, in the limit 
as /? -> 0, depending on whichever is smaller. Using (2.37), the matrix for 5 = 2a 
may be obtained as 




0 1 
— (a^ + /?) —2a 


(2.38) 


and is seen to have both of its poles at —a, for /3 = 0. Similarly, for 5 = 2e, the A^ 


matrix becomes 


0 e/a 

-{ea + ep/a) —2e 


(2.39) 


and for /5 = 0, both the poles of the above Ab matrix are easily seen to lie at — e. It 
follows from Appendix A that the decay of the Lyapunov function (2.24) as obtained 
in (2.33), i.e. for the dynamics in (2.21) with the A matrix in (2.22) having its poles 
as —a and — e, is guaranteed to be better than the decay response in (2.34), which is 
achievable by a Ab matrix as selected in (2.38) or (2.39) which has both of its poles 
at — min[a, e]. This implies that the error dynamics in (2.21) in the reaching phase 
is bounded by the response of the second order Ab dynamics in (2.35) and a robust 
performance in the presence of the matched uncertainties and disturbances is achievable 
even before the onset of sliding. Estimates of the guaranteed times to reach the OTE 
may be obtained using the At dynamics. It is obvious that shorter times to the OTE 
are expected as both a and e are increased. However, the control law in (2.9) will 
also require larger control magnitudes as a and b will also increase as can be seen from 
(2.8). Next the effect of the relative values of these parameters is considered. 
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2.3.1 The Relative Choice of a and e 

Let the poles of the A matrix in (2.22) be selected as —p and — g, where p and q are 
positive real numbers satisfying p > q This could be achieved in the following two 
ways: 

• Case CP’, a = p and e = g. This would result in the SS being e+pe, say Sp. 

• Case CQ: a = q and e = p, with which the SS becomes e + ge, say Sq. 

It follows from (2.8) that a = pq and b = p + q remain the same in both the cases 
resulting in the same A matrix in (2.22). With o and b being the same, R in (2.11) 
and thus K (2.16) are expected to have comparable magnitudes for both the cases. 
Thus the peak values of the control effort u in (2.9) should also be nearly equal. The 
performance obtained in the two cases is compared next. For this, a common Lyapunov 
function is first constructed as follows. 

Claim 2 Let P and Q he symmetric positive definite matrices as specified in (2.25) 
and (2.29) respectively, for the Case CP, i.e. a — p and e = q. Then, there exists a 
sufficiently large positive scalar M such that 

V{y) = y'^Py + M\Sp\ + M\Sq\ (2.40) 

where, y = (e, effi , is a valid common Lyapunov function for Cases CP and CQ. 

Proof: For F to be a Lyapunov function, at least negative semi-definiteness 

(n.s.d.) of its time derivative must be established for the two cases. Detailed derivations 
are given in the Appendix C. Let K, i = p,q represent the time derivative V of the 
Lyapunov function in (2.40) for the cases denoted by the subscript ’i’. The Vp has been 
derived in (C.6) - (C.9) for the Case CP and Vq in (C.12) - (C.17) for the Case CQ 
respectively. The negative definiteness of the time derivatives in the two cases is shown 
for a choice of M as in (C.19). ■ 

The discussion in the Appendix C shows that (C.6) - (C.9) giving Vp for the Case 
CP may be combined to give the time derivative as in (C.IO). Similarly (C.20) gives 
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Vq for the Case CQ. Further, for the Case CQ, it is shown in the Appendix C that 
there exists a 0 < 9 (< p by selection), using which Vq in (C.20) may be rewritten 

as 

V, < -V^Qv - '5,Af[iS,| + li'pll - (p - (2-41) 

From Claim 1, it is known that the Q matrix may be decomposed as given in (2.32) 
as C = 5qP + Hq, where Hq is a positive definite matrix, since 5q < q < min[2g, 2p]. 
Using this decomposition, the Vq above may be expressed as 

V, < -sypy - y'^H,v - + |Sp|) - (p - or 

< -S,V,-{jp-5,)M\S,l (2.42) 

using the definition of V in (2.40) where the subscript ’g’ represents the Case CQ. As 
p > q > 5q, the (p — 5q)M\Sq\ term in the above equation is always non-negative, and 
thus the results of Appendix A imply that the Vq decay response of this system will be 
better than the exponential decay response offered by 

Vq = -6qVg, where 6q e (0,q). (2-43) 

Note that 6q can assume a maximum value oi5q = q in the limit. 

For the Case CP, the Vp is seen from (C.IO) to be given as 

V, < -y'^Qy - M[q\s,\ +p|5,|I - 27,'|Sp|. (2.44) 

The Q matrix in this case is decomposed as 5pP Hp, where 0 < 5p < 2q = min[2g, 2p] 
has to be found out. Following the stef)s as for the Case CQ, the Vp above may be 
rewritten as 

K < -W - y'^H,y - (p - i5p)iW|S,| - (? - S„)M\S,,\ - 2p'|S,|, or 

< -«-v’-J3pP-(p-«Af|S',|-(?-<5p-l-^)M|Sp|. (2.45) 

2r]' 2rf 

Thus, it is seen that if 5p < minlp, (P~<5p)M|5g| + (g— 6p + — )M|5pI 

term in the above equation is always guaranteed to be non-negative. As discussed in 
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the Appendix C, the condition on M for to be a Lyapunov function is M > 0. 
For M — )■ 0, the 2t)' jM term will increase, However, since the decomposition in (2.32) 
is possible only for 5p < 2q = min[2g, 2p], the 5p may be selected to be smaller than 
min[p, 2q]. It thus follows that the Vp decay response of the above system will be better 
than the response of the system given by 


Vp = —5pVp, where 5p G (0,min[p, 2g]). (2-46) 


Thus the upper limiting value of 5p is seen to be given by 5p = min[p, 2q]. Thus it is 
obvious that 5p > 5g. 

With the guaranteed decay responses of Vq and Vp being given as in (2.43) and 
(2.46) respectively, it is clear that Case CP guarantees a faster exponential decay of 
the Lyapunov function (2.40) as compared to the Case CQ. These guaranteed decay 
performances may be used to develop an estimate of the time to reach a specified level 
of OTE in terms of the Lyapunov function. Thus, given an initial value V'(O) at t = 0, 
the time tf required for Kj i =9, p dynamics, in (2.43) and (2.46) respectively, to decay 
to a final value of Vf may be evaluated as 



V(0) ' 

. V, , 


(2.47) 


where, i =9, p, and 4 are as defined earlier and represent the maximum limiting value 
of 5^ for the 5q and 5p as given in (2.43) and (2.46) respectively. For the same V{Q)/Vf, 
it is obvious that tf in the Case CP will be smaller than in the Case CQ, as the right 
hand side of the above equation will be smaller because of the fact that 5p > Sq. The 
tf in (2.47) gives conservative estimates for the decay of the Lyapunov function. The 
actual times tors required, for the V decay responses in (2.42) and (2.45) for the Cases 
CQ and CP respectively, to reach the Vf are expected (see Appendix A) to be smaller 
than the tf, i.e. tore < tf- Thus the a and e parameters should be selected as in the 
Case CP to have a smaller value of tf, i.e. for a given choice of the poles of the A 
matrix in (2.22), e should be selected as the magnitude of the dominant pole, and a as 
the magnitude of the other, or simply a> e. This is summarized as follows. 
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Lemma 1 The choice of a and e parameters, for the constant plus proporiional rate 
reaching SMC law as in (2.9), with a and b selected as in (2.8) and the sliding surface 
as in (2.5), such that 

a > €. (2.48) 

guarantees a faster decay of the Lyapunov function (2.40). 

Note that the above lemma indicates that the choice of a and e as in (2.48) only 
ensures that tf will be smaller. It however does not guarantee that the tors for the 
Case CP will be smaller than the Case CQ, although it is so expected and can actually 
be seen from the simulation results presented in the next section. There may however 
be parameter combinations where this would not hold. 

A smaller value of V implies a smaller tracking error in the sense of the norm (as 
defined by the Lyapunov function). Thus, since the Case CP has a faster guaranteed 
decay performance, it is expected to take lesser times as compared to the Case CQ to 
reach an acceptably small OTE. As discussed earlier, the control effort is expected to 
be almost the same for both the CP and CQ cases. Note that, lower reaching times 
correspond to the case of e > a (Case CQ) as can be seen from the results in the next 
section. However, the above theoretical developments show that shorter times to a 
desired OTE are expected for the cases where a > e (Case CP). The estimates of the 
time to the OTE using (2.47) showing the effect of the relative choices of a and e as 
suggested in Lemma 1 are verified by applying the design method to two mechanical 
systems. The systems and the obtained simulation results are presented next. 

2.4 Simulation Results 

For the illustration of the proposed SMC design method, the results presented in the 
previous section have been used to design controllers for two practical problems. The 
numerical results of these examples ascertain the validity of Lemma 1. 
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2.4.1 Nonlinear Spring Mass Damper System 

The first example concerns a non-linear mass-spring-damper system from Choi et. 
al. [61]. The governing differential equation is 

mx + u{x, t) -f- fi{x, t) — u{t) + d(t), (2.49) 

where, x,x represent the position and the velocity of the mass m, u{.) represents the 
applied input force, and d{.) the disturbance signal. i/(., .) and //(., .) represent the 
damping and spring force term respectively. The system parameters and functions 
are given to be m = 1, f = fi + f 2 + fs + fi, which are given as /i = yUoo;, /z = 
~fj,ix^, fs = —i/qx, and fi = —Uix\x\ with fj.o = 0.5, /ii = 0.5, uq = 0.3, and i/i = 
0.3. The system also has uncertainties expressed as ui/i + a 2 f 2 + 03/3 + ^4/4 with 
ai = —.05 + 0.25sin(57rt), 02 = -.05 -1- 0.25sin(57rt), 03 = —.05 -1- 0.15sin(77rt), and 
04 = — 05 -f 0.15sin(77rt) and disturbance d — 0.05 4- 0.25cos(37rt). The signal to be 
tracked is given by x^ = —0.5 cos(7rt/5). The above system in (2.49) may be expressed 
in standard form as in (2.1). The initial conditions are taken as 0.5 and 0 for the two 
states respectively. 

Simulation was done using the fourth order Runge-Kutta method with a time step 
of 2.5ms. The remaining parameters were selected as P = 1 and 77 = 1. Several 
combinations of a and e were used. 

• Caise 1: Q! = 2 and e = 6 

• Case 2: O' = 6 and e = 2 

• Case 3: O' = 3 and e = 4 

• Case 4: a = 4 and e = 3 

• Case 5: a = 4 and e = 5 


• Case 6: a = 5 and e = 4 



Table 2.1: Comparison of time for OTE and Reaching: Example 1. 


Case 

no. 

a 

e 

Time to OTE (< 1 x 10“^) in (s) 

Reaching time (s) 
(A.7) 

t, (s) (2.47) 
for 

Bound obtained 

from Ah dynamics 

(2.35) 

Actual 

1 

2 

6 

4.24 

3.30 

0.183 

4.605 

2 

6 

2 

4.24 

2.00 

1.472 

2.302 

3 

3 

4 

2.83 

2.35 

0.426 

3.070 

4 

4 

3 

2.83 

2.00 

0.732 

2.302 

5 

4 

5 

2.13 

1.80 

0.439 

2.302 

6 

5 

4 

2.13 

1.65 

0.650 

1.842 


The time to OTE (defined as the time after which the error is less than 1.0 x 10"^) 
obtained for the bounding ^6 dynamics in (2.35), the actual tracking times, the reaching 
times calculated from (A.7) and the tf for Vo/Vf = 1.0 x 10^ as evaluated from (2.47) 
for all the cases are presented in Table 2.1. It can be seen from the table that the 
actual time to achieve the OTE is appreciably smaller for Case 2 as compared with 
Case 1. The Table 2.1 also shows that similar results are obtained for Cases 4 and 6. 
In all these cases a > e, thus verifying Lemma 1. It is noticed from Table 2.1 that, 
for all of these cases, the reaching time is greater than the corresponding Cases 1, 3 
and 5. Although Cases 1,3 and 5 have smaller reaching times, they require about 65%, 
17.5%, and 9.1% more time to the OTE than the time required by the Cases 2,4, and 6 
respectively. Thus, from the table it is clearly evident that parameter selections giving 
lower reaching times actually take longer times to reach the specified OTE. Moreover, 
the time bounds tj calculated from (2.47) are seen to be better than those obtained 
from the dynamics (2.35) in the sense of being able to bring out the eflfect of the 
relative values of a and e. 

Choi et al. [61] reported a time of 2.58 seconds to achieve a zero tracking error 
(OTE). It can be observed from Table 2.1 that in almost all the cases under consid- 
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time (s) 


Fig. 2.1: Tracking responses for Case 1 (’Cl’) and Case 2 (’C2’). 

eration, a lesser tracking time is achieved by this approach. In fact the best tracking 
time achieved is 1.65 seconds for Case 6. In all of the cases the actual time to reach 
OTE is smaller than the bound given by the Ab dynamics and the tf. As qualitatively 
similar responses were obtained, the responses obtained for the parameter combination 
in Cases 1 and 2 are presented in Figs. 2.1, 2.2, and 2.3 respectively. The responses 
are indexed by the case numbers in these figures. Note that Case 1 corresponds to 
Case CQ, and Case 2 to Case CP. Fig. 2.1 shows the tracking response obtained 
in the two cases. It is evident from the figure that the ’C2’ (satisfying the Lemma 
1) trajectory approaches the OTE in appreciably lower time as compared to the ’Cl’ 
trajectory. Fig. 2.2 shows the s response in the two cases. The ’Cl’ response has a 
lower reaching time (see also Table 2.1) as compared to ’C2’. Fig. 2.3 shows the error 
decay with time for both the cases and also for the solution of the At dynamics in 
(2.35), labeled as e^. Although, only a better decay of the Lyapunov function in (2.24) 
is guaranteed for the control law in (2.9), the figure shows that the decay of e obtained 
in the two cases are better than the solution of the Ab dynamics in (2.35). 

The required control effort in the two cases are shown in Fig. 2.4. As discussed in 
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time (s) 


Fig. 2.4; The control effort u required for Cases 1 and 2. 

the previous section, since a = 12 and 6 = 8 from (2.8) are same for the cases ’Cl’ and 
’C2’, thus the control efforts in these cases should be of about the same magnitudes 
This can indeed be seen from Fig. 2.4, thus it verifies that that the proposed method 
uses about the same control effort to achieve the same tracking accuracy much faster. 
However as the reaching time in ’Cl’ is smaller, the input in Case ’Cl’ begins chattering 
earlier as can be seen from the figure. Thus the selection of the parameters can be easily 
done drawing on the experience of the initial condition responses of linear second order 
systems. If the a and e are chosen so as to satisfy Lemma 1, better results are expected. 

2.4.2 Two-link Robotic Arm 

The design method given in Section 2.3 is extended to multi-input case consisting of 
the two-link robot arm problem [37]. The equations governing the dynamics are given 
as 

H{q)q + C(q,q)q = r + d (2.50) 

where, q = (^ 1 ,^ 2 )^ gives the angles of the two links in radians. H{q) is the positive 
definite inertia matrix. The vector t = (r^Ti)*^ is the control torque vector, and 
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d = {di, d 2 )'^ represents the disturbance vector with bounds on their magnitude given 
as |di| < Di, and \d 2 \ < D 2 . The expressions for all the functions in the above equation 
are taken from [21]. The above equation can be rewritten as 

q = H-'^{q)[T + d-C{q,q)q]. (2.51) 


The trajectory to be followed is given as qd = and its time derivatives. The 

error vector is defined as before. 




^ 

ei 


9i — 6id 

62 


9% — 92d 

^62 y 


\ 92 — 92d ) 


Following the developments in Section 2.2, the input torques are selected as 


(2.52) 



H{q) 


+ C(g, q)q. 


(2.53) 


hd - Oiei - hiei - .Disgnfsi] ^ 

92d - 0262 - ^262 - 52 Sgn[s 2 ] j 

where, the quantities Oi, 5i, 02, 62 > 0 are selected by appropriate choices of ai, a 2 and 
Cl, €2- Also Di — Di + qi, and T>2 = T>2 + ^2- Under the application of the control law 
the system decouples into two single input systems. It is easily seen that the above 
control law satisfies the following constant plus proportional rate reaching law 


si < — eiSi — 77isgn[si], and 

h < -£252 - %sgn[s2]. (2.54) 


Then the choice of the parameters can be done based on the discussion in Section 2.3. 

Simulation was done as in the previous example, but with a time step of 500/iS. 
The parameters of the robot, the disturbance signals, and the controller parameters are 
taken from [37] in which the SMC parameters were taken as o: = 1, and e = 15 (Case 
CQ of Section 2.3). The Lemma 1 indicates that the choice of the two parameters 
should be the other way round. P and 77 were selected as in the earlier example. The 
Table 2.2 presents the various time bounds and the actual tracking times for four 
parameter combinations. 
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Table 2.2: Comparison of time for OTE and Reaching: Example 2. (Ll) and 
(L2) denote responses for the Links 1 and 2 respectively. 


Case 

no. 

a 

6 

Time to OTE (< 1 x 10 in (s) 

Reaching time (s) 
(A.7) 

tf (s) (2.47) 
for 

= io4 

Vf 

Bound obtained 

from At, dynamics 

(2.35) 

Actual 

1 

2 

6 

4.24 

3.128 (Ll) 
2.890 (L2) 

0.183 

4.605 

2 

6 

2 

4.24 

1.095 (Ll) 
0.910 (L2) 

1.472 

2.302 

3 

1 

15 

2.83 

5.683 (Ll) 
5.475 (L2) 

0.027 

9.210 

4 

15 

1 

2.83 

0.563 (Ll) 
0.423 (L2) 

4.732 

4.605 


It is seen from the table that the Cases 1 and 3, corresponding to the Case CQ, 
have very small reaching times as compared to the Cases 2 and 4 which correspond 
to the Case CP. Although the Case CP exhibits a much larger reaching time, it is 
expected from Lemma 1 and the results in the previous section that it should exhibit 
a faster tracking performance as compared to the Case CQ. The times to OTE for the 
Case CP is indeed seen to be substantially smaller than the corresponding Case CQ. 

Since the qualitative nature of the responses in the Cases 2 and 4 are similar, as 
was also observed for the the other two cases, the error responses obtained only for the 
Cases 3 (labeled as C3) and 4 (labeled C4) for the first link is presented in Fig. 2.5, 
and that for the second link in Fig. 2.6. As is indicated by the very small reaching time 
for the Case CQ, it was observed from the numerical simulation results that in Case 
CQ, sliding begins almost instantaneously as compared to the Case CP. However, 
from the CQ response in the figures, the time to OTE is seen to be about 5.5s as 
compared to about 0.5s for the Case CP. This again shows that substantially better 
performances are achievable by satisfying Lemma 1. It also substantiates the view that 
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Fig. 2.5; Error responses for the first link for Cases C3 (CQ) and C4 {CP). 

a lower reaching time does not imply faster tracking. Thus, it may be concluded that 
the time to OTE should be used as the criterion to select the controller parameters 
rather than the reaching time. 

The required control efforts T\ and rj (2.53) for parameter combinations in Cases 
3 and 4 are presented in Figs. 2.7 and 2.8 respectively. It is observed that required 
torque levels for both cases are similar. 

2.5 Conclusions 

In this chapter the design of the constant plus proportional rate reaching SMC for 
second order nonlinear systems has been studied. An investigation of the effect of the 
choice of the controller parameters has been done. Using a Lyapunov approach it has 
been shown that the design method imposes a bounding second order dynamics on the 
error variables, and thus a robust performance in the reaching phase is guaranteed. A 
new method of selection of the SMC parameters which results in smaller times to reach 
a given level of OTE, while using almost similar levels of maximum control effort, is 
proposed. It is seen that for the class of systems under study, a small reaching time does 







34 



0123456789 10 

time (s) 

Fig. 2.8: The input torque T 2 . 

not guarantee a fast tracking response. The method of parameter selection should thus 
be based on the time to OTE, rather than on the criterion of the reaching time. The 
design technique has been applied to a nonlinear spring-mass-damper system under 
matched uncertainties and exogenous disturbances. Prom the simulations done it is 
seen that the control laws based on the lower reaching times require approximately 
10 — 60% more time as compared to the proposed control laws to reach the same level 
of OTE, while using the same levels of control magnitudes. The method developed for 
a single input system has been extended to a two input two link robotic arm problem. 
Results of numerical simulation for the parameters used show that significantly faster 
tracking is obtained with the proposed control law, thus establishing the utility of the 
method. However, more work is required to develop better estimates of the time to 
OTE. 




Chapter 3 

SMC for the Longitudinal 
Dynamics of a Fighter Aircraft 

3.1 Introduction 

The aircraft flight control system (FCS) forms the backbone of the current generation 
fighter aircrafts. This chapter discusses the application of the SMC method for the 
design of FCS for fighter aircrafts. The aircraft is conventionally modelled as a 6 
degree of freedom (dof) rigid body [10]. As discussed in Chapter 1, the nonlinear 
inverse dynamics based flight control design has been applied in many works [17, 23, 
24, 27-29, 39] to various aircrafts for varying purposes. However, the issue of robustness 
of the feedback linearization techniques to disturbances and uncertainties still needs to 
be addressed. The SMC, with its robustness properties, offers an alternative method 
for this problem. However, it has not been readily used for the FCS design problem 
as is evident from the lack of a single apphcation paper in the the special issue of 
the International Journal of Control (Vol. 59 no.l, 1994) on aircraft flight control 
design [17, 39]. 

The SMC method has earlier been used for FCS design problem by some researchers. 
Calise and Kramer [40] applied the SMC design method to the linearized djmamics of 
the AV-8A Harrier in a hover mode under matched uncertainties and compared the 
obtained responses with those obtained using a LQ based design. They observed that 
the SMC method outperformed the linear controller in the presence of disturbances 


35 



36 


and uncertainties. Mudge and Patton [64] designed SMC laws with eigenstructure 
assignment for a remotely piloted vehicle and provide frameworks for assessment of 
robustness. Singh [41] defined a switching function which included the integral of 
the tracking error also and led to improved closed loop performance. Control laws 
for large rapid rotational maneuvers of aircraft under the perturbational model for 
aerodynamics and matched uncertainty are developed. Hedrick and Gopalswamy [42] 
show that the undesirable zero-dynamics could be tackled by output redefinition. They 
designed control laws for the nonlinear pitch axis system dynamics with linearized 
aerodynamics under matched uncertainty. The HQ criteria were incorporated into the 
design procedure. They also developed a method of selection of sliding gains for a 
MIMO system. Thukral and Innocenti [43] designed SMC laws for a linearized model 
of the longitudinal dynamics for the stabilization and tracking problems. They also 
presented a review of related work. 

In this chapter, the longitudinal or the pitch axis control problem for the nonlinear 
aircraft dynamics has been studied. The aircraft longitudinal dynamics and the force 
control approach (FCA) are explained briefly in the next section. In addition to the 
uncertainties and disturbances generally considered, an error between the command 
force vector and the actually obtained force vector is taken into consideration. A 
model following (MF) structure is next described and is used to design a full envelope 
SMC law taking these uncertainties into account. Analytic control laws for the tracking 
and the trimming modes of the aircraft are developed and simulations are done using 
the inertia parameters of the F-16 aircraft to show the effectiveness of the MF-SMC 
(model following-SMC) approach. 

3.2 Longitudinal Dynamics and Force Control Ap- 
proach 

The aircraft longitudinal dynamics in the body fixed frame is defined by the set of 
ordinary nonlinear differential equations [10] as given in (3.1). 
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(3.1) 


-QW - go sin(0) 

QU + pocos(0) 

0 

<5 

C/cos(^) + T^sin(0) 

\ U sm(9) - W cos{e) ) 

The system states are U, W, the velocities along the body X and Z axes respectively, 
the pitch rate Q, the pitch angle and Px, Pz the positions along X and Z axes 
respectively from an earth frame of reference. F^, Fz represent the forces along the X 
and Z directions respectively; and M is the pitching moment. Also ^'o denotes earth’s 
gravity, Rm denotes the inverse of the mass of the aircraft, and Cj is the inverse of 
inertia of the aircraft about its pitching axis. 

The forces and the moments in (3.1) are in general nonlinear functions of the state 
and the deflections of the CFG (conventional flight control) surfaces i.e., the elevator 
and the throttle [10]. However, the SMC design technique has mainly been developed 
for systems that are input affine. Most of the FCS applications cited earlier are based 
on input affine or linearized models of the aircraft. Bartolini and Zolezzi [65] extended 
the SMC method to include the case of systems nonlinear in the input by redefining 
the equivalent control, based on the existence of a particular map. This approach has 
been used in this work through the FCA, which is discussed next. 

If the forces Fx,Fz, and the pitching moment M can be independently varied, these 
may be treated as the pseudo input vector, and would render the governing system of 
equations in (3.1) input affine. Independent control of the 3-component force vector 
with the two CFC inputs is however not possible, and therefore, it becomes necessary to 
use additional control effectors. Smith and Meyer [23] treated the desired accelerations 
as pseudo-inputs. Hauser et al. [66] assumed that there exists a one-to-one mapping 
from the input space to the force vector. Snell et.al. [49] used TV (thrust vectoring) 
to augment the standard control surfaces. Franklin [39] used the forces as input for 
the translational longitudinal control problem. The TV concept is gaining increasing 
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popularity and in fact is going to be present in almost all of the next generation 
fighters [6, 8]. It has the advantage that its effectiveness does not vary much with 
the aircraft speed as opposed to the conventional aerodynamic surfaces. Hence, it 
may be used to advantage when the CFG surfaces are ineffective or are damaged in 
combat [46]. The TV-FCS has been flight tested on scaled aircraft models and is 
already equipping the front line USAF F-15,F-16, and YF-22 flghters [7], the X-31 
experimental aircraft [45], and the Russian Sukhoi family of flghters [9]. With TV 
capability, the forces Fx, Fz and the pitching moment M, with the assumption that 
the thrust acts through the center of gravity, are given [46] as 

cO 

Fx = qS\Cx{a,5,{) + ^Cx^[a)^+Tx, 

r cO ^ 

Fz = qS Cz{a,5ei) + -;^Czq{o') +Tz, and 

L ZVt } 

r cO 

M = qSc Sgi) + -^^CmQioc) + Cmrv^v , (3-2) 

where, g, S, a, 6ei, c, Vt represent the dynamic pressure, the wing surface area, the 
angle of attack, the elevator deflection, the mean aerodynamic chord, and the total 
velocity. Also the 5y represents the effective pitch thrust-vectoring angle. The various 
^i(‘! •) give the functional relations of the component of the non-dimensional force 
with its arguments (Refer [10] for further details). Also Tx and Tz are the components 
of the thrust vector in the directions specified by their subscripts, and are given as 
Tx = 7]TaCos{ST,) and Tz = r]TaSm{5^), where 77 is a dimensionless nozzle thrust coeflS- 
cient, and Ta is the actual thrust of the engine [46]. Thus the force vector is seen to be 
a function of three variables 5ei, the and the vectoring angle (5„. Other more conven- 
tional alternatives, such as the air brakes, spoilers, flaps etc. may also be used. Thus 
using a suitable combination of the available effectors, it should be possible to control 
the 3-component force vector. Durham [48] considers the control allocation problem to 
generate a set of desired accelerations. In case, many solutions exist, a suitable com- 
bination of these alternatives may be used to perform other control configured vehicle 
(CCV) functions. For this work, it is assumed that the aircraft has sufficient number 
of control effectors so that the 3-component force vector may be used as pseudo-input. 
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With this, the governing equation (3 1) becomes input affine and can be seen to 
depend on parameters such as Rm and cr Depending on the mission objectives of the 
aircraft project, a fair estimate of these parameters can be made using data from other 
aircrafts of a similar category It is better to define these parameters as a range of 
values they are likely to take This would have to be taken as such even if the aircraft 
parameters were known, because fuel and armaments carried on-board vary with time 
and also from mission to mission These have been incorporated as time varying but 
bounded multiplicative uncertainties as 

Rm ~ Shi^Rfii, and 

Cr = (1 + Sb2)cr (3 3) 

where Rm and cy represent the nominal values of Rm and Cy respectively The time 
varying multiplicative uncertainties represented by ^(,1 and Si ,2 are bounded as |^6i[ < 
Dbi, and |^ft 2 | < Db 2 la this section the aircraft longitudinal dynamics system, as 
considered in this work, has been explained As discussed m Chapter 1, a model- 
following approach has been used for this problem and is presented next 

3.3 The Mo del- Following Control 

A MF structure as shown in Fig 3 1 is used to achieve robust tracking performances 
The nonlinear aircraft state equations (3 1) have already been presented m the previous 
section The method also requires a model response traj’ectory This can be generated 
by the pilot with the inputs available to him namely the elevator, the throttle gear and 
the mode selector The mode selector specifies the outputs to be tracked DiflFerent 
models may be selected, for the various modes and with varymg operatmg regimes, 
as specified by the HQ engineers and/or the pilots It suffices to assume, that one is 
provided with the state trajectory to be tracked and its time derivative The model 
used m the present study is given m the Appendix D 

As can be seen from the figure, the controller gives as output the desired force vector 
Fcom which should be applied to the aircraft An aerodynamic inversion is required 
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Trim Variables 



Fig. 3 1 Schematic Model-Following Control Scheme 
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to get the actual control surface deflections On applying these to the actual aircraft, 
the generated force vector Fact is likely to be different from the Fcom as shown in Fig 
3 1 Most of the works cited earlier [23, 39] assume that these signals are identical 
However, m practice this assumption is likely to lead to errors for the following reasons 

1 Uncertainties in the aerodynamic modeling and data: The data used for 
control design is normally based on some models and/or observations from wind- 
tunnel experiments on the scaled models of the aircraft Also the aeroelastic 
effects are generally not considered 

2 Uncertainties in the aerodynamic inversion block: The nonlinear mversion 
process itself may not be exact 

3 Unmodeled actuator dynamics The actuator dynamics are generally ne- 
glected in a first level control design 

4 Uncertain position of the center of gravity (CG): The moment M depends, 
besides other factors, on the actual position of the CG 

5 Uncertcdnty in the air velocity relative to the aircraft: The forces and the 
moments vary proportionally with the square of the relative air velocity Errors 
in measurement and the presence of wind gusts can lead to this uncertainty 

It is assumed that errors from the above sources are bounded and can be represented 
by time varying multiplicative uncertainties as 

Fact = (/ + Aae)Fcom (3 4) 

where Aae is a 3 x 3 diagonal matrix with the diagonal entries being Saei^), i= 1, 2, 3 
all of which are bounded with |(5ae(^)| < Daei^), a known quantity 

The Section 3 2 and this section completely define the FCS design problem with the 
system equations m (3 1), the uncertainties involved m the aircraft inertia parameters 
and the inputs m (3 3) and (3 4), the desired model response (as m Appendix D), and 
the MF control structure used (as shown in Fig 3 1) The SMC design method for the 
MF problem of a MI nonlinear system as presented above is next discussed 
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3.4 Tracking in Multi-Input Nonlinear Systems 

The aircraft system of equations in (3 1) with external disturbances can be expressed 
in a general form as 

X = f(x) + G(x)u + d, and 

y = y(x), (3 5) 

where x G 7^”, u, y € represent the state, input and output vectors respectively, 
and d represents the external disturbance vector given bounded as \dt\ < i = 
1, ,n Since the number of trackable outputs equals the number of independent 
inputs, the number of outputs and inputs are assumed to be equal Also G(x) is 
assumed to be of rank m over the entire state space The number of differentiations of 
y(x) required for the input u to appear in its expression is called the relative degree 
of the system It is assumed without any loss of generality that all components of y 
have a relative degree T’ 

Under a suitable coordinate transformation from x to z = where the 

superscript T denotes the transpose, and ^ G the above nonlmear system may 

be expressed in a reduced form [37, 52] as 

y = fi(z) + Bi(z)u + di, and (3 6) 

^ = f2(z)+d2 (3 7) 

The (3 6) represents the matched dynamics and the disturbances and/or uncertamties, 
dj, entering this equation are called matched The (3 7) gives the unmatched portion 
of the dynamics and the disturbances and/or uncertainties, d 2 , entenng this equation 
are said to be unmatched In this chapter, only matched disturbances/uncertainties are 
assumed since the SMC guarantees invariant performance m the presence of this class of 
uncertainties The case of unmatched uncertainties is considered in the next chapter. 
The matched disturbance vector di is given to be bounded with [dij < JDi,, i = 
1, ,771 Under the assumptions made above the matched dynamics represents a 

square system of m equations Also Bi(z) is an m x rri matrix having full rank over 
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the state space z thus implying a strong relative degree The arguments of the 
various functions defined above are dropped if no confusion is likel}' to arise Let Bi 
represent the nominal Bi matrix It is assumed that the Bi may be expressed as 
Bi = (/ + Ai)Bi, m which all the matrices are appropriately dimensioned with I 
representing the identity matrix Also the matrix A^, is bounded b^ a known matrix 
Db as |A6(^,;)| < A(hj) 

As explained earlier, u in (3 6) is the actually applied control signal, i e Uact With 
the uncertainty as modeled in (3 4), Uact = (/ + Aae)ucom , the system equations in 
(3 6) and (3 7) can be rewritten as 

y = fi + (J 4- A6 )Bi(/ + Aae)ucom + di, and (3 8) 

e = f2 (3 9) 

The MF-SMC design for the case when Aae = 0 is reported in Slotine and Li [21] and 
Hedrick and Gopalswamy [42] In this chapter the method has been extended to a 
more general case with Aae ^ 0 

Let the trajectory to be followed be specified as jd and jd Define S as a vector 
consisting of the error function for the output vector as 

S = y - yd (3 10) 

It IS desired that S = 0 To this end, taking time derivative of the above equation, and 
using (3 8) gives 

S = y-yd 

= fl + (-f + A6 )Bi(/ + Aae)Ucom + di — y^ (3 11) 

The equivalent control [21] is defined as the input Ucom , which makes S = 0, when 
S = 0 for the nominal plant, i e Ab = 0, Aae = 0, and di = 0 Thus the equivalent 
control becomes 

flcom = Bj;^(yd - fi) (3 12) 

The applied input function should make the SS attractive m the presence of distur- 
bances and uncertainties also This is achieved by selecting the apphed control signal 



44 


as the sum of the equivalent control and a discontinuous function as 

Ucom = - fi - k sgn[S]) (3 13) 

where, k sgn[S] is an m x 1 column vector with all A:i > 0, z = 1, ,m Let 
(k sgn[S])i = feiSgn[S'i] represent the component of the vector The real vector 
k has to be selected so that trajectories starting from nonzero values of S approach 
the SS Putting Ucom from (3 13) in (3 11) gives 

S = fi + (J + A6 )Bi(/ + Aae)B;'^(yd - fi - k sgn[S]) + di - (3 14) 

The above equation may be simplified by defining an m x m matrix At, such that 

(/ + At) = (J + A6)Bi(/ + Aae)Bri (3 15) 

The bounds on the matrix At can be evaluated from the known bounds on At and 
Aae as Will be shown in the next section Let these be given as |At(z,y)| < Dt{t,j) 
With this, (3 14) can be rewritten as 

S = fi + (7 + At) (yd - fi - k sgn[S]) + di - yd, (3 16) 

which may further be simplified to 

S = At(yd - fi) - (7 + At)k sgn[S] + d^ (3 17) 

The SS IS made attractive by selecting k such that 

5.5. < -77,15,1, (318) 

where 77.,z = 1, ,m is a positive scalar quantity to be selected, thus satisfying the 
reaching condition [21] The above equation guarantees (see Appendix A) that the 5. 
response will be better than the response of 

5. = -77.sgn[5.] (3 19) 

Following the algebraic manipulations along the lines of [21, 42], k can be obtained as 
the solution of the following system of linear algebraic equations 

(3 20) 


(7-A)k = ^, 
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where, Dj is the matrix with as its components, and is an m x 1 vector with 

m 

$. = »(, ) = + + • =1. (321) 

J=1 

It is obvious that > 0, i = 1, ,m from the above equation The existence of 
a unique solution k with all Ajj > 0 is easily established using the Forbenius-Perron 
Theorem [21] (See Theorem 1 m Appendix G) and requires that the largest positive 
eigenvalue of Dt be less than 1 0 

The use of a sgn[ ] function in (3 13) gives rise to the chattering problem Various 
methods have been suggested m the literature for its alleviation [21, 56] In this study, 
the concept of a boundary layer as proposed by Slotine [21] has been used to approx- 
imate the discontinuity by a continuous control law Accordingly, for implementation 
purposes, the vector k sgn[S] in (3 13) is replaced by k sat[S/$], whose element 
is given as Aij sat[5i/^i], where $ is a m x 1 column vector of appropriately selected 
boundary layer thickness, given as [^j]"^, i = 1, ,m The control law in (3 13), then 
becomes 

Ucom = (yd - f - k sat[S/$]) (3 22) 

The Si dynamics then becomes 

Si = -riiSat[Si/(f>i] (3 23) 

With the Ucom as selected above, the motion in the sliding mode is guaranteed to be 
inside of thickness $ of the SS [21] Thus, 


\Si\<(j>i (3 24) 

Moreover, starting from a non-zero initial condition So, the dynamics m (3 23) guar- 
antees that its z*^ component Si reaches the boundary layer within a finite time (toj 
given by 

to < (3 25) 

’ Vx 

Thus, a larger value of 77 leads to a smaller reaching time and larger control efforts 
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Thus, the above selection of Ucom m (3 22) guarantees that the outputs are tracked 
with a desired accuracy, however, the stability properties of the unmatched dynamics 
in (3 9) need to be investigated It has been shown [21] that this can be determined 
by studying the stability of the zero dynamics (ZD), which is defined as the dynamics 
of the system when the outputs and all its time derivatives are kept at 0 (zero) by 
appropriate choice of the input The ZD for the system (3 5), with all outputs having 
a relative degree T’ is thus seen to be a special case of the unmatched dynamics in 
(3 9) as 

e = f2([0W) (3 26) 

The stability of the unmatched dynamics is proved if the ZD is at least bounded [21] 
The SMC method as outlined here is applied to the design of control laws for the 
longitudinal dynamics of an aircraft using the MF structure as discussed m Section 3 3 
Therefore, the method is termed as the MF-SMC The results of numencal simulation 
are presented next 

3.5 Numerical Example and Simulation 

MF-SMC laws are synthesized, following the steps outlined in the previous section, for 
the FCS design problem as presented in Sections 3 2 and 3 3 The numerical values 
of the time varying multiplicative parametric uncertainties and disturbances and their 
bounds as defined in the previous sections as also the mass and inertia parameters of 
the aircraft as used for simulation purposes are given in the Appendix D To show 
the applicability of the proposed method, controllers are designed for the following two 
cases 

1 Cl: In this case, the controller design is done assuming that the Fcom and F^ct 
are identical, i e Aae = 0 

2 C2: Here, the controller design takes into account that Aae 7^ 0 

The controllers corresponding to the cases Cl and C2 will be called as the Ml and 
M2 respectively 
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As stated earlier, it is possible to track only three states or outputs using the 
3-component force vector Thus various operational modes for the aircraft can be 
specified, depending on the variables which are of concern during a given flight phase 
However, any trajectory of interest can be tracked using a suitably generated velocity 
signal denoted by Ud, Wd, and Qd where the subscript ’d’ stands for desired trajectory 
Trimming of an aircraft at a level flight condition is of vital importance to the pilot 
The trim variables are also required to evaluate the desired trajectories as can be seen 
m the Appendix D The responses obtained by the two controllers (Ml and M2) 
are then compared for the Tracking and Trimming modes The velocity tracking is 
considered first 


3.5.1 Controller Design for Tracking Mode 


In this case, the state equation (3 1) including the disturbances can be seen to be 

f U\ i -QW-gQSin{e) \ ( RmF^\ ( ^ 

W QU + goCos{9) 

Q __ 0 CjM 

ego 

Px U cos(9) + W sm(9) 0 

\ Pz J \U sin(0) - W cos(6) y \ 0 / \ ^ / 


+ 


du 

dw 

dq 

0 

0 

0 


(3 27) 


where, dj, i = U^W,Q represent the external disturbances and uncertainties entering in 
the channels denoted by the subscript i and their magnitudes are bounded by ld,| < £>i 
It is seen that, with y^ = [17 kT Q] as the output vector, the above equation is already 
in the reduced form, with the first three components forming the matched dynamics 
(3 6) and the remaining three the unmatched dynamics (3 7) The ZD (3 26) in this 
case IS thus given by 

( 9 

px 
Pz 


^0^ 
0 

VO/ 


(3 28) 
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Thus, the states of the ZD are in a neutral equilibrium, implying that the unmatched 
dynamics are bounded-input bounded-state (BIBS) 

For both cases (Cl and C2), the Bi matrix is seen from (3 27 to be 

Bi = diag[^, Rm, cr], (3 29) 

where diag[ ] denotes a diagonal matrix Also, the matrix is seen to be diag[(l 6 i, 

The vector S = [Su, Sw, 5'q]'^ is as defined in (3 10) Following the steps as outlined 
in the previous section, I + At in (3 15) can be simplified as 

(J -I- At) = diag[gi, ^ 2 , Qsl (3 30) 

where, qi = (l-H56l)(l+<^ae(l))! 92 = (l + ^ 6 l)(l + <Iae(2)), and 93 = (1 + ^62)(1 + ^ae(3)) 
The At matrix for the case Cl can be obtained by putting 6ae{i) = 0 for i = 1, 2, 3 In 
this case the Dt matrix, which bounds the At matrix m (3 30), is seen to be given by 
the bounds on 5bi and 5b2 (given in the Appendix D) as 

Dt = diag[l/3, 1/3, 3/7] (3 31) 

For the case C2, the /-|- matrix in (3 30) would attain its maximum (minimum) 
values when both the terms in the product individually take their maximum (minimum) 
values Using the bounds as given m Appendix D, the elements of the diagonal matrix 
/ -h At can be seen to vary within the following bounds 

— :<1 + Af(l, 1)< - 

— <1-1- Af(2, 2) < - 

^ <1 + A,(3,3)< y (3 32) 

from which, the Dt matrix bounding At in (3 30) becomes 

Dt = diag[3/5, 3/5, 5/7] (3 33) 

It is obvious that in both cases, all the eigenvalues of A are less than 1 0, thus guar- 
anteeing the existence of k in (3 20) The positive scalar ry, in (3 18) has been taken to 
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be 1 0, i = 1, ,3 With these choices, the 3x1 vector ^ in (3 21) can be seen to be 


^-1 =^( 1 )= lQ + Du-^Dt{l,l)\Udi-QW + gQsm{6)\, 

'^2 = ^^(2) = 1 0 + Dw + A(2, 2)|Wd - (Jt/ - £rocos(0)|, and 
^^3 =^(3)= 10 + DQ+A(3,3)|(5d| 


(3 34) 


Note that for the two cases Cl and C2, different Df as m (3 31) and (3 33) respectively 
have to be used 

The 3x1 vector k in (3 20) can then be evaluated for the two cases as 


1 Cl: kT = [1 5^1, 1 5^2, 1 75 T 3 ], and 

2 C2. k^ = [2 ST'i, 2 5^2, 3 5 T 3 ] 


The calculated values of k are then used in (3 22) to obtain the control law as 


^fcom 


2com 


M., 


1 

An 

1 

An 


Ud + QW + go sm{6) - A:isat[^] 

^ rU j 


Wd — QU - go cos(^) — A: 2 sat[^^] 1 , and 


4>w 


Qd 


ks sat[^] 
9q ) 


(3 35) 


The boundary layer thicknesses {(f>xj, (f)w and i^g) for the three outputs as used in 
the simulation are specified in the Appendix D It can be seen that the state variables 
Px and Pz are not used in the control law, however they constitute the ZD, and hence 
they have been mcluded in the above development Simulation has been done using 
a fourth order Runge-Kutta algorithm with the time step of 0 25ms The aircraft is 
initially assumed to be in a trimmed level flight with a total velocity of 502 feet per 
second (fys) Till 0 5s, the elevator and throttle gear commands are mamtamed at 
their equilibnum values, and at this time instant the pilot increments the elevator by 
15° and the throttle gear by 20 units and maintains them at a constant value till 2 Os 
The responses were obtained for two sets of uncertainties (UCl and UC2, as spec- 
ified in the Appendix D) for the two cases (Cl and C2) The tracking error of the 
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time (s) 

Fig. 3.2 Tracking error in U (fps) 


outputs {U, W and Q) are shown in Figs 3 2, 3 3 and 3 4 respectively In these figures, 
the solid lines labeled ’A’ and ’C’ correspond to the responses using the M2 controller 
for the uncertainties UCl and UC2 respectively Similarly the dashed lines ’B’ and ’D’ 
show the results using the controller Ml for UCl and UC2 respectively It is observed 
from the figures, that the output tracking error stays well within the specified bound- 
ary layer thickness and good tracking of the desired trajectory is achieved. Prom the 
figures, it is clear that the tracking errors in the case C2 are smaller than those of the 
case Cl for all the outputs Although, 6 is not a tracked output, in the absence of 
unmatched uncertainties, it was found to be well tracked with a similar trend Thus, 
the controller M2 clearly gives better performance as compared to the controller Ml 
Thus it is seen that ignoring the uncertainties involved in the force generation can lead 
to a degraded system performance 

The desired force command and the actually applied force vector, with the assumed 
uncertainty in aerodynamic inversion /modeling is shown in Figs 3 5, 3 6, and 3 7 for 
the UCl uncertainty and case C2 
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Fig. 3.7 The command M {Mcom ) and the actually realized force (Mact ) 


3.5.2 Controller Design for Trim Mode 


The trimming of an aircraft at a level flight condition requires that an equilibrium 
condition is reached at any given velocity Most earlier works require a trim map or 
some algorithm [23] to generate the trimming variables With the proposed design 
method this is not required Given the velocity at which the aircraft is to be trimmed, 
the control law leads the aircraft states to the desired trim variables The pilot is 
provided with a trim switch to trim the aircraft at the current value of the total 
forward velocity Continuing with the simulation results in the previous section, the 
pilot activates the trim switch at f = 2 Os The tnmmed or the equilibrium state is 
defined as 




{Vo\ 

7 

= 

0 

V ^ ) 


1 0 J 


(3 36) 


where, Vo is the desired trim velocity In the above equation Vt = (iP + t^ 2 ^i /2 
represents the total velocity, a = (W/U) is called the angle of attack, and 7 = 




54 


6 — a IS the flight path angle It is assumed that there are appropriate sensors to 
measure these quantities Also define S = [5yt, S-y, Sq]^ as the error function for the 
outputs as m (3 10) To track the outputs as defined in (3 36) the dynamic equations 
governing them need to be defined These are given m standard textbooks [10] or can 
be directly obtained from the definitions of the outputs as above, and using (3 27) 
The state equation in the new set of coordinates, = [14 7 Q ^ Px Pzli may be written 
as 


fv.-' 


^ -^osm(7) ^ 


i?m cos (a) 

Rm sin(a) 

0 

7 


-go cos (7) 


i2mSm(Q;) 

-Rm cos(a:) 

0 

Q 



0 

+ 

0 

0 

C7 

9 


Q 


0 

0 

0 

Px 


14 cos (7) 


0 

0 

0 

[pz y 


^ V[sin(7) y 


0 

0 

0 


/ P \ 


VMy 


+ 


( dvt ^ 

diiy 

dq 
0 
0 
0 


V ^ J 


(3 37) 


where, dvj, dy, and dq represent the matched external disturbances and can be evalu- 
ated as 


dvi = dtrcos(a:) -f dwSin(Q:), 
d~y = df/sin(a) — dvKCOs(a!), and 

dq = dq (3 38) 

Their bounds change to 


Dvt 

= Du -f- jDvv, 


D, 

= Du + D\y, and 


Dq 

= Dq 

(3 39) 


It IS seen that, with = [1^ 7 <3] as the output vector, (3 37) is already in the reduced 
form as before with the first three equations defining the matched d)Tiamics and the 
remaining unmatched dynamics Also the ZD m this case are the same as m (3 28) 
Thus the unmatched dynamics is BIBS 
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For the system (3 37) the input distribution matrix (denoted by B'J, corresponding 
to the Bi in (3 6), is defined in the first three equations in (3 37), and the nominal 
matrix is seen to be given as 


Bi 


Rm cos(a;) Rm sm(Q;) 0 
RmSm{oc) —RmCos{a) 0 
0 0 C7 


(3 40) 


for both the cases It is easily verified that B'^ matrix above is invertible over the 
entire state space and the chosen parameter range The matrix is the same as in 
the tracking mode Following the steps as in the previous section, / + At in (3 15), can 
be obtained as 


(I + A,) 


(1 + <^61)91 |(1 + <^61)92 0 

|(1 + <^61)92 (1 + <^61)93 0 

0 0 (1 + (562)(1 + <5^ae(3)) 


(3 41) 


where, gi = (1 + 5ae(l) C0S^(Q!) + 5ae(2) Sin^(Q;)), 52 = (<^ae(l) - ^ae(2)) sin(2Q;), and 
93 = (1 + 5ae(l) Sin^(Q!) + 5ae{2) cos^(a)) 

Note that for the case Cl with Aae = 0, 51 = 53 = 1 0 and 92 = 0 With this, the 
matrix Dt for this case is seen to be the same as in the tracking mode and is given m 
(3 31) 

For the case C2, with Aae(?) 7^ 0 using the bounds as given in the Appendix D, it 
IS seen that the 51, 52 and 53 as defined above are bounded as 

6 

<9i< 7 

-«2_ g 

< 93 < ^ (3 42) 



Proceeding as in the previous section, the components of the symmetric I + Aj matrix 
in (3 41) are seen to assume values as given by 

^ <[/ + A,](U)< I 
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if + ^ 

^ < 1/ + A,l(2, 2) < ? 

^ <[/ + A,J(3,3)< y (3 43) 

From the above equation, a 3 x 3 bounding matrix Dt can be selected as 

3/5 4/15 0 

A = 4/15 3/5 0 (3 44) 

0 0 5/7 

It can be easily verified that for both the cases, all the eigenvalues of Dt are less than 
1 0 The T vector (3 21) in this mode becomes 

= 1 O + IPi/j +£?oA(l,l)|sin(7)| + <7oA(1,2)|cos(7)|, 

^2 = 1 0 + A +i?oA(2, l)|sin( 7 )| +5'oA(2,2)|cos( 7)|, and 

^3 = 10 + A+A(3,3) (3 45) 

with the time derivatives of the outputs (i e Vt^, 7d, and Q^) taken to be zero Using 
the Dt m (3 31) and (3 44) for the cases Cl and C2 respectively, the k m (3 20) can 
be evaluated in the two cases as 

1 Cl: = [1 5Ti, 1 5^2, 1 TSTa], and 

2 C2: kT = [4 5^1 + 3 0T^2, 3 + 4 5^2, 3 6 ^ 3 ] 

Using the k as above and the inverse of the B'^ matrix in (3 40) m (3 22) gives the 
control law as 

Aeon. = sm{e) - i [ ^1 sat[|^] cos(a) + ^2 sat[|^] sm{a) , 

Rm Rml 9Vt <p7 J 

Aeon. = CQS{0) - i sat[f^] sm(Q;) - h sat[^] cos(a:) , 

Rm Rml <PVt } 

Mcom = -^Mat^, (3 46) 

C7 <Pq 
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where, the boundary layer thicknesses used are specified in the Appendix D It can be 
easily seen that the command forces (-F^com ^zcom ) consist of two terms, the first of 
which gives the equilibrium force component for the nominal aircraft, and the second 
term is a correction factor based on the errors in the output 

At the end of the simulation as shown in the previous section, the trim switch is 
activated Simulation results are presented for an interval of 2 Os from 2 0 to 4 Os 
The tracking errors for the outputs (Vt and 7), the angle of attack a and the altitude 
variation are shown in Figs 3 8, 3 9, 3 10 and 3 11 respectively The labels ’A’-’D’ are 
as explained for the trackmg mode It was noticed that the pitch rate Q attains very 
small values almost instantaneously m both the cases and hence it is not presented 
here It can be seen from Fig 3 8 that the control law M2 takes less than 1 Os to trim 
the aircraft for both uncertainty cases and that the error in Vt remains well inside the 
boundary layer However, the control law Ml is seen to be unable to trim the aircraft 
even till 2 0 seconds, and in fact for the UC2 uncertainty, the error in Vt is seen to go 
outside the boundary layer This can also be seen from the slow response of 7 (Fig 3 9) 
obtained in case Cl as compared with C2 Similarly the a variation (Fig 3 10) shows 
that fast responses are obtained in the case C2 From, the altitude variation (Fig 3 11) 
it IS seen that for the controller M2, the aircraft looses an altitude of about 70 and 
58 feet for the uncertainty sets UCl and UC2 respectively Comparatively the aircraft 
drops approximately 244 and 225 feet with the Ml controller without achievmg the 
trim state in the same time Thus, the M2 controller is seen to give vastly improved 
performance over the Ml controller 

Thus it is observed that the MF-SMC gives fast and robust trimming performance 
It IS evident from the simulation results that the proposed design method is robust to 
the uncertainties that can creep into the aerodynamic inversion block (see Fig 3 1). 
Moreover, the method does not require a trim map or an algorithm to find out the trim 
variables It eliminates the necessity of developing a dynamic trim map and reduces 
the pilot workload of actively tnmming the aircraft, and thus represents considerable 
improvement over the prevailing practice 
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3.6 Conclusions 

f 

In this chapter, the robust nonlinear control law design for the nonlinear longitudi- 
nal dynamics of a fighter aircraft, based on a model following sliding mode approach 
and using a full state feedback, is presented The components of the force vector are 
considered as inputs to the aircraft, with the assumption that a desired force vector 
can be generated using appropriate control surface deflections An error between the 
commanded force vector and the actually obtained force vector is considered, which 
is m contrast to the assumptions made in the literature The MF-SMC method is 
extended to include this uncertainty Simulation results show that ignoring this uncer- 
tainty m the controller design may lead to substantial loss in performance All other 
uncertainties and disturbances are assumed to satisfy the matching conditions Ana- 
lytical expressions for the control laws for two operational modes, namely the velocity 
tracking and the trimming are derived It can be seen from the developed control 
laws for the two modes in (3 35) and (3 46) that they depend on the estimated mass 
and inertia parameters and their expected range of variation only In a preliminary 
aircraft design project, these control laws can thus be used to generate the required 
forces and moments for some critical flight maneuvers which can be placed as targets 
to the aerodynamics and the propulsion groups This would enable the control group 
to have a more active role in such an exercise, and help m achieving the best possible 
performance from the airframe Numerical simulations on a fighter aircraft show the 
robustness of the control law An interesting feature of the trimming control law is 
that it does not require an algorithm or a trim map to specify the trim variables 



Chapter 4 

Multi-Loop SMC Laws for The 
Unmatched Uncertainty Case 


4.1 Introduction 

In the previous chapter, controller design for output tracking with matched uncertainty 
and disturbances was considered As seen there, the number of trackable outputs are 
limited by the number of independent inputs available This leaves out a portion 
of the dynamics, also called the internal dynamics, and is seen from (3 7) to be a 
function of the tracked outputs and the internal states An approach to control these 
internal states was to formulate the aircraft dynamics as a two time scale problem and 
use a multi-loop (ML) strategy The state vector was thus separated into fast and 
slow variables [49, 50] and a dynamic inversion based inner loop control law was then 
designed for the fast variables With the command signals for the fast variables as 
the input, an outer loop control law is designed to track the slow states, which gives 
the command signals to be tracked in the inner loop The inversion based control law 
was reported to give better responses as compared to a gain-scheduled linear law [49] 
Analysis for robustness was, however, not performed Azam and Smgh [50] used an 
integral error feedback m the outer loop for robustness against parameter variations 
Both of the works cited above use conventional control surfaces as inputs and hence 
were hmited to track the angular velocities in the inner loop and corresponding angular 
states in the outer loop However, with the force control approach (Section 3 2), control 
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of the entire state vector is possible, as is shown in this chapter 

The presence of unmatched disturbances and uncertainties, 6.2 m (3 7) in the in- 
ternal states, which was neglected in the works cited above, is a practical problem 
and leads to a degraded system performance This chapter deals with the design of 
SMC laws m the presence of unmatched uncertainties using the ML strategy A basic 
assumption used m this work is that the nominal unmatched dynamics are assumed to 
be acceptable Accordingly, tracking of the position state variables is accomplished in 
the outer loop considering the command velocity signals as the inputs A SMC law is 
designed m the inner loop so that the velocity states track these augmented velocity 
commands generated in the outer loops The ML strategy is shown to lead to a simpler 
propagation of disturbances m the outer loop to the inner loop using which robustness 
properties are evaluated in terms of the bounds on the tracking errors m the inner 
and the outer loops Simulation results for the nonlinear longitudinal dynamics of an 
aircraft are then presented The ML strategy is next applied to a higher order SISO 
system with unmatched disturbances and uncertainties The developments herein can 
be seen to be giving results similar to those of [67, 68] which uses the singular perturba- 
tion methods (SPMs) The SMC design method presented here requires the knowledge 
of the bounds as in the previous works [67, 68] It is shown that the conventional SMC 
design, with the modified way of selecting the switching gams, may be used to handle 
this problem 

4.2 ML Strategy for Mechanical Systems 

Many mechanical systems may be represented in the reduced form as in (3 6) and (3 7) ' 

as 


V = fi(X,V)+Bi(X,V)u + di, (41) 

X = f2(X,V)+d2 (4 2) 


where X, V and u E 7V^ represent the position and velocity state variables and the 
input vectors respectively, and di and d 2 represent the combined influence of the exter- 
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nal disturbance and uncertainty vector and are given to be bounded representing 
the component of the disturbance vector d, is bounded as | < 1),^ , z = 1, 2, j = 
1, , n Let Z = [X'^ Also Bi(Z) is assumed to be of full rank over the entire 

state space The other assumptions in Section 3 4 for matrix Bi hold in this chapter It 
IS clear from the above system equation that with the velocity variables as the output, 
it has a strong relative degree of 1 Also, the above nonlinear system is seen to be in 
the reduced form [37, 52] The arguments of the various functions defined above are 
dropped if no confusion is likely to arise 

As m the Section 3 4, (4 1) represents the matched dynamics and (4 2) gives the 
unmatched portion of the dynamics The ZD, for the class of systems defined above, 
were seen to be a special case of the unmatched dynamics (3 26) These equations 
in such systems generally are integrators of the velocity variables and as seen in the 
aircraft examples m Sections 3 5 1 and 3 5 2 are BIBS It is also seen from the system 
equation in (3 27) m Section 3 5 1 that the term f 2 (X, V) is actually separable in X 
and V, 1 e it may be represented as B 2 (X)V Thus (4 2) may be rewritten as 

X = B2(X)V + d2 (43) 

The earlier control laws (in Section 3 4) were designed neglecting the presence of 
d 2 In this section, a ML strategy for the SMC is presented It is assumed that the 
desired trajectories for the velocity variables are specified (V^, V^) and the associated 
desired position trajectories gven by the nominal X dynamics (4 2) as 

X, = B2(Xd)V, (4 4) 

are acceptable For notational simplicity, B 2 d will be taken to represent B 2 (Xd), and 
similarly B 2 = B 2 (X) from now on It is further required that B 2 be nonsmgular 
everywhere 

The ML control scheme is shown in Fig 41 As can be seen from the figure an 
inner loop is used to synthesize Ucom such that the velocity variables track an aug- 
mented signal Vfl instead of the conventional Yd The outer loop is used to create the 
augmented signal for the inner loop such that the position tracking is within acceptable 
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Xd 

Vd 



Fig. 4 1 Schematic Multi-Loop Control Scheme 

bounds Accordingly bounds on the tracking errors m the inner and outer loops are 
developed as follows 

Define Sx and as vectors consisting of the error function for the position and 
velocity variables as 


= X - Xd, and (4 5) 

S„ = V-V„, (46) 

where, Va is yet to be defined Taking the time derivative of and using (4 3), 
becomes 

S^ = B2V + d2-B2dVd (4 7) 

Using V = Va + from (4 6) in the above equation gives 

S,: = BaVa + B2S„ + d2 - BadVrf (4 8) 

Note that the time derivative of Va would be required for V to track Ya m the inner 
loop, therefore Vo should not have discontinuities and should be selected so as to make 
Sj; as small as possible With the invertibility condition on B 2 , the command signal 

I 

Vo for the inner loop is selected from (4 8) to cancel out the known terms and by 
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introducing a stabilizing feedback as 

V„ = B2-nB2dVd-K,S,), (4 9) 

where, K^Sa; € 72-" with its component as [Ka;Sj;]i = Kx,Sx, and iQ, > 0, ^ = 
1, , n Kx IS taken as an x n real matrix as diag[Kxi, , Kxn\ 

Using the Vq selected above in (4 8) gives as 

Si = 'KiSi + 628^ + d2, (4 10) 

which IS rewritten as 

Si = Ti + d2 (4 11) 

for simplicity, and Ti is obviously defined The 5i, is thus seen to have the dynamics 
given by 

5i, = —Kx^Sx^ + [B2Sa,]j + d2, (4 12) 

The assumption of f2 = B2V has been used in the selection of Vq m (4 9) However, 
only the existence of the following map guarantees that there would exist a unique Va, 
such that 

f2(X, V,) = f2(X,, V,) - KiSi (4 13) 

From (4 12), it is seen that any bounds on the position error Si depend on the 
velocity error S„ Hence, it is first necessary to check the S^ dynamics before any 
estimates of the Si bounds can be made The time derivative of the Vq in (4 9) would 
be required for V to track Vq m the inner loop Snell et al [49] used a band limited 
filter to generate this time derivative Azam and Singh [50] used the nominal system 
equations (4 1) and (4 3) to evaluate the derivative However, the effect of disturbances 
and uncertainties was not directly considered m their study The time derivative of Va 
IS obtained from (4 9) as 

V„ = ^[Bj'](B2dV, - KiSi) + Bji^psdV, - KiSi] (4 14) 

The definitions and notations used in the present work to calculate the time deriva- 
tives of scalar, vector, and matrix functions of vectors are given in Appendix E A well 
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known result in matrix theory [69] gives that for a nonsingular matrix B 2 

(4 15) 

Using the notation in (E 8) from Appendix E, the time derivative of B 2 in the above 
equation is seen to be 

|[B2) = *,[4®X1, (4 16) 

where <g) represents the Kronecker product From (4 2), it is observed that X = £2 H-d 2 
Substituting this m the above equation gives 

^[B2] = JbAIu ® £ 2 ] + JbAIu ® d2] (4 17) 

Using (4 15) and (4 17) and + d 2 from (4 11), Vq in (4 14) becomes 

v„ = -B2-‘^[B2lBj-'(BMVj-KA) + Bj-'(^[B2dVjJ-K.S,) 

= B^*1(|[B2jVJ - - K,d2) 

® f2l + Jb,% ® d2])Bj-‘(B2Wi - K.SJ), (4 18) 

which may be rewritten as 

Va = BJ' - JsAln ® £ 2 ]B 2 -i(B 2 <iV, - K,S,)^ 

-B 2 1 (K,d2 + JbAIu ® d2]B2-i(B2dVd - K,S,)) (4 19) 

The above equation may be expressed as 

V^ = Ta + d^, (4 20) 

where 

T. = B2-*^^(B2W<il-K.T.-JB,[/„®f2lBj*(B2Wj-K,sA,(4 21) 

anddi = -BJ-qK.d2 + J'j,,(/„®d2lB2‘(B2Wd-KA)) (4 22) 

It is seen from the above equation that d^ contains the effects of disturbance terms d 2 , 
which was not considered m [50] 
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The dynamics of S„ may now be considered The time derivative of m (4 6) is 
seen to be = V — Va, and using (4 1) and (4 20) it becomes 

S, =fi + BiU-Ta+di-d^ (4 23) 

Thus, it IS seen that a ML strategy passes on the effect of the disturbances d 2 in the 
outer loop to the inner loop as d^ (4 22) as in the above equation, which can then be 
incorporated m the design process Using the augmented disturbance signal, the SMC 
as in (3 22) may still be designed by following the steps outlined in Section 3 4 
can thus be made to satisfy the reaching condition as in (3 18) and the response is 
guaranteed to be better than 

Sv, = -?7t,.sgn[5„,] (4 24) 

as m (3 19) As in Section 3 4, the sgn[ ] term is replaced by a continuous term 
sat(5v,/0u.)) where represents the component of the selected boundary layer 
thickness vector ($„) With the reaching times to, as given in (3 25), all the 5„/s will 
be bounded after some finite time tom (— max[to,, i = 1, ,n], where max[ , ., ] 

gives the maximum of its arguments) as 

(4 25) 

With the velocity error S„ bounded as above, the effect of the [B 2 S„]i term in the 
Sx, dynamics in (4 12) is seen to be bounded by the sum and may be 

considered as an additional disturbance term for all t > tom can be shown that the 
error is globally uniformly ultimately bounded (see Appendix F) as 

IS.. I < 

Thus it is seen that the and S„ errors defined in (4 5) and (4 6) remain within 
bounds as specified in (4 26) and (4 25) respectively The bounds on (4 25) is 
directly governed by the choice of The bound on Sj, (4 26) may be lowered by 
choosing a higher value for the Kx,, however, it is observed that the bound on 
does not change It directly depends on the size of the unmatched uncertainties D 2 
and the selected boundary layer thicknesses 
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Another point of interest is to find out the magnitude of the difference between the 
Vd and Vfl signals From (4 9) it is seen that 

Vd - Va = [/ - B2-'B2d]Vd + (4 27) 

The difference Vd— Va is seen to be equal to 0 if and only if = 0, since it also implies 
that Bj ^Bad = I However, as is only guaranteed to be within some bounds, there 
will be a finite error between the two signals 

Thus it IS seen that the SMC method can be applied to the unmatched uncertainty 
case using a ML framework and bounds for the errors in the loops may be derived as 
shown here The method developed for a general mechanical system is next validated 
by applying it to the aircraft longitudinal dynamics 


4.3 SMC Design for Aircraft Velocity Tracking 


In this section, the ML controller design for the unmatched uncertainty case as dis- 
cussed m the previous section is applied to the aircraft example m Section 3 5 1 The 
state equations for this case are already reported in (3 27) The matched dynamics 
representing the velocity dynamics m (4 1) remain the same as earlier Thus, the V 
in this chapter represents the velocity vector [U, W, QY However, the unmatched un- 
certainties, which were not considered m Chapter 3, appear in the X dynamics as in 
(4 2) Rearranging the last three equations of (3 27) with the corresponding states 
being given by X = [pj., Pzi ^ governing equations may be rewritten as 


^ Px ^ 


^ UcoB{9) + Wsm[e) ^ 




Pz 

= 

U sm{e) -W cos{e) 

+ 

d2z 

(4 28) 

K ^ ) 


\ Q ) 


^ d^e ) 



The f 2 as defined m (4 2) is seen from the above equation as 



cos(0) sm(0) 0 



^2 = 

sm(0) — cos(0) 0 


W 


0 0 1 


{ Q J 


(4 29) 



69 


From the above equation, it is seen that the position state dynamics X as given m 
(4 28) is in the form of (4 3) It is easily verified that the matrix B 2 m the above 
equation is non-singular As in Section 3 5 1, the nominal ZD (3 28) is seen to be 
BIBS 

The error Sj, (4 5) m the position variables is defined as [5i, S^, and is given 

by 

■S'z I = I Pz-Pzd 

e-dd 


/c \ 


^ Px- Pxd ^ 


\SeJ 




(4 30) 


y 


Note that S^, represents a vector quantity, whereas the is a scalar quantifying the 
X position error Similarly, represents diagfifa,, K^, Ke] While Kx is a scalar 
quantity representing the feedback gam of Sx, the is a 3 x 3 matrix Similarly, the 
velocity error = V — Vq m (4 6), where the augmented signal = \Ua-i Wa, Qa]'^ 
will be chosen as in (4 9), represents the column vector [Su, Sw, 

With the B 2 matrix as defined in (4 29), the dynamics in (4 12) may be evaluated 
to get 


Sx = -KxSx + Su cos{6) 4- Sw sin(0) -f d 2 x, 

Sz = —KzSz + Su sin(^) — Sw cos{9) 4 - d 2 z, and 

Sg = —KgSg + Sq -b d2g, (4 31) 


and the augmented signal Va in this case can be derived from (4 9) as 


^ Ua ^ 
y Qa j 


^ UiZQs{Sg) - WdSixi{Sg) - KxSxCos{9) - KzSz sm{e) ^ 
Ud Sm{Sg) + Wd COS{Sg) - KxSx sm(9) + KzSz cos{9) 

^ Qd — KgSg J 


(4 32) 


Its derivative m (4 19) may be evaluated as shown below As B 2 in (4 29) for this 
problem is a function of 9 only, it is seen that for a column vector C = [Cx, Cz, Cg]^, 
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JsJ/ ® C] as required in (4 21) and (4 22) may be evaluated as 

— sin(0) cos(^) 0 

JBA^®C] = CeMiB = Ce cos{e) sm{d) 0 (4 33) 

0 0 0 

Moreover, the following result may be easily established by matrix multiplication for 

the system under study 

- = Mie ' (4 34) 

Combining (4 33) and (4 34), it is seen that 

- BJ ' Jsj/n ® C-jBji = CeM^e (4 35) 

From (4 15) and (4 16), it is seen that ^[Bj = — BJ ^7 b2[4i ® ^ Here C = X 

Thus, using Cq = 9 = Q -{- d 29 from (4 28), the above equation gives ^[BJ^] as 

\ 

^[B^-'l = -BJ ' Jb 2 [In ® XjBji = {Q + d 2 e)Mie (4 36) 

The quantity Ta m (4 21) and in (4 22) may be simplified using the equality in 
(4 35) to get 

T, = B2-i(^[B2rfV,]-K,T,) + (5Mi,(B2rfVrf-K,S,), and 

d' = -B2-'K^d2 + d2^Mifl(B2dVrf-K^S^) (4 37) 


The various terms in Ta may now be evaluated as follows 




Ud cos{Se) - WdSin(5fl) + 6dUdSva{SB) + 6dWdCos(SB) 
UdSmiSg) + W<fCos(5,) - 9dUdCos(SB) +9dWdSm{Se) 

Qd 


(4 38) 


From the dynamics m (4 31), and the definition of (4 11), one can obtain 

^ cos{9) + K^S, sm{e) - K^Sucos^{9) - K,Sus%n'^{9) 
—KxSw sin(6) cos(9) + K^Sw sin(0) cos(0) 

- Bs'^KaTa = KlSx sin(^) - iQS, cos{9) - KxSwSin^[9) - K,Swcos'^{0) 

—KxSu sm(0) cos(0) + iiT^ -S'ly sin (0) cos(0) 

^ KISb-KbSq 

(4 39) 



The third term of Tq is obtained as 


QMi0(B2dVd - K,S,) = 


Q 


-Ud sm{Sg) - Wd cos{Se) + sm{9) - K^S, cos{9) 

UdCos(Se) - ]T'dSin(5'0) - KxSxCos{9) - KzSzSm(9) 

0 


(4 40) 


The complete T^ m (4 21) may be obtained by summing up the contributions of the 
three terms given in (4 38), (4 39), and (4 40) above The Ta is seen to be a 3 x 1 
vector whose components will be represented as [Ta^,T^,T^|'^ 

Similarly the two terms in in (4 37) may be evaluated as 




f -Kxd2x cos{9) - Kzd2z sin((9) 
-Kxd2x sin{9) + K;zd2z sin{9) 
\ —Kgd2B 


, and 


(4 41) 


d2BMuiB2dYd - K^S,,) 


/- 


d2e 


Ud sm{Se) - Wd cos(Sg) + KxSx sm{9) - cos{9) 
Ud cos{Sg) - Wd sm{Sg) - cos{9) - sm{9) 

0 


\ 


(4 42) 


A conservative bound D 2 on the disturbance term d^ can be obtained, by taking the 
sum of the absolute values of the two contributing terms in (4 41) and (4 42), as 


D 


/ 

2 


Ka:D 2 x\ cos(6')| + KzD 2 z\ sin(0) | 

+T>2ait^<iSin(5fl)| + \WdCos{Sg)\ + K^\Sx\+K,\Sz\] 

K^D2x\sixii9)\+KzD2z\sm{9)\ 

+D2g [|17dCos(5e)| + IPTdSiKS'fl)! + ii:,|5,| +ir,|5,|] 

KgD2g 


(4 43) 


The time derivative of the augmented signal Va m (4 32) may be obtained using 
the above developments and Va = Tq + d^ from (4 20) With d^ and its bound D 2 
(4 43) as obtained above, the total disturbance term m S„ dynamics (4 23) is seen 
to be bounded by Di + D 2 The i = 1, 2, 3 (3 34) would now include the terms 



72 


arising due to D2 also However, the gams ki, i = 1, 2, 3 will have the same functional 
dependence on as in the Section 3 5 1 Thus the control law in (3 35) is modified 
with the new Va signals, and gives the new control law as 

i -QU-g« cos(9) - . and 

Afcom = i (4 44) 

With the above control law, the associated tracking errors S„ are guaranteed to stay 
bounded inside the boundary layers = [cjiu, 4>w, ^ discussed earlier The 

results obtained from the simulation using the parameters of the aircraft problem 
considered in Section 3 5 1 are next presented 

4.3.1 Simulation Results 

The numerical values of the unmatched uncertainties and disturbances and their bounds 
are given in Appendix D The variation in the reciprocal of aircraft mass Rm and inverse 
of the pitching moment of inertia C7 as considered for the simulation are as shown in 
Fig 4 2 The pilot maintains the equilibrium value of the control inputs till 0 5s 
and then applies an elevator doublet of magnitude ±15° and a simultaneous throttle 
doublet of ±25 units with a total time period of 1 Os The inputs are then restored to 
the equilibrium values at 1 5s and maintained thereafter The Aae uncertainty as given 
for the case UCl in the Appendix D is used in the simulation For comparison purposes, 
responses were obtained using the control law in (3 35) which was designed under the 
assumption of d2 = 0, to show the effect of the presence of these disturbances These 
responses are termed as Rq m the following discussion The responses were obtained 
for four sets of K^, z = 1, ,4 (the corresponding responses obtained are accordingly 

termed f^) to find out the effect of variation of the individual feedback gains in the 
outer loop Note that Aj’s are vectors given as [Ax, Aj, Kg]J The gam values are 
taken as Ai = [0 1, 0 1, 10 0]^, A2 = [1 0, 1 0, 1 0]^, A3 = [5 0, 5 0, 5 0]^, and 




1 95 


c^ X 10® 


1 85 


1 75 


0123456789 10 


time (s) 


Fig. 4.2 Variation in the Hm and the c-j parameters of the aircraft 


= [10, 1 0, 10 0]'^ Also, IS taken as [0 5, 0 5, 5 0 x 10“'*]'^ and r 7 „, = 10 
Numerical simulation was done using the fourth order Runge-Kutta algorithm with a 
time step of 50/iS 

With the boundary layers for S„ selected as and B 2 as defined in (4 29), a 
conservative estimate of the bounds on |Si| may be derived from (4 26) as 


|5.l < 
IS.I < 

|59| < 


(f>u + 4>w + D-ix 

K. 

(pu + 4>w + D^z j 

Y. • 

<f>Q + Die 
Ke ' 


(4.45) 


Using numerical values, the bounds on iSi, 5*, and Se may be seen to be 20 0, 20 0, 
and 0 03 respectively for the gams in case i2i, 2 0, 2 0, and 0 3 in case 0 4, 0.4, 
and 0 06 for Ra, and 2 0, 2 0, and 0 03 for R 4 , Note that Rq does not satisfy any such 
bounds 

From (4 45) it is clear that the Sx can be made arbitrarily small by an appropriate 
choice of Kx This improvement in position tracking is obtained at the cost of a 
decreased velocity tracking performance in the sense that the controller is tracking Va 
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(4 32) instead of the desired Vj Therefore, it is important to look at the difference 
between the two signals It folloivs directly from (4 27) that 


^ Ud-Ua ^ 

lid - Wa 


^ 17d[l - cos(5'0)] + WdSin(S'e) +-fCE5iCos(^) + -fir25'zSin(0) 
-C/dSin(5fl) + Wd[l - cos ( 55 )] + K:;oSxSm{9) - iir252COs(^) 


\ — Qa J 




KgSe 


(4 46) 

Although the can be made small as discussed earlier, the [Ka:Sa;]i terms as in the 
above equation are seen from (4 45) to be bounded Under the condition of small Sg 
(cos{Sg) = 1, and sin(5e) = Sg) for an appropriately large Kg, and using the bounds on 
the Sj from (4 45), — Va in the above equation can be simplified and approximate 

bounds on its variations are obtained as 


\Wi - w^\ 

^ \Qd ““ Qa\ j 


( \Wd\\Sg\+2{(})u + (j)w) + D 2 x + D 2 z 
+ (pw) + D2 x + D^z 

\ 4>q-^ D2e j 


(4 47) 


For aircraft trajectories in the low angle of attack regimes Ud » Wd, and hence the 
above equation indicates that the error Wd — Wa is expected to be larger than Ud — Ua 
Note that even if 5^ = 0 0, there will be a finite error between the desired and the 
augmented velocity signals leading to a compromised velocity tracking The magnitude 
of this error is seen from the above equation to depend on the selected boundary layer 
thicknesses and the bounds on the unmatched disturbances and uncertainties as 
given by D 2 

The tracking error with respect to the desired response of the velocity outputs 
{U, W and Q) are shown in Figs 4 3, 4 4, and 4 5 respectively for the five responses 
{Rq to R 4 ) With the inner loop control law, the velocities were seen to be tracking the 
augmented signals Va within the specified boundary layers ($„) Thus the tracking 
errors V — are presented in these figures instead of the S„ m (4 6), as the actual 
tracking errors are of direct interest It is observed from these figures that velocity 
tracking m Rq has the best response in all these figures However, as the augmented 
velocity signals are tracked for the various gam cases all the other cases show a com- 
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promised velocity tracking As discussed earlier the errors W — (~ Wa — Wd) in 

Fig 4 4 are seen to be much larger than the U — Ud Ua — Ud) signal m Fig 4 3 

Figs 4 6 , 4 7, and 4 8 show the errors in 9, Px and from their desired values 
Xd Here, as expected, it was found that Rq shows the worst response In fact the 
errors in the 6 , px, and Pz at t = 10 Os are seen to be about 0 4rad, 161ft, and 1048ft 
respectively Responses for the other cases are a lot better, and as would be expected 
the best responses are obtained for the higher feedback gams The position variable 
responses for the Rq case are not shown in the figures as errors in these variables are 
very large 

It IS seen from Fig 4 6 that increasing the gam Kg results in better 9 tracking and 
smaller Sg as is also e\udent from (4 45) As the Kg gam is increased from 1 0 m iF 2 to 
5 0 m Xa and finally to 10 0 m cases Ki and better responses are obtained Also 
seen from the figure is the fact that the responses Ri and R 4 with the same value of 
Kg = 10 0 almost overlap An increase m Kg would give smaller deviations jVd — Va| 
for the translational velocity signals (U and W) as would be expected from (4 47) 
This IS also seen from the V — plots for the two translational velocities m Figs 4 3 
and 4 4 The responses Ri and R^ m these figures are seen to be much better than the 
others and are almost overlapping However, the better 9 tracking would be achieved 
only with an associated mcrease in the Qa — Qd error Since Sq = Q — Qa is guaranteed 
to be smaller than (I)q = 5 x 10 “^, this may be seen from the Q — Qd response m Fig 
4 5 also The figure shows larger deviations from the Qd as the Kg gams are increased 
from K 2 to Kz to the maximum value of 10 0 m the gams Ki and 1^4 Here also, the 
Ri and R 4 responses are seen to be coinciding The \Qd — Qa\ 111 (4 47) is bounded by 
<j)Q + D 2 e = 0 3005 from the numerical values used All the responses m Fig 4 5 are 
seen to satisfy this bound 

Also as the Kx and Kz feedback gams are mcreased from Ki to K 2 to Kz, a better 
position tracking of the remaining vanables, 1 e the px and pz is seen from the Figs 
4 7 and 4 8 It is observed m these figures also that as the corresponding Kx and Kz 
gams are same m the gam sets K 2 and K^, the responses R 2 and R 4 are seen to be 
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Fig. 4.3 Tracking error U — U^ (fps). 

almost same, with their curves almost overlapping Thus it is seen that, the gams 
have a decoupled effect on the related position tracking error, i e a change only in 
Ke effects almost primarily the Se, similarly a change in changes the Sxy and so 
on In all cases an increase in the gains is seen to lead to better position tracking as 
was expected However, it is observed that as the gam Kq is increased to a maximum 
in the gam sets Ki and with the corresponding gams and remaining the 
same, the deviation of the actual translational velocities from their desired signals are 
minimized as is expected from the discussions above Thus, a feasible choice of the 
gams would have a large Kq, with Kx and selected possibly smaller to have an 
acceptably large Sx and Sz position errors The high value of Ke would lead to a 
lesser deviation m the generated Ua and Wa signals from their desired values and still 
achieve an acceptably small 3$ error, by marginally higher errors between the Qa and 
Qd signals, while allowing for slightly larger errors m the p* and p^ 

In all the cases the bounds on the position errors estimated using the (4 45) are 
seen to be satisfied m the Figs 4 6, 4.7 and 4 8 Thus, it is seen that with the 
outer loop feedback, better position tracking may be obtained at the cost of velocity 
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Fig. 4.6 The position tracking error Sg =■ 0 - 6^ (rad) Responses and 
i ?4 are overlapping 


tracking accuracy This then sets up a tradeoff for the control law designer, and the 
actual gain values may be appropriately selected for the various operational regimes 
and requirements 

The commanded forces and moments along with their actually achieved values m 
the i ?4 case are shown in Figs 4 9, 4 10, and 4 11 respectively Thus, the results of 
the numerical simulation show that the error bounds, as calculated using (4 45), are 
satisfied and validate the the SMC control laws in (4 44) For the general mechanical 
systems and the aircraft problem, as considered in Sections 4 2 and 4 3, with velocities 
as the primary outputs, the method has essentially been applied to systems with a 
strong relative degree of 1 The method is next extended for systems having higher 
relative degrees 


4.4 ML Analysis for Higher Order SISO Systems 

In this section, the ML strategy is applied to a higher order single input (SI) system with 
unmatched uncertainties The controller design for the unmatched uncertainty case 
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Fig. 4.7 The position tracking error Sx = Px — Pxi (ft) Responses R 2 and 
i ?4 overlap 
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Fig. 4.8 The position tracking error — Pz^ (ft) Note that and 

i ?4 are overlapping 
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Fig. 4.11 The command and achieved M 

has been considered by many researchers Spurgeon [51] and Spurgeon and Davies [52] 
however fail to address the effect of the disturbances and uncertainties on the ZD of the 
system They also require an assumption on the decomposability of uncertainty into 
matched and unmatched portions Hui and Zak [70] using similar assumptions show 
the practical stability of the system using both state and estimator feedback Under 
the above condition on the uncertainties Liao et al. [53] design continuous control laws 
to guarantee global stability and uniform boundedness of the tracking error The effect 
of the disturbances on the ZD is considered However, the verification of some of the 
assumptions on the disturbances is not trivial Qu [71] introduced equivalently matched 
uncertainties which satisfy some Lyapunov related bounds and designed continuous 
control laws to achieve global asymptotic stabilization Qu [72] relaxed the standard 
matching conditions to generalized matching conditions and a recursive procedure to 
design a continuous control law was proposed Cong and Landers [73] design control 
laws under some conditions to be satisfied by the input distribution matrix The effect 
on the ZD is not considered though Li et al [68], using singular perturbation methods 
(SPMs) design a high gam control law using stable combinations of variables, which 
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they call the stable combined variable perturbation method (SCVPM), under conditions 
to be met by the disturbance terms affecting the ZD They relax the requirement of the 
decomposability assumption Li et al [67] generalize the results of their earlier work 
to give more results about the types of Lyapunov stability achievable under varying 
conditions Wu and Chou [74] use a time varying singular perturbation method In 
the next section, a ML framework is applied to higher order nonlinear systems and 
a SMC law is developed It is shown that the method degenerates to a conventional 
SMC design, and bounds on the tracking errors in the various loops are obtained The 
proposed method differs m the selection of the gam of the switching term 

4.4.1 ML Analysis 

The system under consideration is given as [ 68 ] 

X = f(x) + g(x)u + Af(x) + Ag(x)u, and 
y = /i(x), (4 48) 

where, x 6 and u,y €.71 The functions in the above equations are smooth vector 
fields with appropriate dimensions The arguments of the functions are dropped for 
brevity The nominal system is the system defined above with the uncertain terms 
Af, and Agu = 0 The relative degree of the output, or equivalently the control 
characteristic index is assumed to be p Under a suitable coordinate transformation 
from X to z = where ^ = [yi,y 2 ) € TIP, and 77 € T?.""'’, the nominal 

system in (4 48) may be expressed in a reduced form as 

yt = yt+i, 2 = 1,2, p-1 

yp = a(^, 77 ) + 6 (^, 77 )u, and 

V = (4 49) 

with yi as the scalar output The functions 0 , b, q are as defined in standard text- 
books [21] b IS assumed to be invertible for all 2 : In the presence of uncertainties in 
(4 48), some uncertain terms will be present in the above set of equations Li et al [67, 
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68] require the knowledge of the characteristic indices of the uncertainties in the actual 
system as their design procedure is based on it For the following analysis uncertain 
terms are assumed present m all the equations, and the actual system is then given as 

Vi ~ J/i+l “h * 1)2, 1 P 1 

Up = a{z) + h{z)u + dp+ I3{z)u, and 
7] = qiz) +qi(z,u), (4 50) 

where, di, t — I, , p are taken to represent the disturbances as well as the functional 
uncertainties This may also include the effect of Agu Their bounds are assumed to 
be known and are given as |dt| < Similarly \I3\ < B q and qi are assumed to 
satisfy conditions as specified m [67, 68] Specifically, it is assumed that the system 
IS minimum phase, q is Lipschitz m ^ and the effect of the disturbances on the ZD 
(gi) IS bounded Note that Li et al [67, 68] assume that d, = 0, i = 1, ,(T, where 

cr > 1 represents the characteristic index of the uncertainty A/ Their control design 
procedure depends on the value of a In the method proposed here, this assumption 
IS not required Since all other assumptions on the ZD are taken as in [67, 68], the 
corresponding results hold true here also Therefore in the following analysis, the ZD 
has not been considered Thus only the first p equations of the system dynamics (4 50) 
are used 

The desired tracking signal is assumed to be given as = [yid, y 2 d, , ypd\^ Fol- 
lowing the ML structure as shown in Fig 4 1, augmented signals yta and its time 
derivative y^a, * = 1, ,p will be defined as was done in the Section 4 2 Convention- 
ally the errors are defined as 


yi i/zd) ^ I) ) p 


(4 51 ) 


New tracking errors with respect to the augmented signals are defined as 


Ct — 2/t 2/ta) * — I) ) P 


(4 52 ) 


The analysis for the case of p = 3 is presented here It can however be similarly 
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extended for p > 3 Prom (4 50), the ^ dynamics m this case becomes 


yi 

— y2 + di. 


(4 53) 

2/2 

= 2/3 + d 2 , and 


(4 54) 

2/3 

= a{z) + b{z)u + dz+ P{z)u 

(4 55) 

Define y^ = yu and yu = yu 

— y 2 d 



Consider the dynamics in 

(4 53) The error is ei = 

- yi — yia Note that 

ei = ei 

Taking its time derivative Ci = 

■■y 2 + di-y 2 d Using y 2 = 

2 / 2 a + e 2 from (4 52) it 

becomes 


61 = 2 / 2 a + 62 + di — y 2 d 


(4 56) 


The command signal y 2 o for the inner loop is selected, following the method of selection 
of Va (4 9) as in Section 4 2, to cancel out the known terms and introducing a stabilizing 
feedback term as 

y2a = y2d - Kieu (4 57) 

where Ki is a positive real number With j/aa as chosen above, ei m (4 56) becomes 

&\ = — + 62+^1 (4 58) 

The time derivative of y-za niay be derived using the above equation as 

y 2 a = yzd-Kiei, or 

= yzd + Kle,-Kie2-K,d, (4 59) 

The above equation may be rewritten as 

y2a = 2/2 4- d2c (4 60) 

where, d 2 c = —Kidi represents the disturbance terms and ^2 includes all the remaining 
terms in (4 59) The inner loop will be designed so that yaa and its time derivative are 
tracked 
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Next consider the dynamics in (4 54) The error is defined as 62 — y 2 ~ y 2 a 
Using ys = yza + 63 from (4 52), the time derivative of 62 is obtained using (4 54) and 
(4 60) as 

^2 = yza + 63 + <^2 ~ J/2a 

= yZa + 63 + <^2 “ ^2 — '^2c (4 61 ) 

The command signal y^a for the next inner loop is selected as was done m the case of 
yza as 

Vza = 2/2 - ^262, or 

= yzd + Klei-{K,+K 2 )e 2 , (4.62) 

where K 2 is a positive real number With yza as chosen above, 62 in (4 61) becomes 

62 = —K 2&2 + 63 + (i2 + .^i^^i (4 63) 

The time derivative of yza may be derived using the above equation and (4 58) as 

yza = y^d + Klei - (Ki + K 2 )e 2 , or 

yza = yAd-Klei + (Kl + K,K 2 + Kl)e 2 -iKi + K 2 )ez 

~{Ki + K 2 )d 2 — KiK 2 d\ (4 64) 

The above equation is rewritten as 

yza = 2/3 + dzc (4.65) 

where, dzc denotes all di, ^2 terms and the remaining terms represented by 1/3 Using 
y 2 a from (4 57) and the definition of Cj from (4 51), 62 may be rewritten as 

02 = 2/2 - 2/2d + or 

= 62 + KxEx (4 66) 

The yza (4 62) and its time derivative in (4 64) should be tracked by the next inner 
loop For * = 3, the yz dynamics (4 55) is considered next Here 63 is defined as 
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63 = ?/3 ”■ 2/3(1 Its time derivative yz — yza is obtained using (4 55) and (4 64) as 

63 = a{z) + b{z)u dz P{z)u - yza, ox 

= o + bu + 6^3 + ^v, — 1/3 — dzc (4 67) 

From (4 64), it is seen that dz — dzc — dz {^Kx ■^2)*^2 d" KxK^dx — dl-^ Let the 
bound on dg be represented as D3 Note that 6 !^ includes the effect of disturbances 
and uncertainties in the outer loops, le di and d^ Also, let \ldb~^\ < a < 1 This 
condition guarantees that 6 + /3 is invertible over the entire 2; space Then a SMC law 
may be designed as m Section 2 2 as 

“ = --1^363 -A4Sgn[e3]], (4 68) 

where As, ^4 > 0 can be appropriately selected With u as selected above, 63 in (4 67) 
becomes 


63 = -Kztz -{1+ 13b ^)iF4Sgn[e3] + ^b ^{y'^- a- KzBz) + dg (4 69) 

Following the developments in Section 2 2, K 4 is selected as 

> *^1^3 - Q - KzCzl +Dz+r]K 
~ 1 — a 

where rjK > 0 With this choice, the 63 dynamics in (4 69) can be guaranteed (using 
developments in Appendix A) to have a decay response better than 

63 = —KzCz — 77/cSgn[e3] (4 71) 

The choice of u m (4 68) is thus seen to give a proportional plus constant rate reachmg 
law and the time to reach 63 = 0 starting from a non-zero initial condition 63(0) can 
be evaluated as in (A 7) 

Using yza from (4 62) and 62 from (4 66) and noting that ei = ei, 63 may be 
simplified as 

(4 72 ) 


63 = 63 + {Ki + ^2)62 + K1K2O1 
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This IS the conventionally chosen SS with the poles at —Ki, —K 2 The newly defined 
tracking errors Cj are obtained as in (4 66) and (4 72) These are rewritten collectively 
to get 


Cl — ei, 

^2 = e 2 -{-Kiei, and 

63 = e3 + iKi+K 2 )e 2 + KiK 2 ei (4 73) 

The Cj’s are seen to be stable combinations of the conventional errors Note that [68] 
also uses stable combinations of variables with SPMs The Ci dynamics here is seen 
from (4 58), (4 63), and (4 69) as 

Cl = —KiCi + 62 + d\, 

62 = -K2e2 + es + d2 + Kidi, diJid 

63 = -i^ses - (1 + Pb~^)K4,sgn.[e3] + Pb~^(y'^ -a- K^ez) + dg, (4 74) 

where the 63 dynamics is shown to have a decay characteristic better than (4 71) 
Using yg as defined in (4 64) and (4 65), and as defined in (4 73), the yg — Kze^ 
term in u in (4 68) may be simplified to get 

y' - Kzez = yu - KiK 2 K 3 e^ - {K 1 K 2 + K 2 K 3 + - [K^ + K 2 + (4 75) 

which IS seen to be the error and its derivatives feedback term m the conventional 
SMC design (as m Section 2 2) Thus it is observed that the ML approach applied 
to the unmatched uncertainty case leads to the conventional SMC law m (4 68) , while 
incorporating the effect of the unmatched disturbances and uncertainties through the 
dg = dg — dzc term in (4 67) which is bounded by Pg, which is reflected in the selection 
of the gam of the switching term in (4 70) 

Note that the analysis can be continued for ^ > 3 in a similar fashion with similar 
results Hence the analysis has been shown up to i = 3 Also the es dynamics is 
known to be better than the dynamics specified in (4 71), and hence, it is obvious that 
63 — )■ 0 as it satisfies a sliding condition To alleviate chattering, a sat[e3/ (l>\ term is 
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used instead of the sgn[e3] term, where ^ is a selected boundary lajer thickness With 
this |e3| IS ultimately bounded by <f), and once 63 reaches the boundary layer, the 62 can 
be shown to decay and be ultimately bounded with ie2| < (refer Appendix 

F) With 63 and 62 bounded as above, ei is also seen to be ultimately bounded with 
kil < Thus it IS seen that 

kal 

|e2| 

kil 

From the definition of the Cj, i 
clear that first the SS 63 is reached, then the trajectory starts movmg along surfaces 
with lower dimensions m the e space in a sequential manner, 1 e a representative point 
would first try to reach 63 = 0, when it is sufficiently small then 62 is made small and 
so on 

The bounds for higher order systems p > 3 can be similarly denved From the 
above analysis, it is seen that if the conventional SMC design is performed, with Ki as 
selected in (4 70), then the output tracking error ei = ei will be ultimately bounded as 
in (4 76) foi the unmatched uncertainty case Given positive real numbers Ki and K 2 , 
there is a unique surface 63 = 0 (4 72) which satisfies the slidmg condition, however 
conversely, given a surface 63 = 0 with poles at —Ki and —K 2 , the correspondmg 
dynamics (4 74) cannot be uniquely defined as it depends on the order of choice of Ki 
and K 2 The output tracking error bound however remains the same The choice of 
Ki and K 2 can be done so as to keep the |ei| under some desired value For example, 
m the case of p = 3, it is seen from (4 76) that a smaller |ei| will be achieved when 
Ki and K 2 are sufficiently large and K 1 +K 2 < K 1 K 2 The area above the curve in 
Fig 4 12 gives the region of the Ki - K 2 space in which these conditions are satisfied 
Similar results may be obtained for the higher order systems too 

Under the assumptions made on the t] dynamics in (4 50), and the fact that the e, 
errors are ultimately bounded as in (4 76), the ultimate boundedness of the states and 


< <!>, 

(f> D2 + KiDi 
< ^ and 


< 


K 2 

0 + D 2 d" {Ki d~ K2^Di 
^ 1^2 


(4 76) 


1, 2, 3 in (4 73) and its dynamics in (4 74), it becomes 
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Fig 4.12 The Ki - K 2 space 

thus the practical stability [70] of the complete system (4 50) may be easily established 
using Lyapunov method as is reported in [21, 67, 68] Thus it is concluded that, the 
conventional SMC design can be done for the general uncertainty case with a new 
method of switching gain selection as in (4 70), and the tracking error bounds can be 
derived as abo\e The method developed here is applied to an example problem in the 
next section 
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^ 0 5exp(xi) 


( 0 

A/(a:) = 

0 5 exp ( 0 : 3 ) 

, £\g{x) = 

0 5 sm(a; 2 ) 


0 5exp(x4) 

0 5003(2:3) 


1 0 ) 


1 0 J 


The output y = Xi The above set of equations may be rewritten in a reduced form in 
the new set of coordinates [67] given by ^ = [2:4, a^i +0:4, Xi+x^ and r] = X 2 -xl 
as in (4 50) The control characteristic index of the nominal system is seen to be p = 3 
The resulting disturbance vector is obtained as 


d = 


^ 0 ^ 
0 5exp(a:i) 

^05 exp(a;i) +05 exp(2;4) +05 cos(a;3)u j 


(4 79) 


The SS was selected to have both of its poles at -10, 1 e 63 = 63 + 20e2 + lOOei 
Also (j), Kz, and rjK were taken as 0 05, 5, and 1 0 respectively The initial condition 
was X = [2 0, 2 0, 2 0, 2 O]"^ and the desired output was taken to be sin(7rt) It is 
easily established that a < 1 for this example 63(0) is calculated to be 223 17 With 
the choices of and tjk as above the reaching time may be evaluated from (A 8) as 
1 359s 


Simulation was done using the fourth order Runge-Kutta algorithm with a time 
varying step size This was necessary as the system equations were found to be ill 
conditioned and application of a constant step Runge-Kutta algorithm required a pro- 
hibitively small time step These issues and an adaptive step selection procedure are 
discussed in Zhao and Utkin [75] In the present work, a simpler method involving the 
doubling and halving of step sizes was used The 63 is seen to go mside the selected 
boundary layer m a very small time, approximately 0 076s which is well smaller than 
the estimated reaching time, and then stays there Note that di = 0 from (4 79) Since 
4”! = 2:4 should track sin(7ri), and ^2 = +^4 should track tt cos(7rt), an estimate of the 

bound on the magnitude of aii = tt cos(7rt) — sin(7rt) was taken as tt + 1 and is used to 
evaluate the bound on in (4 79) D 2 — 31.45 approximately Once sliding begins, 
the error 62 is expected to decay and be bounded as specified in (4 76) by 3 15 This 
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Fig. 4.13 Output tracking error ei(= Ex) 

IS achieved by t = 0 23283s, following which the ex should decay and remain bounded 
by 0 315 This bound is seen to be satisfied for all t > 0 38127s Fig 4 13 shows the 
time variation of the conventional error signal Ex It is seen that ex decays and stays 
within Its bounds. The variation of 62 signal defined in (4 73) is shown in Fig 4 14 
Thus, all the error signals are seen to be varying between specified bounds 


4.5 Conclusions 

In this chapter, robust nonlinear control law design using a ML strategy for the un- 
matched uncertainty case is studied The outer loop tracking errors are used to augment 
the command signals for the inner loop dynamics A SMC law is used m the inner 
loop to track the augmented signals The method is directly applicable to almost all 
mechanical systems with velocities as the primary outputs The ultimate boundedness 
of the tracking errors in the two loops is established Control laws are derived for the 
aircraft longitudinal dynamics and simulation results for the velocity tracking mode 
are presented. The choice of parameters may be done based on a tradeoff between the 
errors in the inner and outer loops The ML method is then applied to higher order 





Fig. 4.14 Output tracking error ea 

single input systems (with relative degree > 1) with unmatched uncertainties The 
ML analysis is seen to give control laws equivalent to the conventional SMC design, 
however, the selection of the switching term gam is slightly modified Bounds on the 
defined tracking errors are developed Thus, it is concluded that a conventional SMC 
design can be implemented for the case of unmatched uncertainty also and the bounds 
for the tracking errors may be obtained Minimization of these bounds and the method 
proposed in Chapter 2 can guide the selection of the SMC parameters The method is 
verified by applying it to an example problem taken from published literature. 




Chapter 5 

SMC for Systems with Input 
Dynamics 


5.1 Introduction 

As discussed earlier, the SMC uses a discontinuous/switchmg control law, because of 
the sgn[ ] teim in (2 9), and the associated invariance to matched disturbances and 
uncertainties accrue because of the assumption of exact switching The SMC method 
can therefore not be applied to systems where discontinuities in the input are not 
permissible, as in chemical systems where the inputs such as flow rates cannot be 
varied discontinuously Also, any practical system will have its inherent delays and a 
finite bandwidth, which rules out the possibility of an infinite switching frequency in 
practice Thus it is not possible to ensure that motion will be restricted to the SS 
(ideal sliding), however motion close to the SS (real sliding) may be guaranteed [54, 
55] This was done by using continuous approximations to the discontinuous control 
laws [56] The use of a boundary layer [21] as m the earlier chapters is one such method 
and the motion is guaranteed to he within a boundary layer of the SS An adaptive 
method of alleviation of chattering was proposed by Chang et al [76] 

Elmali and Olgac [77] used a first order filter to smoothen out the response of the 
switching function (s) and then used the signum of the filter output to design SMC 
laws The method using a second order sliding mode results m less frequent switchmgs 
than the conventional method Emel’yanov et al [54] and Levant [55] used higher 
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order sliding modes for this problem Chiacchiarmi et al [60] using a similar approach 
obtained a continuous control action which is the output of a senes of integrators 
Ramirez [57, 58] using a differential algebraic approach, expressed the system in a 
generalized controller canonical form involving the inputs, the states and their higher 
order derivatives, and synthesized SMC laws with applications to nonlinear chemical 
processes [57] Ramirez et al [59] augmented the system equations with the input 
dynamic equations and represented the extended system in a canonical form This 
method was applied to the vertical flight regulation of a nonlinear helicopter model 
Most of these works require higher order derivatives of the outputs or the a:-dynamics 
(5 2), which involves derivatives of the uncertainty or disturbance terms present therein 
Thus, the analysis of robustness properties becomes rather involved and has been 
addressed in very few works [60] 

In this chapter, the ML strategy as outlined in the previous chapter, is applied to 
this problem As seen earlier in Section 4 2, the effect of disturbances m the outer 
loop is easily passed on to the inner loop, and thus the propagation of disturbances 
and as a result the robustness studies are more easily done in this procedure, as will 
be shown m this chapter The method however requires higher order derivatives of the 
desired response signals The inputs applied to any system are generally outputs of 
an actuator system with a finite bandwidth This is modeled in the present work as 
a first order dynamic system This dynamic system is added to the original dynamic 
system forming an augmented state space Following the developments in Chapter 4, 
the tracking of position variables is done m the outer loop by generating augmented 
velocity command signals, which are tracked in the inner loop by synthesizing aug- 
mented input signals, which are tracked in the innermost loop The resulting control 
law uses the full augmented state feedback, and has a discontinuity m the first denva- 
tive of the input It thus gives a continuous response m the system state space and 
the discontinuity is limited to the controller state space Control law is sjmthesized 
for a general SI mechanical system and applied to the nonlinear spring mass problem, 
considered earlier in Section 2 4 1, and simulation results are presented Time domain 
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bounds for the state and input errors are developed The method is then extended to 
the MI aircraft longitudinal dynamics problem of Section 4 3 Bounds on the position, 
velocity, and input errors are derived Simulation results using the developed control 
laws are presented The effect of the feedback gams is studied qualitatively 

5.2 ML Analysis for a SI Mechanical System 

The specified first order actuator dynamics are added to the general SI mechanical 
system in (2 1) The complete system, may then be expressed as 

u = fi{x,x,u) + bi{x,x,u)w-\-du, &nd (5 1) 

X = f 2 {x,x) +b 2 {x,x)u + d:„ (5 2) 

where, w represents the actuator input Let the state vector be represented as x = 
[xi, 3 : 2 ]"’' = [a, As a general case the input dynamics have been taken to be 

functions of the system states (x) also The d^, d^ denote the disturbances and uncer- 
tainties in the channels specified by their subscripts and are known to be bounded by 
Dx, Du respectively Let = 6i(l-l-55i), where represents the nominal input function 
and |55i| < A, z = 1,2 Moreover, bi and 62 are assumed to be invertible Also x is 
supposed to track Xi^, X 2 d = Xu and xu = X 2 d The above set of system equations is 
seen to be similar in structure to (4 1 ) and (4 2 ) The ML scheme as shown m Fig 4 1 

IS applied to this problem In the outer loop the x dynamics (5 2 ) is used to generate 

a signal Ua and Ua which should be tracked by the inner loop, which m this case is the 
u dynamics (5 1) 

Define e = x - xu and e = x - xu Appropriate error functions m a; and 11 are 
defined as in (2 5) and (4 6 ) as 

Sx — e + Ae, and (5*3) 

Su = U — Ua (5-4) 

where, A > 0 and Ua is yet to be defined It is assumed that the input u is available 
for feedback Following the steps in Section 2 2 , Si = e -t- Ae and using (5 2 ), it may be 
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rewritten as 


— /2 + b 2 U + rfx — Xzd + Ae, or 

= /2 + hUa + h2Su + 4 - Xu + Ae, (5 5) 

using u Ug, + Sii from (5 4) As in Section 4 2, Ug would be required for the inner 
loop control law design, the control signal Ug is required to be continuous Thus, it is 
selected by replacing the signum term in (2 7) by a term as in (4 9) and is of the 
form 

Ua = b2^{x3d- f2-ae-be-kxSx), (5 6) 

where, a,b > 0 are selected as in Section 2 3 with the poles selected as -A, -e Thus, 
a = Xe and b = X + e from (2 8) Thus, the —ae — be — k^s^ term m the above equation 
may be written as —{kx + e)sx — Ae Using Ug, as chosen in the abo\e equation, Sx in 
(5 5) becomes 


Sx = 5b2{x3d - 12 + Xe) - (1 + Sb 2 )ikx + e)sx + 625u + 4 (5.7) 


It IS seen that from the above equation that calculation of bounds for Sx requires 
knowledge about the s„ bounds, and that is developed next 

The inner loop control law is to be designed so that u tracks Ua,ie,Sg—^0 As in 
the Section 4 2, the time derivative of Ug in (5 6) is required It may be denved as 


di 

= -^[b2^]{x3d - f2 


d 

ae- be- + ^2 ^ - f 2 -ae- be- fcx^x] (5 8) 


From the above equation it is seen that a higher order derivative of xsd would be 
required The remaining terms, 1 e b 2 , f 2 , Sx and their derivatives are functions of the 
state \ector x and x Accelerometers in mechanical systems are quite commonly used 
Hence the assumption of x being available for feedback is quite justified and will be 
used in the following developments However, this is not a necessary condition. In case 
the acceleration signals are not available for feedback, the above derivation of Ug would 
have a disturbance term as was the case in Section 4 2 (4 20) This disturbance term 
would then be passed on to the inner loop as was done in the previous chapter The 
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control law design for both the cases is outlined in the example problem considered 
in the next section For the rest of this section, it is assumed that the accelerometer 
signals arc available for feedback 

Since xu — o-^ — be-kxSx — 2:3^ - ae - 6e - - /ciAe = xzd - a)e - {kx-\~b)e, 

its time deiivative may be easily obtained as 


dt 


[xsc/ - ae- be - - (^i + b){x- xzd) - {k^X + a)e 


(5 9) 


Using (E 3) and the notations in Appendix E, the derivative ^[^2^] 8) may be 

do 

d 

expressed as -b^^Rb^x Similarly ^^.y be expressed as R/^x Using these time 

CLTf 

deiivatives and (5 9), Ua in (5 8) may be obtained as 


Ua = -b^'^i^zd - f2-ae-be~ kxSx)Rb2^ 

+b2^[xid - (kx + b){x- xzd) - {kxX + a)e- %x], or 

Ua = b2^[x4d ~ + b)(x — Xzd) - {kxX + CL)e 

—Rf^x — 62 ~ ue — be — kxSx)Rb2^] 

Thus with the signals Ua and Ua synthesized in the outer loop as defined in (5 6) 
and (5 10), the design of the inner loop control law is done so as to make Su m (5 4) 
go to 0 Taking the time derivative of s^, it follows from (5 1) that 

Su = fi + hw + du-Ua (5 11) 

The input w is selected as was done m Section 2 2 to satisfy a constant plus proportional 
rate reaching law as 

w = b^^ (ua - A - k^Su - A:dSgn[s„]) (5 12) 

wliori', ha, kd>0 The control law thus requires that all the states x and the input u 
be available for feedback With w as selected above, the s„ m (5 11) becomes 

s„ = 6bi{ua - A + Ksu) - kuSu - (1 + Sbi)kdSEn[su] + 4 (5 13) 


If kd IS selected following (4 70) as 

, PllUa 

kd = 


A + k gSul + Du + Vu 

YZJ, ’ 


(5 14) 
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«ith I,,. > 0, then following developments m Section 2 2, s. can be shown to have decay 
chaiactciistics better than 

Su = —kuSu — 77aSgn[s„] jgj 

Thus 0 in a finite time uhich may be derived using (A 7) and will then stay 

theie Fen implementation purposes, the sgn[] term is replaced b, a sat[s„/^J as m 
Section 3 4, whoie (J)^ is an appropriately selected boundary lajer thickness Note that 
this lca\es w as a continuous signal with discontinuities in its first derivative only 
Thus, using a continuous u; makes the actually applied u continuous up to its second 
del native' In this case, will be ultimately bounded with 


^ 4^11 


(5 16) 


within a finite time which may be evaluated as in (A 8) 

I hen th(’ tenn in (o /) will be bounded by \h 2 \ 4 u With the further assumption 
that 5b ui ~ f'i in (5 7) is bounded as P2\I^B ^ 3 d ~ /a + Ae|, it can be shown 
as in Appi'iuhx F that sv will be ultimately bounded with 


s- I < ~ /a + Ae)| + |b2|<^u + -Pj 

(1 — /?2)(^i + e) 


(5 17) 


Thus, the I'lrois di'linecl in (5 3) and (5 4) are ultimately bounded as shown above 
Ihe (h'signed control law is applied to the nonlinear spring mass problem considered 
in Section 2 4 1. 


5.2.1 A Nonlinear Spring Mass Damper System 

The system equations aie as given in Section 2 41 The u dynamics is assumed to be 
a fust onl(*r actuator system with a bandwidth of lOrad/s, 1 e u = — lOu + lOio + du 
d,t IS taki'u as 0 15 + 0 09cos(37rt) bounded hy Bu = 0 25 The controller parameters 
weie taken as follows* A = 6, € = 2, and kx, k^ = 10 4u = 0 1 and = 1 0 were used 
m the simulation. With this a — 12 and 6 = 8 As 62 = 1 and 02 = 0, Ua in (5 6) 
becomes 

Ua=xsd-f2~0e-12sx (5 18) 
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Ua IS next calculated from (5 10) as 

Ua = — 18e — 72e — Rf^x 

= - 18(a; - Sad) - 72(s - X2d) - Rf^x (5 19) 

Rf,x can be evaluated fiom the functions given in Section 2 4 1 as (pio - Z^ix^)x - 
{I'o + 2!^i 1 1 |sgn[s]) r Thus, becomes 

Ua = t w + 18 cjd + 72x2d - {fiQ - Zfiix'^ + 72)x + (i/q + 2z/i|s|sgn[s] - 18)s (5 20) 

Note that if the acceleiometei is not available for feedback, the x term m the above 
eciuatum mav instead bo wiitten using (5 2) as s = /a + b 2 U + b 2 Sb 2 U + 4 In this the 
hist, two teuiis arc’ exactly known and the remaining two terms being uncertain will be 
passed onto the inner loop as a disturbance term (as m Section 4 2) The Ua may then 
h(‘ expu'ssed as 


Ua ~ [2<'4d + ISlCsd + 72X2d ~ (jUfl ~ + 72)x 

+(1^0 + 2i,'i|xlsgn[x] - 18) (/a + 62 w)] 

+(1^0 + 2i/i|x|sgn[x] - 18)(62^62U + 4) (5 21) 

Th(> Ua as (U'lived m (5 20) or (5 21) m the presence and absence of accelerometers 
iesp('ctiv('ly may be lewiitten as Tu + d'^, where is the term m the above equations 
not Inuing any unceitamties or disturbances, and the contains the effect of the 
distui bailees’ anti unetutamties Thus it is seen that for the case of acceleration signals 
being jnesent, dy = 0, arid the right hand side of (5 20) gives For the other case, 
dy IS seen to be. the last term in (5 21). It is assumed that a bound on may be 
estimated as £>;. Thus for the case where the accelerometers are absent, Sy m (5.11) 
becomes 

Su = fi + - d'J - Tu (5-22) 

Tilt' input w for this may be designed as m (5 12) Note that this approach requires 
higher switching gains, because of the presence of two disturbance terms This is 
expected to give rise to more chattering effects 




Fig. 5.1 The tracking responses for Cases % and a\ 

Simulation was done using the fourth order Runge-Kutta algorithm with a time step 
of 1ms In the following figures, the responses labeled oq are obtained in the case of 
accelerometers being not used The other case is labeled Fig 5 1 shows the tracking 
response obtained for the selected tracking signal A good tracking performance is seen 
in both the cases The bound on the tracking error as given m (5 17) is seen to depend 
on the functional like ( 3 : 3 d - h + Ae), and thus numerical values of these bounds are 
difficult to be evaluated The figure qualitatively shows that these are indeed bounded 
and quite small. Figs 5.2 and 5 3 show the u dynamics (5 1) with the signal Ua to be 
tracked for cases Gq and oi respectively It is observed that u in both the cases starts 
tracking u. within about 0 5s. In both cases the input error was seen to lie inside the 
specified boundary layer. It was however observed from the two figures that the input 
tracking in the case ao is slightly faster than the m response This may be attributed 
to the fact that higher switching gam m this case leads to better Ua tracking Thus it 
IS seen that the ML strategy can be used to design continuous SMC laws In the next 
section, the method developed for SI systems is extended to MI systems 
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5.3 Aircraft Longitudinal Control with First Or- 
der Input Dynamics 


The method applied to a SI system m the previous section is extended to a MI aircraft 
control problem considered m Section 4 3 In this section it is assumed that accelerom- 
eters to measure the V vector are present, i e U, W, Q are available for feedback It 
IS also assumed that the dynamics of the force vector is decoupled Augmenting the 
input dynamics equations to the MI aircraft system of equations in (4 1) and (4 3), the 
complete set of governing equations becomes 

u = g(u,X,V)+Bo(u,X,V)w + do, (5 23) 

V = fi(X,V)+Bi(X,V)u + di, and (5 24) 

X = B2(X)V + d2 (5 25) 


Hoic w € 7^" represents the actuator input, and the matrices B„ i = 0,1,2 are 
assumed to be invertible over the domain of interest, and do represents the disturbances 
and uncertainties in the input dynamics and is given to be bounded by Do Let 
2 = The other assumptions/bounds in Section 4 2 are taken to be valid 

Pi Decoding as in Section 4 2, the error functions for the position, velocity and input 
vaiiablcs are defined as Sj, S,, and Si, as 

Sa; = X Xrf, 
s„ = V-Va, and 
S„ = U-Ua 


(5 26) 
(5 27) 
(5 28) 


wlierc, V„ u, are yet to be defined The procedure .s developed using the ML strategy 
as outlmed m Section 5 2 The control law for an inner loop is therefore designed 
based on the error dynamics of the immediate outer loop Thus S, dynamics is set by 
synthesizing V. and its time derivative V„ for tracking m the inner loop Smularly 
to obtain a desirable S. dynamics, and u. signals wiU be generated as shorn m 
the sequel. The actuator input w for the input dynamics (5 23) is designed so t at u 

tracks Ua in the inner most loop 
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Following the developments m Section 4 2, may be derived as m (4 9) For this 
choice, Sj, has been earlier shown to satisfy the dynamics in (4 10) as 

= -K*S^ + B2S„ + da, (5 29) 

V(i IS lequiicd for tracking in the inner loop as m Section 4 2 and is derived as in (4 19) 
The Sy dynamics can be obtained using (5 27) and (5 24) as 

S„ = V-V„, or 

= fi+Biu + di-V, (5 30) 


Fiom (5 28), u = + S^, and using Va = Ta + d'j from (4 20), the above equation 

nia\ be rewiitten as 

S„ = fi +BiUa + BiSu + di - Ta - d2 (5 31) 


The coniniand signal Ua should be selected so as to make the Su errors small Therefore, 
It IS seh'cted to cancel out the known terms and introduce a stabilizing feedback as m 
Si'ctioii 3 4 Since Ua would be required for tracking m the mner loop, the discontinuous 
bgn[ ] teiiii in (3 13) is replaced by a K„S„ term as m (4 9) Ua is thus selected as 


K. = diag[A'„. Kw, A'gl. all of which are positive real aumbers B, represents 
the nominal B, matrix as defined m Section 3 4 Note that u, represents the desired 
command forces which is different from the actually applied forces as discussed m 
Section 3 3 The effect of this uncertainty is incorporated mto the design procedure 
(by defining the Ai matrix m 3 15) as was done m Section 3 4 Usmg u. as chosen m 
(5 32), S„ m (6 31) may be obtained as m (3 17) as 

S„ = -[/ + A,]K,S. + A,(T, - fi) + BiS„ + d, - d' (5 33) 


The S„ dynamics thus involves the BiS„ term, and thus the dynamics of S need 
to be considered before any bounds on the S. errors can be developed For this, the 
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time derivative of u, m (5 32) is next developed For the tracking mode, since Bi in 
(3 29) IS a constant matrix (Section 3 5 1), is seen to be 

d 

= Br'-[T,-fi-K,S„], or 

= (5 34) 

d 

where — [fi] = Mf,Z, using the notation given in (E3) in Appendix E Note that 
S„ = V - Va fiom (5 27) Thus, using Va from (4 20), S„ = V - - d '2 Putting 

this m the above equation, Ua may be rewritten as 

Ua = Br'[T, K„(V - T,) + K„dy (5 35) 

Now It leinains to evaluate Ta This may be derived starting from the definition of Tq 
in (4 37) Its time deruative may be evaluated as 

+4[QA/„I(BmVj - KA) + - K.SJ (5 36) 

at CL^ 

— [Bj^] IS alieady derived in (4 36) Also S^; = T^; + 62 from (4 11) In the above 

equation J^iQ^ho] may be evaluated as QMig + QJu^ei^n ® X] From Mig as defined 
in (4 33), it IS easily obtained for the column vector C as used m (4 33), that 

- cos(0) - sin(^) 0 

Jmis [I ®C] = CgM 2 g = Cg - sin(0) cos(6») 0 (5 37) 

0 0 0 

Thus, ® X] becomes {Q + ^ 20)^25 as C here is X and = 0 = g + ^20 

from (4 28) Thus, ~[QMig] = QM^g + + {d 2 gQ)M 2 e Using these equalities in 

di 

(5 36), Ta may be obtained as 

T„ = (20iW„ + <i„M.,)(|[B,M-K,T.)+B,-'(^[B,,V,l-K.TJ 

+(QMu + + (d 2 sQ)M 2 ,)(B 2 iVi - K,S,) - QAfijK.d 2 . or(6 38) 
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d\ KajTj;) + B2 ^(^[B2(iVd] — KjTj;) 

HQMie + Q"M29)(B2dVd - K*S^) - QM^gKA^ 

+d26Mu{~[B2dVd] - K:,T^) + {d20Q)M2g{B2dYd - K^SJ (5 39) 

Also needed in the above equation is -K^T^ Prom its definition T^; = -K^Si + 
B 2 SU m (4 11), Ta; = — KiSi + — [B 2 ]S^ + B 2 S„ Since Vq = Tq + d^ from (4 20), 

J 

Sy may be obtained from (5 27) as V - Also, using (4 16) and (4 33), — [B 2 ] 

dtij 

may be written as {Q + d 2 e)Mig Using this, S,, as developed, and Sj, from (4 11), T^, 
may be expressed as 


Ti — — Ka;(Ti + d 2 ) + {Q + d2B)MigSy 

+B 2 (V-T,-d') (5 40) 

Using the definition of d^ from (4 20), B 2 d 2 = -Ka;d 2 + d 2 gB 2 Mig(B 2 dVd — K^Sx) 
Using the reverse equality from (4 34), le replacing Mig by — m this 
equation B 2 d 2 becomes equal to — Ka;d 2 — d 2 eMi 0 Bj^(B 2 dVd— Ka;Sa;), which simplifies 
to — Kxd 2 — d 28 MigYa from the definition of Ya m (4 9) Using this m (5 40) and 
substituting Ya + S„ by V, as is evident from (5 27), gives 

T, = -Ka;T,, + QMigSy + B2(V - T,) + d2BMigY (5 41) 

Using Tx as above in (5 39) gives Tq as 

T. = 20Mi,{||B2dVJ-K.T,)+B,-'[A(B2<Vjl + K=T, 

-QK.MijS.] - Bj ‘K.B2(V - T.) + (QMu + Q^U2e)(B,iVi - K^S^) 

d 

-”(3Mi^Ka;d2 + d 2 eMi(}(—[B 2 dVd] ~~ 

'i^{d2oQ)^^2d{’^2ci^d ~ — <^20^2 (5 42) 
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1 ^ 

The term B 2 -^\B 2 dy d\ m the above equation may be evaluated as 


( cos{Sg)[Ud + - UiQl + - sm(5fl)[Wd - lU^Qa - WgQl - UdQd] ^ 

cos(5^)[Wd - 2UdQd - WdQl - UdQd] + sm(Sg)[Ud + 2WdQd - UdQl + WdQd] 

y Qd j 

(5 43) 

Note that the second derivative terms Ud and so on as required in the above equation are 
evaluated as in (D 9) in the Appendix D These signals are seen to include higher order 
derivatives of the pilot inputs Up Thus the so designed control law will be sensitive 
to the rate of application of pilot inputs However, this must be passed through a low 
pass filter to avoid any pilot induced oscillations 

Thus Ta may be completely evaluated using (5 43) in (5 42) Using fj as defined 
in (3 27), the term Mf^Z in Uq (5 35) may be obtained as 



^ -QW -QW- go cos(d) ^ 


^ go cos{d) 

Mj,Z = 

QU + QU- gosm(d) 

— d2g 

go sm(d) 


[ 0 J 


[ 0 ; 


= fia — d2giid (fi 44) 


Using (5 42) and the above equation m (5 35), Ua may be evaluated For simplicity, it 
IS rewritten as 

u, = T„ + d;, (5 45) 

where, the first term is completely known and the d'u contains the effect of all the 
disturbances present Thus, the term may be obtained as 

+K^Ta; — QK^MigSv) — [Bj’'Ka:B2 + K„](V — Ta) 

+{QMxg + (5'M2e)(B2dVd - K,s^) - fij (5 46) 

Similarly, the total disturbance component of Ua is evaluated from (5 35), (5 42), and 
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(5 44) as 

< = Br'[-QM»K,d2 + d2,Mi,(4[BMV,]-K,TJ 

+{<ktQ)M2e(^2i'Vi - KA) - d2«Bj-'K,M»V + djsfu + K.d;i (5 47) 

A conservative bound (D„) on this term may be found as m Section 4 3 

Now that Ufl (5 32) and its time derivative (5 45) have been evaluated, consider 
the S„ (5 28) dynamics Its time derivative may be expressed as = u — Uq, and 
using (5 23) and Ua as defined in (5 45), it is rewritten as Su = g+BoW+do — 

The input w may be selected, following the procedure in Section 3 4 to ensure that 
satisfies a constant rate reaching law, as 

w = B^^(T^j-g-Kj,sgn[Su]), (5 48) 

where K,,sgn[Su] is a n vector with its component as ii',j,sgn[6'uJ as m Section 3 4 
Kr,^ IS selected as Do, + + r]„ where all With this, the Su, response can be 

shown (as m Section 3 4) to be better than the response of 

Su, = -77tSgn[5„J (5 49) 

The sgn[ ] term m the control law (5 48) is replaced by an equivalent s^t\Su,/^u^ 
term as in (3 22) Also let represent the vector of boundary layer thicknesses as 
1 = 1, ,n Thus as in Section 3 4, the Su, is guaranteed to be inside a 

boundary layer of thickness ^^u,, i e 

IS’u.l < K p 50) 


within a finite time as given in (3 25) 

With Su shown to be bounded by $u, the bounds for the S„ and Sx errors may now 
be derived Note that the dynamics in (5 33) involves d^ which is bounded by Dj 
(4 43) This bound is seen to depend on Sx Also, the Sx dynamics as derived in (5 29) 
depends on S^ Thus, the and Sx dynamics are seen to be mter-dependent on one 
other Under the assumption that S* can be shown to be bounded, and (T^ - fi) term 


J 
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IS bounded, and noticing that Dt is a diagonal matrix as given in (3 33), it follows from 
the Sy dynamics (5 33) that the error satisfies 


K|< 


[|Bi|$. + Di + D^ + A(|T,-fi|)], 


(5 51) 


With |Su| bounded as above, the error can be shown to be bounded, as in (4 26), 
by 

|5^J < (5 52 ) 

The velocity and the position tracking errors S„ and errors are expected to be 
bounded as shown above As discussed above, the |S„| and iSij are dependent on each 
other, resulting in simultaneous linear equations as shown in (G 5) with the bounds 
as variables These bounds should be positive to be practically meaningful It is 
shown m the Appendix G that this requires that K^, should be sufficiently large, or 
more specifically Ku, Kw » 2 As The actual bounds may be derived by solving the 
linear equation set in (G 5) It may be easily shown that these errors are uniformly 
ultimately bounded as in Appendix F with respect to these bounds However, the jA, I 
depends on functionals like [jT^ which cannot be easily estimated and hence the 

numerical values of these bounds cannot be derived, and thus detailed derivations are 
not presented It suffices to know that these quantities are bounded Thus, the ML 
sti ategy is easily applicable to a MI problem as shown here Simulation results for the 
aircraft longitudinal dynamics is presented next The qualitative nature of variation 
of the error bounds for different sets of gams are also studied 


5.3.1 Simulation Results 

Simulation was done on the aircraft example considered in Section 4 3 In this example, 
the pilot maintains the equilibrium values of throttle gear and elevator till 0 5s, and 
then applies one complete cycle of the sinusoid elevator and throttle gear commands 
simultaneously The elevator command is assumed to be A A(l) = ^ei = 20 sin(2t-l 0) 
while the throttle command is taken as AUp{2) = Stg = 8sin(2t - 1 0) The pilot 
then maintains the control inputs at the equilibrium values after tt + 0 5s Note that 
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the Up commands are continuous, however, Up is not The input dynamics (5 23) 
foi each component of the force vector was taken to be a first order system as Ui = 
— 20 Uj + 20 wj + do. The disturbance vector do was taken to be 

— 10 5 + 8 5 cos(77ri), 

~ 10 5 - 8 5 cos( 67 ri), and 

du — 10 5 + 5 5 cos( 57 ri), (5 53) 

with the bound £)„ = 20 0 used for in Section 4 3 was taken as = 
[1 0, 1 0, 10 0]'^ and were taken as 1 0 and 50 0 respectively For compari- 
son purposes, the response that would be obtained by neglecting the presence of the 
input dynamics, i e by using the control law in (4 44), is obtained and is labeled 
as Ra in the following discussions and figures The responses for different choices 
of K„ = [Ku.Kw.KqY are presented next is chosen as [10 0 , 10 0 , 10 0 ]^, 

K „3 = [10 0, 30 0, 20 Of, and K ,3 = [20 0, 30 0, 20 0]^ The corresponding responses 
are termed Ru R 2 , and R 3 respectively Prom Bi as defined m (3 29) and the nominal 
parameters as given in Appendix D, it is clear that will have large diagonal entries 
leading to a large d'^, as can be seen from (5 47), and thus the gams K,, are expected 
to be quite large Thus a smaller time step for the Runge-Kutta algonthm needs to 
be taken Here it was selected as 50/iS for Rq, 5^s for Ri, and R 2 , and 2 5/rs for the 
higher gam case of K „3 

Figs 5 4, 5 5, and 5 6 present the Su = V - Va responses obtained for the various 
gam cases K„ for the velocity variables U,W, and Q respectively Note that with the 
ML strategy, the velocities V are supposed to track the Va signals as discussed m 
the previous section For the Ro response, which is the case where no input dynamics 
are present, the S„ errors should be bounded by the boundary layer thicknesses = 
[0 5, 0 5, 5 0 X 10-'^]'^ as selected m Section 4 2 That this condition is actually satisfied 
is easily seen from the Ro response in these figures The S„ responses for the gam cases 
z = 1, 2, 3 m the Figs 5 4, 5 5, and 5 6 show that, with the specified input 
dynamics, larger deviations are seen at the smaller gam values and these are seen 
to decrease with higher K„ gams as is seen from the R 2 and R 3 responses These errors 
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creep into the system because of the additional input dynamics Thus, it is observed 
from the figures that an increase of the K„ gams leads to smaller S„ as was expected 
from (5 51) 

It IS also clear from the figures that an increase in Ku leads to a smaller value of 
Su only This is seen m the U — Ua response m Fig 5 4 The Ku gam has the same 
value in the gam sets and while the other Kw and Kq gams are increased 
The U — Ua response for these cases are seen to be almost same as from responses Ri 
and i ?2 as m Fig 5 4 and are found to overlap However, as the Kw and Kq gams 
are increased from to the corresponding R 2 responses m the Figs 5 5 and 
5 6 show a marked improvement over the R\ responses, with almost no effect on the 
Su Similarly, as the Ku gam is increased from to Kyj, Rz m Fig 5 4 shows a 
distinct improvement over the and R 2 responses Here, the gams Kw and Kq are 
the same and the R 2 and Rz responses for the associated velocity variables W and Q 
show almost no improvement with these responses almost overlapping as in Figs 5 5 
and 5 6 Thus the feedback gams Ky have a decoupled effect on the related velocity 
errors It is thus observed from the figures that as the magnitudes of the gain Ky 
are increased, the system performance approaches the Rq performance m the limit, as 
would have been expected, since the Rq effectively has ’unity’ as the specified mput 
dynamics (ucom = w), implying a infinity bandwidth dynamic system 

As discussed earlier, as the Kq gams are increased, better Qa tracking is expected 
and this would result in better 9 tracking Similarly, for the higher Ku and Kw gams 
better U and W tracking of the augmented signals are observed This unplies that 
larger Ky gams should result m better position tracking performances The obtained 
position error (S. responses are shown in Figs 5 7, 5 8, and 5 9 for the considered gam 
cases A distinct improvement m the tracking performances is observed as the gams 
are increased due to the better velocity tracking performance and is seen to approach 
the Ro response The errors m the high gam cases of R 2 and Rz are seen to be 
almost overlapping m these figures The effect of the K,, gams on the S. errors has 
been discussed in Section 4 3 1 However, as discussed m the previous section, the 
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The Su = U -Ua response Note that responses Ri and R2 overlap 



Fig. 5.5 The Sw = W -Wa response The R2 and R3 are seen to be overlapping 
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Fig. 5.6 The Sq-Q-Qa response The i ?2 and are seen to be overlapping 

and S„ errors are inter-dependent (as shown in Appendix G) and the effect of varying 
the gains and K„ simultaneously needs further study 

The tracking error of the input dynamics (u — Ua) is as shown in Figs 5 10, 5 11, 
and 5 12 It is seen from these figures that the Su error (5 28) remains inside the 

selected boundary layer thicknesses (<^u, = 50 0) 

It may be concluded that as the K„ gam terms are increased the S„ error, introduced 
due to the input dynamics, decreases Thus it is seen that the input dynamics can be 
easily incorporated m a SMC approach with an attendant increase in velocity and 
position tracking errors, which can however be minimized by choosing appropriate 
gams as shown above The simulation results validate the control design method for 
the MIMO case The various loop errors are seen to be bounded though by unknown 
quantities A detailed parametric study of the involved feedback gams needs to be 

done for their proper selection 

This chapter so far considers the development of a SMC design method for the case 
of specified first order input dynamics However, the method is not limited to this class 
of systems only and may be easily extended to the case of higher order input dynamics 



Fig. 5.8- The Pz - Pz, responses Note that the J?2 and R 3 responses overlap 
















116 


using the ML analysis as was done in Section 4 41 The derivations for in (5 32) 
and Ua = T„ + d;, as in (5 46) and (5 47) would still remain the same as in Section 
5 3 In case of specified second order input dynamics, for e g , another derivative of 
Tu would be required and the term would be passed on to the inner loops as was 
done in Sections 4 2 and 4 4 1 From its definition in (5 46), T„ would require another 
dernative of the desired model responses Thus the simple disturbance propagation 
method offered by the ML approach may be utilized to handle higher order input 
dynamics also 

5.4 Conclusions 

In this chapter, a robust control law design method, using a ML strategy, for a nonlinear 
system with specified input dynamics under unmatched uncertainties is developed The 
method does away with the requirement of higher order derivatives of the disturbances 
and uncertainties, leading to a simplified robustness analysis A discontinuous SMC 
law or its continuous approximation is designed in the inner-most loop so that the 
control inputs track a desired input signal generated in the outer loops Thus, the 
method results in continuous control inputs with discontinuities in its first or second 
derivatives only Since, it has been developed using a model following (MF) framework, 
higher order derivatives of the model (desired) responses are needed The method uses 
a full state feedback, including a measurement of the applied mputs The method 
is first developed for a SI nonlinear spring-mass-damper system It is shown that 
accelerometer signals are not necessary, although if present, they can be made use 
of Simulation results for the problem establish the applicability of the method The 
method is then applied to the MI aircraft problem with unmatched uncertainties as 
considered in Section 4 3 The bounds for the position and velocity errors are developed 
Results obtained from numerical simulation for different gam sets are presented and 
compared The ML strategy based SMC law synthesis method for a system with a 
specified input dynamics thus offers a alternate approach to the problem Also, the 
imolved robustness issues are much easily tackled in this framework as compared to 
the other methods 



Chapter 6 


Conclusions Sz Scope for Further 
Work 


The thesis is a study on some practical application aspects of the SMC design method 
The nonlinear longitudinal dynamics of a generic aircraft with the force-moment vector 
as the input has been treated as a benchmark problem The capability of generating 
a desired force-moment vector, given the airframe and engine characteristics, needs 
further study The control laws are developed under the assumption that an mverse 
map from the force space to the actual control surface deflections exists Although 
contiol laws have been developed for the longitudinal dynamics only, the method may 
easily be extended to the complete aircraft dynamics The main results of the study 
are as follows 

1 A new method of selection of parameters for the constant plus proportional rate 
reaching law based SMC based on the criterion of the time to a pre-specified level 
of output tracking error (OTE) is proposed It is established the method guar- 
antees a robust performance in the reaching phase also and gives faster tracking 
performances as compared to controllers based on the reaching time criterion. 

2 The SMC has been used to design controllers for the aircraft dynamics which 
IS nonlinear in the input This is achieved by using a force control approach 
(FCA), treating the force-moment vector as pseudo-inputs It is assumed that 
there are methods to generate any forces and moments as commanded by the 
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control law The error between the commanded and the generated forces is 
modelled as a bounded multiplicative uncertainty The conventional MIMO SMC 
design method is extended to include this uncertainty Analytical control laws 
are developed for the tracking and trimming modes for the longitudinal aircraft 
dynamics m the presence of matched uncertainties 

3 The multi-loop (ML) strategy has been utilized to propose a method of designing 
SMC laws m the presence of unmatched disturbances and uncertainties The 
position and velocity tracking errors are derived and shown to be uniformly ulti- 
mately bounded Analytical control laws are synthesized for the tracking mode 
and verified through numerical simulation using the parameters of the F-16 air- 
craft 

4 Using the better disturbance propagation property of the ML approach, a pro- 
cedure for the selection of the gain of the switching term is proposed which 
permits extension of the conventional SMC design methods to the unmatched 
uncertainty case The switching surface parameters can be selected using the 
developed bounds on the output tracking error 

5 The ML structure has been extended to generate SMC laws for systems with 
specified input dynamics in the presence of matched and unmatched uncertain- 
ties and disturbances Robustness analysis has been carried out and the ultimate 
boundedness of the position, velocity and input tracking errors are established 
Simulation results are obtained for a nonlinear spring mass system and the air- 
craft longitudinal dynamics and establish the efficacy of the proposed approach 

Thus, the present work establishes that the ML (multi-loop) strategy can be used 
to advantage m designing SMC laws for mechanical and in particular aircraft systems, 
considering some of the practical issues involved With the FCA as proposed in this 
study, and an estimate of the aircraft mass and inertia parameters, the developed 
control laws can be used to generate force-moment specifications to be met by the 
aerodynamics and propulsion groups at preliminary design stages of a development 
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project for the aircraft to have desired the HQs The explicit model following structure 
allows for easy alterations in the desired model responses at any design stage or even 
in flight 

The capabilities of the SMC method need to be extended for it to become a powerful 
alternative to the conventional PCS design methods For this, further studies along 
the following directions will be beneficial 

1 The work relies on the existence of a map from the force space to the effector 
space The capability of varying all the components of the force-moment vector 
independently needs further study It will be interesting to find out the minimum 
number of control eff'ectors required for this purpose and their design This issue 
needs further detailed study, before the FCA (force control approach) can actually 
be applied to the full envelope FCS design problem 

2 The control laws developed m this work use a full state feedback The case of 
output feedback and the design of nonlinear observers based on the SMC approach 
to generate estimates of the remaining states should be studied 

3 It IS desirable to extend the present study to include bounded control input and 
rate saturations 

4 The specification of model responses based on the handling quality (HQ) require- 
ments and its selection in the case of an existing airframe with a given limited 
force generation capability should be studied The advantage of a model follow- 
ing approach may be utilized by specifying different linear model responses for 
different operating conditions, which could depend not only on the input and rate 
saturations but also on the need to limit some particular states of the aircraft 
These model responses may then be scheduled instead of the controllers as m the 
gam-scheduled FCS design approach This could also be achieved by specifying 
the desired response through a single nonlinear system 

5 The control laws synthesized in this work are analog Most of the controller im- 
plementations, hovever, are now done on digital computers Thus it is desirable 
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to consider SMC design for sampled data feedback systems 

6 The results for second order systems (c f Chapter 2) would be of use for many 
mechanical systems, for e g , robot applications Thus, further study needs to be 
carried out to develop better estimates of the time to reach a specified level of 
output tracking error 
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Appendix A 


The S response 


The s dynamics is shown to be given in (2 18 ) as 

s = —es — i7'sgn[5], (A 1) 

where, 77' depends on disturbances and uncertainties and is thus an unknown function 
of the time variable t However, it is known to be bounded and 77' > 77, where 77 is 
a constant, for all t Let the response of the above system be indexed by the gain 
of the switching term, 1 e as s,,/ Here, the decay performance of s,,/ is studied In 
particular, consider the case of s decaying from an initial positive value of Sq at t = 0 
The above reaching law shows that s = 0 is an equilibrium point of the system and 
further that once s reaches zero, it will stay there for all future time instants Starting 
from a positive Sq, the s will always be non-negative The sgn[s] term is thus replaced 
by 4-1 in the above equation Thus, the s dynamics becomes 

s = —es — T]', or 

= —es - 77 4- (77 - 77' (t)) (A 2) 

The r}-rf{t) = n{t) term is seen to be non-positive The solution of the above equation 
may be written [ 69 ] as 

(t) = + n{T)]dT, or 
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< 5oe~‘* + ^ (A 3) 

since n() < 0 for alH For negative initial values of Sq, it can be similarly shown that 

Sr)'{t) > + f e~^^^~'^^[rj]dT, (A 4) 

J 0 

The right hand side of the above inequalities (A 3) and (A 4) represents Sr,, i e the 
initial condition response of 

s = —es - T^sgn[s] (A 5) 

These inequalities may be rewritten as a single inequality as 




(A 6) 


which implies that if r)'{t) > rj for all t, the decay of s exhibited by (A 1) is guaranteed 
to be better than the response of (A 5) as given in (A 6) Therefore, the reaching time, 
1 e the time taken to arrive at s = 0 starting from Sq will be bounded as 


e V ^ . 


(A 7) 


which is the reaching time for (A 5) [37] The (A 6) also implies that increasing the 
magnitude of 77 leads to a faster decay of s and lower reaching times 

If, the sgn[s] term is replaced by a continuous approximation as sat[s/<?i>] in (A 5), 
where (p represents an appropriately selected boundary layer thickness, then the s 
dynamics is given by s = -es - r}S&t{s/(l>) The time taken from an initial condition of 
So at t = 0 to reach the boundary layer may be evaluated as 



/elsol+^\ 

\ 6(1) + T] ) 


(AS) 



Appendix B 


Proof of Claim 1 


The following proof shows how a 5 > 0 may be selected, under the assumptions in the 
Claim 1 and starting from the premise that the (2 32) is true, i e , 

Q = §P + H (B 1) 


The P and Q matrices m the above equation are as given m (2 25) and (2 29) respec- 
tively 

Proof: 

From the P and Q matrices as given, the elements of the symmetric matrix H = 
Q - 5P from (B 1) are obtained as 


/ill 

= 2a^e- 6(0^ + J3) 


hi2 

= 2ae — ^ — 5a 


/l22 

= 2e — 5, 

(B 2) 


where, represents the element H{i,3) It is required that H be at least a posi- 
tive semi-definite matrnx The positive definiteness of this matrix would require that 
hn,h 22 > 0 and det[K] > 0, where det[] represents the determinant of its argument 


The two former inequalities gives 


6 < 2e 
6 < 2e 




a^ + 13' 


and 


(B3) 
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It IS seen that these inequalities are in agreement with one another For the last 
condition, the det[Fr] may be easily obtained from its elements as defined in (B 2) as 
hnh 22 — h \2 and after some algebraic simplification it may be written as 

det[ir] = /3[d2-2(a + e)(J + 4a€-/?] (B 4) 

The positive definiteness of P and Q imposes the condition that 0 < y? < 4Q;e as has 
been discussed in Section 2 3 Since /? > 0, the term factored out m the above equation 
IS dropped in the further developments The above equation thus shows that det[i?]//? 
is thus a quadratic function (<5/) m d as 

Qf{5) = 5'^ -2{a + e)5-\-iae~ (5, (B 5) 

and the constant term in the above equation is seen to be positive from the bounds to 
be satisfied by /3 as given in (2 30) The roots ri,r 2 of this quadratic equation may be 
easily calculated as 

’’1, r2 = {a + €)±^J{a-e)^ + P (B 6) 

Note that \J{a — e)^ + /3 can always be represented as (a — e) + di where di > 0 may 
be uniquely found Thus, the roots in (B 6) may be simplified to get 

ri,r2 = 2Q; + di,2e-di (B 7) 

Note that a small positive value for /? would imply a small di It can be shown that 
even for values of /? as large as permissible by its bounds, the r 2 = 2e — di > 0 Since 
the coefficient of the second order term in Qf in (B 5) is seen to be positive, it implies 
that the Q / has a minimum, and the corresponding value of 5 (= 5m) which gives the 
function minimum can be evaluated as 

d,7j = a + €, (B 8) 

which being the arithmetic mean of ri and r 2 will always lie between the two roots. 
Combined with the fact that Q/ = 0 at the roots 5 = n and ra (B 7), it is seen that 
the Qf < 0 IS possible only for values of 5 between the two roots and it would be 
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positive everywhere outside of this set Moreover, at <5 = 0 the value of Q/ is seen to 
be 4o;e — /?, which is positive as discussed earlier 

Thus it is evident that a positive 5 can always be selected from the sets given by 
(0, 2q;] or (0, 2e), whichever is smaller The two may be combinedly expressed by letting 
7 = min[ 2 o;, 2 e], where min[ , ] is taken to give the minimum of its arguments Thus, 
the permissible values of 5 may be specified by an open set as 

<^e(0,7) (B9) 

Moreover, in the limit of /? -> 0, a choice of J -> 7 would be achieved It is thus 
established that a positive <5 always exists which would result in H being a positive 
definite matrix, thus satisfying the Claim 1 ■ 



Appendix C 


Proof of Claim 2 


The arguments of the candidate Lyapunov function V{, ) in (2 40) are omitted when 
no confusion is likely to arise Note that this candidate function is not differentiable 
on the planes Sp = 0 and S', = 0 However, it is possible to use the generalized 
Lyapunov theorem, which allows non-differentiability of the Lyapunov function in a 
set of measure zero [41, 60] It is obvious that H is positive defimte function with 
7(0,0) = 0, and positive everywhere else in the e-e plane The at least negative 
semi-definiteness (n s d ) of 7 is established next for the two cases 

Proof* Using (2 31), The time derivative of (2 40) may be derived following the 
steps as outlined in Section 2 3 and gives the V here as 

7 = —y^Qy + "h ^ [^>SpSgn[iS'p] + 5',sgn[iS',]j , (C 1) 

except for the two planes where the function is not differentiable Note that q is 
different for the Cases CP and CQ as they have different switching surfaces S, and 
S, However the P is common for both the cases It has been selected as F m (2 25) 
for the Case CP Thus, the /Fg may be rewritten as S,K, where h, = -i/sgnlS.l 
and the subscupts t = p.g represents the cases CP and CQ respectively Moreover 
,t follows from q' < m that the maximum value of |h,l = Vm The S, and S, as 
required in the above equation and are evaluated as follows With h. as defined above 
the e = -ae- be + K from (2 20), for the case represented by the subscript and 
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a = pq,b = p + q from (2 8), the Sp = e + pe may be simplified as 

Sp = —ae — (b — p)e-\-hi 
— —pqe — qe + hi, or 

Sp = ~qSp + hi (C 2) 

Similarly, Sg can be obtained as 

S, = -pSg + hi (C3) 

Thus the time derivative of the candidate Lyapunov function as in (C 1) may be rewrit- 
ten for the case using the y'^Pq = Sph^ as defined above and Sp and Sg as in (C 2) 
and (C 3), as 

K = -y'^Qy + 2Sphi -f M [(-g5'p + /i,)sgn[5p] -H {-pSg + hi)sga[Sq]] , or 
= — -h 25phi - M[g|5p| TplS',!] + M[sgn[5'p] + sgn[5g]]/i, (C 4) 

It is seen from Fig C 1, that the two planes, 5p = 0 and Sg = 0, divide the e — e 
plane into four regions K as evaluated above is checked for sign definiteness in each 
of these regions for the two cases First the Case CP is considered 
Case CP 

In this case (i = p) and hp = -7?'sgn[S'p] Thus, it is obvious that 

hp > 0, when Sp < 0, and 

hp < 0, when Sp > 0 (C 5) 

Let Vp represent V for this case Note that, the —y'^Qy — +p|‘5'g|] term m Vi 

in (C 4) is always negative definite The sign definiteness of the remammg terms are 

established in the four regions 

1 Region 1 In this region 5p > 0 and Sq> ^ Thus, [sgn[5p] + sgn[5g]] — 2 and V 
m (C 4) reduces to 

Vp = -y'^Qy - M[q\Sp\Pp\Sg\] + 2Sphp -h 2Mhp 
With hp < 0 (C 5), it can be seen that < 0 in this region 


(C6) 
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Fig C.l The e-e plane 

2 Region 2 In this region Sp <0 and 5, < 0 and V becomes 

Vp = -y'^Qy - M[q\Sp\ + p|5,|] + 2Sphp - 2Mhp (C 7) 

In this region, hp > 0 which makes Vp <0 

3 Region 3 In this region Sp > 0 and 5, < 0 and with [sgnfS'p] + sgn[6'5]] = 0, y 
may be obtained as 

= -y^Qy - M[q\Sp\ +^15,1] + 2Sphp (C 8) 

From (C 5), /ip < 0 ensures that Vp is negative 

4 Region 4 In this region Sp <0 and S', > 0 and V becomes 

Vp = -y'^Qy - M[q\Sp\ +p15,1] + 2Sphp (C 9) 

hp > 0 in this region, and thus the time derivative of V can be seen to be negative. 

The above checks establish that V in (2 40) is a Lyapunov function for Case CP Note 
that the Sphp term m the above equations may be rewritten as -T]'\Sp\ which is always 
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negative The Vp m the four regions as in (C 6) - (C 9) may be jointly expressed by 
the following inequality 

Vp < -y^Qy - M[j|5p| + p|S,|] - (C 10) 

It IS evident that V in (2 40) remains positive definite with the Vp above being negative 
definite for all M > 0, and thus the V is seen to be a valid Lyapunov function for the 
CP case 

The candidate function is next checked for the Case CQ 
Case CQ- 

In this case with t = q, hq is given as — 77'sgn[5'g], and thus, 

hq > 0, when Sq < 0, and 

hq < 0, when Sq> 0 (C 11) 

Let Vq represent V for this case 

1 Region 1 In this region S'p > 0 and > 0 V in (C 4) reduces to 

Vq = —y^Qy — M[q\Sp\ +p|5'g|] + 2Sphg + 2Mhq (C 12) 

With hq<0 from (C 11), it can be seen that Vq<Qm this region 

2 Region 2 In this region Sp < 0 and Sq<0 and V becomes 

Vq = —y'^Qy — P^[<l\Sp\ d-pl'S'^O + 25p/ig — 2Mhq (C 13) 

In this region, hq> 0 which makes Vq < 0 

3 Region 3 In this region Sp> 0 and Sg <0 and V becomes 

V, = -/Os - sMISpI - pM\S,\ + 2Sph, (C 14) 

Note that the first three terms are negative, however the last term is positive as 
hq> 0 (C 11) The above equation may be rewritten as 

Vq = —y^Qy — pM\Sq\ — {qhi — ^hq)\Sp\ (C 15) 
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It IS known that 77 ' and therefore |/i,| is bounded Thus if {qM - 2h,) > 0 , then 
Vq becomes negative in this region 

4 Region 4 In this region 5'p < 0 and 5, > 0 and F becomes 

y, = -y'^Qv - qM\S^\ -pM\Sq\ + 2S^hq (C 16) 

As in the Region 3, here also the first three terms are negative, but with /i, < 0 
(C 11 ), the last term is positive Rewriting the above equation as 

K = -y^Qy - vM\Sq\ - {qM + 2/15)15),) (c 17) 

If + 2 /I 5 ) > 0 , then the time derivative of 17 can be seen to be negative Here, 
hq IS negative and hence the inequality is the same as m the previous region 


Thus it IS seen that, if 

q ~ q ' 


(C 18) 


then V (2 40) becomes a Lyapunov function for Case CQ, since \hq\ < rju In partic- 
ular , if M IS selected as 

M=^. (C19) 

where 0 < 5q < q, then qM — 2hg term in (C 15) and (C 17) may be replaced by SgM 
The Vq m (C 12 ) - (C 17) for the four regions may then be rewritten jointly as 


Vq < -y'^Qy - pM\Sq\ - SqM\Sp\ (C 20 ) 


This proves the claim 



Appendix D 


The Aircraft Parameters, 
Disturbances, and Uncertainties 


In this work, the model responses are generated by 


AA = + BmAUp, 


(D 1) 


where A A = A — Xeq, where A = [Vt, a, Q, dY represents the vector of desired 
state trajectories and AC/p = Up — U^q, where Up represents the pilot mputs, namely 
the elevator stick and throttle gear deflections The X^q and Ueq denotes the quantities 
at the equilibrium condition Note that the aircraft states in the chapters are taken 
as [U, W, Q, 6Y Thus, appropriate transformations for the first two states needs to 
be performed The remaining states (namely, Px and p^) not specified in the above 
equation evolve according to (31) The matrices Am and Bm were taken as 


A 


m 


2moE + 0 

- 30237T; - 3 

- 24369T^ - 2 
000005; + 0 


482545; + 3 

- 11753E + 1 

- 126635; + 2 
OOOOOE + 0 


- 350395; + 2 
854885; + 0 

- 588055; + ! 
100005; + 1 


- 609605; + 3 

- 815185; - 2 

- 456275' + 1 
000005 + 0 
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402395 + 0 170395 + 0 

_ - 315705 - 4 - 215075 - 2 

tSm — 

000005 + 0 - 182435 + 0 

000005 + 0 000005 + 0 

The above model gives desired trajectories in the Vt- a terms where as these are 
required to be given m the 5 - frame in the chapters The quantities may be 
transferred from 1^ — a reference frame to the U — W using the following formulae [ 10 ] 

U = T^cos(a), and 

W = V^sin(a), (D 3 ) 

or inversely as 

Vt = U^ + W^, and 

tan(a) = W/U (D 4) 

The Ud and Wd required in the work are evaluated by taking the derivative of D 4 
which on simplification using (D 3) gives 

UdUd + WdWd = Vt,Vt,, and 
-WdUd + UdWd = Vtlad (D5) 



The above equation may be solved for Ud and Wd, in terms of the known Vtj, ad, and 
therefore Ud, Wd and the time derivatives Vt^ and ad to get 


Ud = 

= ^^- 7 - - Wdad, and 
n 


Wd = 

a. rr ^ 

= T/ UdOid 

n 

(D 6 ) 


The second derivatives of the model response are developed by taking a time deriva- 
tive of (D 1 ) and simplifying it to get 


AX = AmAX + BmAUp, or 
AX = Al,AX + AAAU + BmAUp 


(D7) 
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Thus, the second derivatives Vt^ and are easily calculated Note that they contain 
the deiivative of the applied pilot inputs It now remains to transform these into Ua 
and ll'd as are required in (5 43) For this, taking another time derivative of (D 5) and 
rearranging terms, it is seen that 

UiUd + WaWi - and 

-WdUd + UdWd = 2Vt^Vt^ad + Vt^ad (D 8 ) 

The above linear equations may be solved for the unknowns Ud and Wd to get 

= -^Pd{Vi^ + VtiVt^~Ul-Wl)-Wd{2Vt^Vt^ad + Vlctd)], 

= ^mVi+V„V,,-U^-WS) + Ud2V„V„o„ + V=a,)] (D 9) 

The remaining second order derivative in (5 43, i e the Qd is directly available from 
(D 7) 

The uncertainties and disturbances considered considered in this work are 

• Uncertainty in the Mass and Inertia: The nominal values of reciprocal of 
aircraft mass, = 1 5 x 10 “^, and inverse of the pitching moment of inertia, 

C 7 = 1 75 X 10~® have been used The range of their possible variations have been 
taken as 1 0 x 10 "^ < Ki < 2 0 x 10 “^ and 1 0 x 10 "® < C 7 < 2 5 x 10 "® 
The actual values of these parameters as used in the simulation are reported in 
Fig D 1 The bounds Dbi and Z ?52 can be evaluated as 1/3 and 3/7 respectively 

• Uncertainty in Aerodynamic Inversion and/or Modeling; It is assumed 
that a maximum 20 percentage of error can occur, which means that \5ae\ < 0 20 
For simulation purposes two different sets of variation m iJoe has been used as 

UCl Uncertainty set 1 UC2 Uncertainty set 2 

5ae(l) = 0 05 + 0 145 COs(37rt) (iae(l) = -0 19, P 

<^ae(2) = — 0 02 + 0 15cos(57n{), and <iae(2) = 0 175, and 

^ae(3) = 0 1 + 0 08 COs(37rt) ^ae(3) = 0 185 

for the 3 -component force vector 
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1 95 


C7 X 10® 


1 85 


1 75 


Fig D 1 Variation m the Rm and the cr parameters of the aircraft 

• External Disturbances in the governing differential equation For simu- 
lation purposes the external disturbances were assumed to be 

djj = 0 05 -I- 0 25 cos(37rt), 

= -0 05 -hO 25cos(77rt), and 
dq = 0 05 -t- 0 15 cos(57rt) 
dg = 0 04 -1- 0 25 cos(37rt), 
da; = -0 74^0 + 0 25 cos(77rt), and 
4 = 0 05d0 + 0 85 cos(57rt) (D 11) 

with Du = 0^, Dw = 0 3, and Dq = 03 Note that the other disturbance terms 
are unmatched Their bounds are taken as Dg = 0 3, Dx 10, and Dz 10 

The boundary layer thicknesses have been selected as (^ = [0 5, 0 5, 0 005] for the 
tracking mode and = [0 5, 0 01, 0 01] for the trimmmg mode 




Appendix E 


Time Derivatives of Functions 


The definitions and notations of time derivatives of smooth scalar, vector and matrix 
functions of a vector of variables used in this work is given here The vector of variables 
IS taken as x = [ 2 : 1 , 0 : 2 , , XnY 6 and x is assumed to be available Also arguments 

of the functions are omitted where they are not likely to create confusion 


Derivative of a Scalar Function 


Let /(x) be a scalar function Then its time derivative is given usmg Lie derivatives [21, 
22] as 

rif Hf Hf 

L^f (E.l) 


df df df 

~ j 1 ~ Q ■*" a 

dt uXi C/Xfi 


In this work, the above Lie derivative will be expressed as 

df 


dt 


Lxf — -SyX 


(E2) 


R f represents a 1 x n element row vector formed from / as defined above 
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Derivative of a Vector Function 

Let f(x) [/i, / 2 , ) /m] be a 77T. X 1 vector function Then using the above notation, 

its derivative is given as 


i?/jX 


= MfX, 


where M{ may be rewritten using (E 2) and the above equation 


dfi 

dfi 

dfi 

dxi 

dh 

dX2 

dxji 

af2 

dxi 

dx2 

dxji 

dfm 

dim 

dfm 

dxi 

dx2 

dXn 


(E4) 


Mf represents a m x n matrix formed from the components of the vector function f as 
defined above 


Derivative of a Matrix Function 

Let ^(x) represent a m x p matrix function By definition, 

d[A(x)] _ dotj(x) 
dt dt ’ 


(E 5) 


wheie, aij(x) represents the element in the row and column of A Since Cy is a 
scalar function, its time derivative may be written using (E 2) as 

^ = i?„..x (E6) 
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The derivative of A (E 5), using the above equation can be expressed in a matrix 
product form as 


-Ran -Ratp X 0 0 

d[A] Rail Raip 0x0 
dt 0 0 

. Ramp J L 0 Ox 


(E7) 


where 0 in the above equation represent a n column zero vector Using the definition 
of a Kronecker product, and defining the first matrix in the above product as Ja, the 
above equation may be rewritten as 

(E8) 

at 


where, Jp denotes a p x p identity matrix, igi represents the Kronecker product In the 
Ja as defined in (E 7) it is clear that, since Ra,^ represents a 1 x n row vector, the 
matrix will be a m x np matrix Also, the [Jp ® x] matrix will be a np x p matrix by 
definition Thus the time derivative of A results m a m x p matrix The Ja matrix as 
defined above is seen from (E 7) to contain the partial derivatives of with respect 


to X as 


dan 

dan 

dxi 

dX2 

da2i 

da2i 

dxi 

dX2 

dajni 

daffii 

dxi 

dX2 


^Qii ^Qi2 dai2 daip 

dxn dxi ~d^ 

OCL21 Oa22 

dXn dxi 

dam! 90'm2 

dXji dXi ^^2 


(E 9) 



Appendix F 


The Uniform Ultimate 
Boundedness of Systems 


Definition 1 (Global uniform ultimate boundedness) [52] 

An uncertain feedback system with state x(t) is said to he globally uniformly ulti- 
mately bounded with respect to a set S C 'RA, if 

1 for each uncertainty realization and for each (to»x(to)) e'JZx'RA there exists at 
least one solution x( ) [to, ti) RA, ti > to, 

2 given any real number S > 0, there exists a real number d{5) > 0 such that, for 
any solution x( ) with ||x(to)|| < 5, ||x(t)|| < d{5) for all t 6 [to,ti) all solutions 
can thus he continued over [to,co], 

3 for every x(to) € RA, there exists a non-negative constant T{x{to),S) € RA such 
that, for every solution x( ) with to arbitrary, x(t) e S for all t > to+T(x(to), 5') 

Using the above definition, the global uniform ultimate boundedness of the followmg 
system is established The dynamics in (4 12) may be written m general form as 

X = —kx + d, (F 1) 
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where x 6 A: > 0, and d represents the unknown uncertainty term which is given to 
be bounded with (d| < D The subsequent developments follow the steps as outlined 
in [52, 70] Select the candidate Lyapunov function as 

= 5^’ (F 2) 

Its time derivative may be calculated as V = -kx^ + xd With the bound on d as 
specified above, it is seen that 


V < —kx"^ + D\x\ (F.3) 

Note that |a| = y/W from (F 2) Thus, V < -2kV + DV^ It is evident from this 
equation that V < 0 if V > ri—e, where e > 0 is a small constant, and ri = e+D'^/{2l^) 
Define the set S = x ^x^ < ri From the fact that V < 0 if F > ti - e, it is seen that 
if x(to) £ S then x{t) € S for all t >to The (F 3) may be rewritten as 

V < -2kV + DVW (F.4) 


The time taken for V to change from Vo at t = to to at t = ti may be evaluated as 


- to < ^ In 


k^/W^-D 

k./Wx-D 


(F 5) 


Thus, if x(to) ^ 5', then the trajectory will reach 5 m a finite time given by 




k^x^-D 
k\fS\- D 


All conditions for uniform ultimate boundedness are thus satisfied 


(F.6) 



Appendix G 


The Error Bounds 


The IS shown to satisfy the bounds as in (5 51) Putting = [|Bij?i>„+Di +I?t(|fi - 
TflDjj, it may be rewritten in a simpler form component-wise as 

ajj + D2U 


l-S'crl < 


Ku[l — Dt{l, 1 )] 


|5q|< 


aw + D2W 

Kw[^ — Dt{2, 2)] 

+ ^2Q 

Aq[1-A(3,3)] 


(G 1) 


Similarly, the IIS'!, I ^ire given to be bounded as in (5 52) and using the bounds on the 
A matrix in (4 29) , these may be expressed as 


Kx\Sx\ < Sxj -{■ Sw + D2X1 
i^z\Sz\ < Sxj -{■ Sw D2Z1 3 'D-d 

Kg\Se\ < Sq + D 2 B (^- 2 ) 

The bounds Dj as required in (G 1) are given in (4 43) and are seen to involve the 
bounds on the position errors, le ISe] Using the bounds on the position errors as 
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above in (4 43) gives 


K^D2^\cos{9)\ + K,D2,\stn{6)\ 

-‘rD 2 e [\UdSin{Se) \ + \ WdCos{Sg) \ + D 2 x + D 2 z] 

+2A4l'^a| + 

Kj;D2x\stn{9)\ + KzD2z\sin{9)\ 

+D 29 ]\UdCos{Sg)\ + \ WdS%n{Sg) \ + D 2 X + D 2 x] 
-\- 2 D 2 g[\Su\ + I'S'iyl] 

KgD 2 g 

I3u + 2D29[\Su\ + \Sw\] ^ 

/3vv + 2 D 2 g[\Su\ + |*5'iv|] ) 

KgDoa I 


(G3) 


where, /?[/, Pw are positive quantities and are obviously defined It is seen from the 
above equation that does not depend on any of the velocity error bounds Thus 
the bounds and the 15^1 may be directly evaluated as defined in (G 1) and (G 2) 
respectively However, the |5[/| and the bounds are implicitly defined and may 
be evaluated as follows 

Using the above developed bounds for D 2 (G 3), and putting a' — ctj+y^,, t = W, W, 
qi = Ku[l - A(l, 1)], and 52 = Kw[l — A(2, 2)], the IS'nl and 15^1/1 in (G 1) may be 


rewritten as 


\Su\ < 


a'xj + 2D2g[\Su\ + liS'v 


I'S'wl £ 


a'^ + 2A^[|‘S'i7| + 1‘^wl, 


(G4) 


The above inequalities may be expressed as set of simultaneous hnear mequalities in 
matrix form as 


1 — 2D2b/q.i ~2D2g/qi I ^ / Oi'ulqi 

—2D2e!Q2 l~2I?2e/92 y I'S'wl j 

Note that a', i — U,W and qi,q 2 are always positive The solution l^trl and [S'wl of 
the above equation should be positive to be practically meanmgful This may be easily 
established using the following theorem 
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Theorem 1 (Forbenius-Perron) [21] Consider a square matnx A with nonnegative 
elements Then the largest eigenvalue pi of A is nonnegative Furthermore, consider 
the equation 

(/ - p~^A)j = z, (G 6) 

where, y, z are vectors, I is the identity matrix of appropriate dimensions and all 
components of the vector z are nonnegative 

If p > pi, then the above equation admits a unique solution y, whose components 
Pi are all nonnegative 

From (G 5), it is seen that the z = [ay-/gi has all its components positive, 

and the A matrix, as m the above theorem, is seen from (G 5 to be 

2 D 29 I qi 2Z?2i9/9i 
2D20/92 ‘^ T ) 2 Blq 2 

The eigenvalues of the A matrix are easily obtained as 0 and 2D2b / (91 + 92) According 
to the above theorem, if the largest eigenvalue is less than 1 0, then the (G 5) is 
guaranteed to have a positive solution for the norm of the error bounds Now, by 
definition 91 + 92 = Ku[^ ~ A(l,l)] + Kw\}- ~ A(2, 2)], and since D 2 B is a known 
bounded quantity, sufficiently large positive values for Ku and Kw can be selected to 
make the non-zero eigenvalue less than 1 0 Thus, the error bounds may be obtained 
by solving the simultaneous linear equation set in (G 5) Further, once the |S„| are 
guaranteed to be positive, the positiveness of |Sa:| is obvious from its definition m 
(4 45) Similar results are presented in [70] using a different approach 



(G7) 



APPENDIX 1 
ERRATA 


The following corrections are to be included in the thesis 

1) Ref: Pg-6: 5th Line- 

" lining this result, is proposed' should be corrected to 

This result has been used to propose a new method of selection of controller parameters^ 
which gives a lower time to OTE 

2) Ref: Pg-9: 15th Lme- 

" Application of the unmatched uncertainty” should be corrected to 

The ML approach is then applied to a third order nonlinear system with unmatched 
uncertainties and disturbances A new method of selecting the gam of the switching term is 
proposed Further it is shown that the conventional SMC may be extended to tackle the 
unmatched uncertainty case with the proposed gain selection 

3) Ref: Pg-72: 18th Lme- 

"The Aae in the simulation” The 'Aae ' in this line should be corrected to 'A,,c' 

4) Ref: Pg-94: 14th Line- 

”As seen earlier in this chapter” should be corrected to 

As discussed in Section 4 2, the effect of disturbances in the outer loops is easily passed on 
to the inner loop This simple propagation of disturbances results in a simplified robustness 
analysis, as will be shown in this chapter 

5) Ref: Pg-94: 21st Line- 

” Following the developments in the innermost loop” should be corrected to 

Following the developments in Chapter 4, the tracking of position variables is done in the 
outer loop by generating augmented velocity command signals These augmented velocity 
signals are tracked in the inner loop by synthesizing augmented input command signals A 
SMC law is then designed to track the commanded input signals in the innermost loop 

6) Ref: Pg-H6: 10th Line- 

”In this chapter, is developed' should be corrected to 

In this chapter, the ML strategy has been used to develop a robust control law design 
method for a nonlinear system with specified input dynamics under unmatched uncertainties 
and disturbances 
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